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ABSTRACT 


A new method of stabilizing laser gas discharges through the 
use of a specially profiled magnetic field is described. The magneto- 
gasdynamic (MGD) stabilization technique utilizes crossed electric 
and magnetic fields to create and maintain a rapidly rotating plasma 
of large volume. This rapid recirculation creates a homogeneous normal 
glow discharge which suppresses the bulk thermal instabilities. In 
addition, large secondary flows were Eadie exist which serve to 
convect the hot neutral gas out of the active discharge volume. Sig- 
nificant increases in specific discharge power loadings have been 
achieved without the benefit of external gas convection. The gas dis- 
charge performance was found to be largely insensitive to the gas 
composition and pressure. 

Theoretical investigations of this new technique were carried 
out in order to predict the complete flow pattern and the gas discharge 
characteristics. This analysis was based upon a one-fluid MGD model. 
A computer code, MAGIC, has been developed to numerically solve this 
system of time-dependent, nonlinear, coupled partial differential 
equations. The algorithm used to solve the compressible MGD equations 
is based upon an iterative, alternating direction implicit (ADI) 
temporal advancement scheme. This method avoids the limitation in 
time step size imposed by explicit methods thereby substantially 
reducing computational cost. 

Numerical results from the computer code revealed a high speed 


rotational flow dependent upon the applied Lorentz driving force. 
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Furthermore, the secondary flows formed a convective cell for the 
neutral gas. Finally, the one-fluid MGD model indicated that the 
spatially nonuniform electric and magnetic fields within the gas 
discharge caused the plasma and neutral gas velocities to be strongly 
sheared within the cathode fall and positive column regions. This 
Sheared rotation serves to minimize plasma micro-instabilities near 
the cathode surface. These results were confirmed by direct obser- 


vation and by detailed diagnostic measurements. 
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determinant of the metrix tensor 


degeneracy factor of the Jth rotational level 


lineshapping function 
Planck's constant 


scale factor for the x coordinate 


Hartmann number 


initial laser intensity 


total discharge current 
unit tensor 

current density 
Jacobian tensor 
Boltzmann's constant 


collisional rate coefficient 
gas permeation constant 


iteration index 
inductance 


scale length 
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neutral, electron, and ion mass respectively 
average molecular mass of a gas mixture 


Mach number 


Molecular weight 

neutral, electron, and ion particle density respectively 
negative and positive ion particle density respectively 
Avogadro's number 

number of degrees of freedom of the zth gas 

population density of the /th level 

neutral, electron, and ion kinetic pressure respectively 


power of a convection and diffusion cooled laser 
respectively 


static and stagnation gas pressure respectively 


polyvinyl chloride 
gas throughput 


collisional cross section between the ry and s gas species 


radius 


Larmor radius 


universal gas constant 


plasma impedance 
magnetic Reynold's number 


radio frequency 
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pumping speed 


ionization source term of the jth process 
collision integrals 


time 

spontaneous transition time 

neutral, electron, and ion temperature respectively 
transversely excited atmospheric 

mobi lity 

solution state vector 


ultra-violet 


Alfvén speed 
mass average velocity of the s species 


voltage or volume 
voltage-current 


electron, negative-ion, and positive-ion drift velocity. 


particle velocity vector 


excitation probability of zth state 
axial direction coordinate 


ionization number for the s species 


first Townsend ionization coefficient 
absorption coefficient 


electron and ion Hall parameters 


ratio of specific heat capacities 
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y(v) 


optical gain coefficient 
particle flux density 


energy loss factor for the s species 
Kronecker delta function 


incremental change 
convergence criterion 


free space permittivity 
Levi-Civita permutation symbol 


bulk viscosity coefficient 
azimuthal direction cordinate 


conservative molecular property of the s species 
thermal conductivity 


characteristic wavelength 


mean free path 
Debye length 


Shear viscosity coefficients for charged and neutral 
fluids respectively 


frequency of radiation 
electron-electron energy exchange collision frequency 


average collision frequency for momentum exchange between 
ry and s particles 


electron energy 
relaxation energy 
partial pressure tensor of the s species 


neutral, electron, and ion mass density respectively 
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electrolyte fluid condictivity 


electrical conductivity tensor 


parallel, perpendicular, and transverse components of 
the electrical conductivity tensor with respect to the 


magnetic field 


characteristic convective gas cooling time 
characteristic diffusion time 

characteristic time for the jth instability process 
discharge rotation time 


viscous stress tensor 
angle of rotation between coordinate systems 


plasma frequency 
electron and ion cyclotron frequencies respectively 


vorticity vector 


. collision frequency tensor 


perpendicular 

parallel 

transverse 

infinity 

gradient vector 

denotes a tensor quantity 


denotes an averaged quantity 
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CHAPTER 1 


INTRODUCTION 
1.1 PROJECT MOTIVATION 


The study of naturally occurring and laboratory created plasmas 
presents an enormous challenge to our fundamental knowledge of this 
fourth State of matter. Various successful investigations have led to 
the modern day development of high voltage thyratron switching devices, 
magnetohydrodynamic (MHD) power generating systems, and many other 
discharge devices. In addition, the analysis of ionized gases has 
prepared a background for theoretical fusion studies. 

In recent years, gas discharges have also provided an efficient 
pumping mechanism for a new generation of high powered gas lasers. 
Experience has shown that the performance of these electric discharge 
lasers is highly dependent upon the quality of the glow produced within 
the laser. The demand for larger volumes, higher pressures and complex 
gas mixtures has placed severe limitations on present day aas dis- 
charge technology. Specifically, an initially uniformly distributed 
discharge is observed to undergo a glow-to-arc transition whenever the 
input power density exceeds a certain instability threshold level. 
Numerous other plasma instabilities can also produce local current 
constriction which in turn results in a very nonuniform excitation and 
excessive local gas heating. Due to these circumstances, the scaling 
of larger devices to higher energy densities has become increasingly 


difficult. 
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In order to circumvent these laser discharge stability problems, 
various new stabilization techniques have been developed. Most of the 
attention has been focused upon the basic principle of convective gas 
cooling and its effect upon a glow, sustained by an auxiliary ioniza- 
tion source. Devices employing both these techniques simultaneously, 
have proved to be most successful thus far. However, these procedures 
have still been unable to overcome the scaling problems associated with 
the development of large volume, high power lasers. 

The major objective of this research endeavor is to continue the 
study of laser discharge instability processes. In an attempt to 
eliminate the effects of these undesired processes, a new magnetogas- 
dynamic (MGD) stabilization approach has been devised. 

In this chapter, a brief history of high power lasers is presented. 
Included are some of the fundamental physical processes which determine 
useful laser operation. Some initial concepts, concerning a plasma 
under the influence of a magnetic field, are also discussed so as to 


provide the necessary background material for this project. 


1.2 BACKGROUND HISTORY 


Early designs, of electrically excited molecular lasers, con- 
sisted of either sealed off tubes or long tubes containing slowly 
flowing gas mixtures [1.1]. These gas discharges were cooled by 
molecular diffusion and consequently, vibrational energy was lost to 
the tube walls. For this type of laser discharge, the energy balance 


or resulting power density can be approximately expressed as 
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Here, : is the gas discharge current density vector, E is the applied 
axial electric field, and 7 is the resulting gas temperature. This 
equation clearly shows the thermal constraints imposed upon the gas 
discharge. Since He is the major constituent in C0. lasers, the thermal 
conductivity, Kn? of the gas mixture is effectively constant. As a 
result, higher power densities could only be attained by maintaining 
very cold walls and by limiting the tube radius to about 1 cm. 

In these diffusion cooled lasers, the maximum power level is 
determined by the rate at which waste heat is removed from the system. 
For a discharge column of diameter d, the characteristic diffusion 


time is 


where os is the mean free path of the C0. molecules and e is the 
thermal molecular speed. Since the achievable power is inversely pro- 
portional to the characteristic cooling time, it can be shown that 
hae 
Fave eS 
where o is the gas density. Consequently, this axial flow design could 
only produce power densities of a few tenths of a Watt/cm?; thus 
forcing the construction of very long discharge tubes. Any further 
increases in input power led to higher gas temperatures and a sub- 
stantial reduction in gain [1.2, 1.3]. 
Recognition of the important role played by gas temperature and 


transport processes led to the development of a convective laser. In 
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this design, excess heat is removed convectively by the rapid flow of 
the gas mixture. If the gas is flowing at a speed v, then the character- 


istic cooling time is given by 
i i ayo. . 


The flow of fresh material into the active region provides more efficient 
laser pumping and reduces the thermal population of the lower laser level. 
In order to describe the factors which promote higher power levels 
and system scaling, the translational and vibrational energy equations 
were examined. The resulting studies revealed an energy dependence upon 


gas density as well as flow velocity. In other words, 


A comparison of achievable power levels between a diffusion cooled 
laser and a convectively cooled laser, with equal active volumes, was 
simply proportional to the ratio of their characteristic cooling times, 


i.e. 


Since the ratio was much smaller than unity, it was obvious that con- 
vective cooling was a superior design for a laser. 

Although the utilization of convective cooling increased laser 
discharge power density by a factor of ten, most high power molecular 
lasers could only operate at levels substantially below those dictated 
re thermal considerations. This was primarily due to the onset of 


plasma instabilities. 
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As closed-cycle lasers became commonplace, it was found that 
instabilities limited the discharge power to levels significantly below 
those characteristic of open-cycle lasers. A new technique involving 
the application of RF power to the DC discharge was then introduced 
[1.4]. It was observed that the radio frequency source improved 
discharge uniformity and hence stability. In addition, the enhanced 
uniformity kept the current density to a minimum, which in turn, 
minimized the local buildup of arc-like filaments. The RF supply was 
also thought to produce a rotation in the current density vector. This 
produced a damping of unstable wave modes, provided the RF supply fre- 
quency exceeded the growth rate of the unstable wave. This process 
increased the discharge power density to some extent. However, as the 
RF power increased further, the applied DC power peaked and then 
decreased due to thermal loading effects in the gas. Hence, this pro- 
cedure has only limited capability and use. 

Further improvement in gas discharge stability was obtained 
through the use of external control of the electron density. This led 
to the development of the electron beam sustained laser discharge [1.5]. 
In this technique, uniform eeatetuen is created throughout the 
active volume by a powerful electron beam. The external ionizing agent 
decouples charge production in the discharge from the electric field. 
This decoupling reduces a primary source of instability caused by the 
loss of electrons due to volume recombination of electron-ion pairs. 
Though effective, the electron beam controlled laser has some limita- 
tions with regards to gas discharge current, pulse duration, and the 


frequent occurrence of foil damage by energetic electrons. 
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Another common technique for conditioning the discharge is the 
use Of photoionization [1.6-1.9] in the form of uv light sources. In 
this configuration, the high energy electron beam is replaced by an 
inexpensive low energy photon preionization source. Here, the applied 
electric field no longer sustains the plasma, and the ratio Z/n can be 
independently adjusted in order to optimize vibrational excitation of 
the upper laser levels. This method can also provide a spatially uni- 
form background for avalanche development in the main discharge. How- 
ever, if the preionization electron density is less than the equilibrium 
electron density produced by the avalanche discharge, the preionizer 
loses control after the avalanche develops. The stability of the main 
discharge plasma is then again restricted by the occurrence of a glow- 
to-arc transition. 

A unique technique that has had some success deals with aero- 
dynamic gas conditioning [1.10, 1.11]. The introduction of turbulence 
upstream from the gas discharge has a noticeable stabilizing effect at 
Tow pressures. Turbulence affects the plasma in two important ways. 
First, it can change microscopic properties such as the particle 
velocity, energy, and density distributions. Secondly, macroscopic 
features can be altered such as the overall plasma shape. Experiments 
have shown that a weakly ionized plasma cannot significantly alter the 
turbulent field of the neutral gas, although jion-neutral collisions 
readily transfer turbulent motions from the uncharged to the charged 
species [1.11]. As a result, these small-scale eddies of turbulence 
enhance neutral particle and plasma transport processes. Such enhanced 


transport properties cause the critical wavelength of most instability 
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modes to increase. This allows more time for the gas flow to convect 
any current filaments out of the discharge region. 

Associated with this turbulence is an increase in particle loss 
and a rise in electron temperature. Under certain conditions this can 
favourably alter the conditions for an instability occurring in the 
negative-ion and ionization modes. However, the observed improvements 
in discharge stability have been difficult to interpret because of a 
lack of detailed measurements and by complications created by other 
possible effects. 

Even with the successes achieved through the use of such inno- 
vative techniques, problems still remain which prevent the development 
of even larger volume, higher energy lasers. In particular, maintaining 
a stable, uniform plasma of high optical quality with dependable long 
term qualities has been difficult to achieve. Consequently, the need 


for new dynamic, uncomplicated methods of stabilizing discharges exists. 


1.3. POSITIONAL STABILIZATION OF A PLASMA COLUMN 


Arc positional control and extinction techniques by means of 
magnetic fields have been considered for several decades [1.12, 1.13]. 
These typical interactions considered between arc, gas flows, and 
magnetic fields cry dealt with two situations: those in which 
the arc was moved through a gas by an externally applied magnetic field, 
and those in which the arc was held stationary by a balance of forces. 
Such balanced arcs involve a match between aerodynamic and magnetic 


forces. Initial experiments determined that the magnetic field 
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strength necessary to balance the arc column was proportional to the 
Square of the applied transverse gas velocity. 

Such studies have led to a technique of magnetically stabilizing 
the position of a low density plasma column. In early low pressure 
convectively cooled axial discharge C0. lasers, the plasma column was 
strongly bowed downstream due to the influence of the transverse gas 
flow. For efficient laser operation; the discharge has to be straight- 
ened and aligned with the optical axis of the laser resonator. The work 
done by Buczek et. al. [1.14], involved the use of a tapered transverse 
magnetic field mutually perpendicular to the axial electric field and 
the gas flow velocity. Positional stabilization of the discharge column 
was accomplished by the balance between gas flow velocity and the drift 
velocity of the plasma caused by the resulting Lorentz force. The 
interaction of the transverse magnetic field and the charged particle 
drift velocity Vg» created a force which counters the flow velocity 
force. Thus, the position of the discharge column could be accurately 
placed. 

A tapered magnetic field was designed so that it would force the 
discharge to flow on the desired axis. Essentially, this created a 
potential well in which the discharge was held against velocity and 
current fluctuations. The size and depth of this well was dependent 
upon the shape and strength of the transverse magnetic field, and upon 
collisions with flowing neutral molecules. Thus, using a properly 
designed magnetic field, the plasma column could be confined. By 
aligning the potential well with the laser's optical axis, most of the 
electric energy deposited in the plasma column was available for laser 


output. 
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These and other studied have indicated that for low density gas 
discharges, the magnetic field necessary for positional stabilization 
varied linearly with pressure and flow velocity. These important facts 
allow the idea of system scaling to be considered. Thus, it should be 
possible to scale up the power density by increasing gas flow velocity 
and magnetic field proportionally, and still maintain a high operational 
efficiency. 

Present day convective lasers utilize a transverse electrode 
design rather than the cross-field electrode design of Buczek. In 
these new laser devices, the gas discharge is transverse rather than 
parallel to the optic axis. Since the distance between electrodes is 
now much shorter in the transverse design, the problem of the discharge 


bowing downstream is significantly reduced. 


1.4 ROTATING PLASMAS 


Recent investigations of magnetized fluids have explored pos- 
sible applications of rotating plasmas [1.15, 1.16]. Various studies 
have produced a wealth of information dealing with plasma confinement 
and momentum balance of such rotating plasmas. Several research lab- 
oratories are studying the strong centrifugal force created by the 
rotation and are exploring its possible uses. They have found that 
this force can be externally controlled and used to simulate the effect 
of a large thermal pressure. It also seems likely that special rotating 
plasma devices could be built to predict various transport coefficients 


in a magnetic field. Rotating plasmas are also important in fusion 
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research since the centrifugal confinement should help reduce the end 
losses in magnetic mirror devices. 

Ionized particles, within a gas discharge, begin to gyrate under 
- the influence of an electrodynamic force J x B, produced when an 
electric current flows through a magnetic field. This force produces 
a torque only when the current density and magnetic field components 
are mutually orthogonal. Thus, macroscopic fluid motion is created. 
Such motions may also arise from the interaction of time-dependent 
electric and magnetic fields. 

The most significant aspect of a rotating plasma is the ability 
of the J x B force to produce a strong transmission of momentum in a 
direction orthogonal to the magnetic field, thereby creating a Lorentz 
force. In addition, this large force has proven to be an effective 
. method for removing impurities from within the plasma bulk. This idea 
led to the development of plasma centrifuges for element and isotope 
separation [1.17-1.21]. In such a device, ionized particles are sent 
into rotation by the Lorentz force. Most devices of this type use two 
ring electrodes of different radii, thereby creating a radial electric 
field. Rotation can be observed when this radial current density is 
combined with an axial magnetic field. Plasma centrifuges can attain 
rotational velocities several orders of magnitude higher than their 
mechanical counterparts. Consequently, light and heavy ionic and 
atomic isotopes can be eens ated more efficiently. 

These proposed plasma centrifuges can utilize one of two types 
of working fluids, a low density collisionless plasma or a high density 


collision dominated plasma. The low pressure, fully ionized plasma 
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requires a large electrical power input to produce the high degree of 
ionization necessary for practical isotope separation. However, the 
highly ionized gas centrifuge has reduced velocity losses at its walls 
due to the absence of ordinary hydrodynamic boundary layers. This 
accounts for the extremely high rotational velocities observed. 

In the high density gas centrifuge, only a small fraction of the 
isotopic particles need be ionized. The Lorentz force drives the ions 
which in turn drives the neutral particles through inter-component 
viscous forces. The large number of collisions and the viscous boundary 
layers at the containment walls interact to reduce the rotational 
velocity somewhat. Due to the lower input energy, the gas temperature 
will be substantially lower and this will have a favourable effect on 
isotope separation. This will also lead to a more energy efficient 
device. 

At present these centrifuge devices are being built and theor- 
etically analyzed. The analysis of these gas discharges are based upon 
magnetogasdynamic equations. Thus far, only simplified two-dimensional 
solutions have been attempted. With these first attempts some progress 


has been made to confirm the usefulness of these types of devices. 
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1.5 MGD STABILIZATION CONCEPT 


As reported earlier, gas discharges operating in high power 
lasers suffer from numerous types of bulk and micro-instabilities. It 
was proposed that a laser gas discharge could be subjected to a Lorentz 
force to create a high velocity rotating plasma. This technique would 
then be employed to rapidly convect thermal instabilities out of the 
discharge volume by reducing the characteristic gas cooling time ie 

The viability of this concept would only prove successful pro- 
vided that a proper magnetic field profile be attained. In addition to 
a rotational driving force, a theoretical model for a confining force 
was developed. With the proper configuration of magnetic field, an 
elaborate magnetic potential well could be constructed to confine a 
rapidly rotating discharge plasma. In this manner, the tendency for a 
glow-to-arc transition at the electrode extremities could be suppressed. 
Thus, the need for specially shaped or profiled electrodes would be 
eliminated. 

It is expected that viscous forces will effectively couple the 
fast moving ionized particles to the slower moving neutral particles. 
However, at the outer edges of the potential well, the driving force 
terms would be weaker. This would then allow a decoupling of these 
particles and result in more free movement for the neutral particles 
in this region. The neutrals, now only under the influence of a cen- 
trifugal force, should tend to fly off tangentially. This concept would 
create an ideal situation in which ionized particles are captured by 
the magnetic field, while neutral particles would be free to travel 


through the containment region. 


ms 


$1 este? Fant bing ies 40! en 9 
wins ls & 92 butondue set bits sevedosth Hip asst «i 
bt pow suptndsey anny. ‘ire! iutasion Si potoy doh ‘ys 

aly Fo die gets i idetent an towns Esha bey 


sia ‘ 
ar 


<< 


9 Omn2, BAT TooD Zap nite Segoe ery: ont #8 oe 
cn lwWiersoove sverg yin bilew HsSWO9 as ene ee 
ot mOhtthbs 1 ‘fone $16 °390 af} Fearn bis! 31 ~_ aint i 
aswel griniinos 4 yo? j,Sp0m iso snbeth By i aobyt fs | 
16’, Ofat? srgeapem to noriswpr te ~gtloty, Eo aay | 

5 snrined of Sereo sd biuus Tia rer Instod soe rod 

s “ot vansbrist anit; annem eras ae master ‘spugiozth aca 
2ermawe ad’ bt yaa. eshte yy, stiotyost EL Ng hs. ott aie 9 be 


ad bluow esbenfoats hel hong 10 ‘haber LAR Pstoste.’ 70 boon ot 8 


7 
G 
y | i! Fi _ r a 


" o% 


| ads slayun wtayeandate “fit tw Le et audovk sais pesaaiis oF Pa 
estates Tertuan pel Mow newigyhe, eid ad, asipiasy besians wise Sa 
gone? yatwith ey . (Tow reidigton: gilt 10 esebe vatuo art +5 4 evan 
searld to guTiquodsb 6 Wollis aed. biwow aid? <Vsolaiaw sd) bTuvom —p 
esi ai rreg Peatuen oft vot SismBVeR 9547 avon np Sfuzsy Beas, wre 
“nao 6 to aoneGTani ond “be ine bh stiwen ont aad ihe abit “7% “ 
bivew Jqanno2 ataT ul a an Te Vit oF bred. bhvbre So agus” " 
Yd beyuioss evs ealotitsg ogni dotiw wt Hor deud 2 teabh ne oseenD es 
loysit of a37) od. bfuaw zealot tieg . Lengwen, 8 {'TiW, shia? or jengem ont inl 


Pa tuganibetnos arid stougaas iy 


The goal of this project is to produce a highly stabilized laser 
discharge, utilizing the concepts of a rotating plasma within a pot- 
ential well. This new dynamic stabilization system should be much 
Simpler than any other existing approach. If the rotating gas fica. 
velocity is sufficiently high, there should be little or no need for 
an external blower system. With the appropriate heat exchanger design, 
the gas may be allowed to circulate on its own accord. 

This unique approach to gas discharge stabilization seems to 
Satisfy the necessary requirements for producing a uniform glow dis- 


charge that should be highly immune to various types of instabilities. 


1.6 THESIS PROLOGUE 


The primary objective of this thesis is to define the fundamental 


principles for a new generation of high power gas discharge lasers 


based upon a MGD stabilization concept. It has long been known that 


magnetic fields can significantly alter the characteristics of a plasma. 


Furthermore, it is evident that such fields could be utilized to 
provide a stabilizing mechanism in which a uniform glow discharge could 
be maintained. It was then essential that the different types of gas 
discharge instability processes be analyzed to determine their signifi- 


cance and role. As will be explained in Chapter 2, a gas discharge 


contains a wide variety of instability modes that can develop with time. 


Each mode has its own particular driving mechanism and characteristic 
critical wavelength. Many of these modes are coupled to one another, 


thereby complicating a detailed analysis. 
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Some of the major components used to test the theory of MGD 
stabilization will be discussed in Chapter 3. Due to the harsh oper- 
ating environment found within a high power laser, the physical com- 
ponents are forced to undergo severe thermal and structural tests. In 
addition, the water-cooled electrode and siecemnnanes assemblies are 
required to operate within a vacuum chamber. Using proper materials, 
these problems were overcome. However, to achieve a uniformly distrib- 
uted glow discharge, the electrode structure underwent several design 
changes before a satisfactory system was developed. 

In order to understand the complex processes occurring within 
the gas discharge, a computer code has been developed. The computer 
program MAGIC (MAgnetoGasdynamic Implicit Code) is designed to evaluate 
time varying changes in pressure and velocity. In MAGIC, a system of 
MGD equations was solved by the alternating direction implicit (ADI) 
technique on a variable Eulerian mesh. Chapter 4 developes the MGD 
equations necessary to form an accurate model of a partially ionized 
glow discharge. Furthermore, the gas transport coefficients, found in 
these equations are described. 

The MGD equations are finite differenced in orthogonal curvi- 
linear coordinates to allow program symmetry and to accommodate a 
variety of coordinate systems. The ADI method described in Chapter 5 
allowed the use of longer timesteps, reduced computing costs, and 
improved program stability. After time and spatial differencing, the 
nonlinear set of MGD equations are simultaneously solved using a Newton- 
Raphson iterative procedure. The set of equations are iterated until a 


convergence criterion is satisfied. 
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The results of the computer similation on a transverse electrode 
geometry are presented in Chapter 6. These computed values are then 
compared with previously determined models and experimental data. In 
addition, the performance of the MGD stabilization technique in producing 
a uniformly distributed discharge is evaluated by using high speed 
motion piCbures. To demonstrate the improvement in discharge power 
density, the MGD stabilization technique is compared with conventional 
gas discharge systems. The results are graphically illustrated for a 
variety of parameters. 

A detailed explanation of the physical processes occurring within 
a gas discharge is discussed in Chapter 7. Also analyzed are the changes 
in the structure of the glow discharge under the influence of a mag- 
netic field. Furthermore, an explanation of the role that the non- 
uniform magnetic field performs in the MGD stabilization concept is 
carried out. 

The final chapter also summarizes the many other possible aspects 
that need to be explored. Moreover, it describes some problems still 
remaining to be resolved before the potential of this new discharge 


device is maximized. 
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CHAPTER 2 
DISCHARGE INSTABILITY ANALYSIS 
2.1 INTRODUCTION 


The analysis and understanding of a low density, non-magnetized 
glow discharge is an extremely complex problem which involves a series 
of energy transfer and particle conservation interactions. Examination 
of these interactions is of special interest to those in the laser dis- 
cipline who Scrutinize instability processes. Their efforts have been 
focused upon discerning the roles played by different gas mixtures, 
external ionization sources, temperature and density fluctuations, and 
gas transport processes. This chapter will attempt to examine some of 
the roles and determine how they affect glow discharge behavior. 

Plasma instabilities are the result of small amplitude pertur- 
bations that occur within the gas discharge. These fluctuations in 
temperature and/or electron density excite a number of different wave 
modes. These modes and their stability have been thoroughly analyzed 
by Haas [2.1]. Associated with the instability modes are characteristic 
times which collectively span the time scale 10°/° to 10°*s. Instab- 
ilities produced in a large volume, high energy discharge plasma are 
Strongly connected to these modes and are easily excited by a variety 
of mechanisms. A comprehensive list of all processes and characteristic 
times ts can be found in Table 2.1. The symbol ie represents the 
particle density of the jth species. Also represented are the rate 


coefficients of the jth collisional process Ks . This term is defined 
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as the normalized electron-molecule collision frequency, that is 


The main processes that can occur within the discharge can be 
grouped into four categories: space-charge relaxation, plasma kinetic 
processes, transport processes, and collisional energy transfer. An 
analysis of the space-charge relaxation process by Haas [2.1] has 
determined that this mode poses no overall stability problems. The 
-remaining groups will be briefly examined in order to investigate 
their relative strengths and cumulative contributions to the onset of 


instabilities. 


2.2 PLASMA KINETIC PROCESSES 


The characteristic times associated with the plasma kinetic 
‘processes measure the temporal evolution of particle production and 
loss processes, and are typically in the 10 °-10 ‘s range. These 
modes along with their individual rate coefficients are: electron 
impact ionization of molecules (k.) and electrically excited species 
(ke), electron-molecule attachment (kK); detachment by neutral impact 


(k.,), electronic excitation (k"), quenching of electric species (ke) 


7 
two-body electron-ion recombination (Ki) and positive-ion negative-ion 
recombination (kr). 

Working from the Boltzmann equation, an energy balance can be 


obtained for electrons within a spatially uniform gas in the presence 


of a steady electric field E. This electron energy equation compares 
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the rate of energy deposition from the electric field (electron Joule 
heating) with the rate at which electrons exchange energy through 
collisions with molecules. The resulting normalized electron-molecule 


collision frequency v_g/"; can then be defined as [2.2] 


Whe: 5 po 
CE (22 
Ge ra [ cf E)@, «(546 - (2.1) 


Here, the various rate coefficients are determined by averaging the 
proper cross section Q,;(5) over the electron energy distribution 
function f(&). The symbol € represents the electron energy in eV, 


through the definition 


Other charged particle collision frequencies or rate coefficients are 
determined by a similar analysis using the appropriate distribution 
function. 

Since the electric field provides the force which accelerates 
the electrons, the ratio Z/n becomes proportional to the electron 
energy. Thus, the rate coefficients described in Eqn. 2.1 are strongly 
dependent upon the parameter £/n. 

In a subsequent chapter, the effect of a magnetic field on both 
the electron distribution function and the rate coefficients will be 
investigated. 

Due to electron impact dissociation of C0., and No » a variety 
of minority molecular species are created. The formation of such 
products as CO, 0, NO and a host of positive and negative ions have 


little effect on electron-molecule vibrational coupling or on 
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vibrational-translational relaxation. However, they can influence the 
production and loss of electrons through recombination, attachment, 
and detachment. These dissociative products can accumulate to signif- 
icant levels, even exceeding the electron density within typical fluid 
residence times. Some of the more studied reactions occurring in C0, 
laser discharges are listed in Table 2.2. 

Of the primary species found in most C0., laser mixtures, C0. 
has the largest rate coefficients for electron impact ionization and 
dissociative attachment. Since both these coefficients are strongly 
dependent on electron temperature, the electron-C0,, kinetics dominate 
ionization and attachment processes in laser discharges. For 
“comparison, several important C0, rate coefficients are plotted 
against electron temperature in Fig. 2.1. Through experimental 
analysis, Nighan [2.4] has determined that most discharges are 
dominated by detachment effects. Accordingly, the indicated discharge 
operating conditions were found to be quite insensitive to the electron 


energy distribution, fractional ionization, or positive-ion species. 


2.2.1 Effects of Detachment Kinetics 


A survey of the literature has revealed that detachment kinetic 
processes are largely responsible for determining the behavior of a 
self-sustained gas discharge. In particular, if the density of CO jis 


large, then the reaction 
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TABLE 2.2 Important Chemical Reactions in CO,-N,-He Discharges and 
their Rate Coefficients [2.3, 2.4]. 
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Fig. 2.1 Effect of Electron Temperature on the Rate Coefficients in 


a C0. Laser Gas Mixture. 
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will keep the negative-ion density small and positive-ion negative-ion 
recombination will be overshadowed by detachment kinetics. As a result, 
the characteristic time for the negative ions to respond to a disturbance 
is much less than the time for the electron density to respond. Hence, 
during an effective detachment process electrons are released from 
negative ions so that n/n, << 1. Consequently, for this case, negative- 
ion kinetics will have little or no influence on either steady-state 

or transient electron production and loss processes. An examination 

of the electron continuity equation then reveals that ionization 


balances electron-positive-ion recombination, that is, 


e 
kK = gS: , 


In this situation, the discharge can noramlly be sustained at a lower 
value of electron temperature as indicated by point B in Fig. 2.1. 

Examination of the situation when detachment processes are 
negligible, reveals that under self-sustained operating conditions, 
the electron production due to ionization must balance electron loss 
due to attachment, that is, k. es: K . This occurs for an electron 
temperature of 1.7 eV (point A) for the conditions of Fig. 2.1. In 
this situation, the negative-ion density can become comparable to the 
electron density [2.4]. Such a coupling of electron and negative-ion 
kinetics will likely lead to an attachment induced ionization 
instability. 

The normal region of operation for a self-sustained gas dis- 
charge, lies between points A and B of Fig. 2.1. Since CO is an 


effective detaching species, easily produced by dissociative 
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recombination or attachment [2.3], the discharge operating point is 
Significantly influenced by its fractional concentration. For the 
case of sealed or slow-flow discharges containing C0, >» in which CO 
has had sufficient time to accumulate, detachment processes are 
common. Thus, negative-ion processes are likely to be insignificant. 
However, in fast-flow convection dominated discharges, the brief gas 
residence time (< 10°*s) does not allow the CO fraction to become 
Substantial. Consequently, the negative-ion density starts to rise 
and thereby begins to influence plasma stability within these non- 


recirculating systems. 


2.2.2 Ionization Instabilities 


Ionization instabilities involve the production of two possible 
modes. The first is identifiable with the electron production and loss 
kinetics (ionization mode), and the second is associated with the pro- 
duction and loss of negative ions (negative-ion mode). Disturbances 
in the electron and negative-ion densities are normally controlled by 
electron-ion recombination and by negative-ion detachment processes. 
Both these effects exert a damping influence; recombination reduces 
the magnitude of electron density fluctuations, and detachment reduces 
the negative-ion concentration. However, these damping processes are 
not always successful in their attempts to eliminate ionization instab- 
ilities. As illustrated in Fig. 2.2, several different scenarios are 


possible depending upon the dominant rate process. 
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As shown in Fig. 2.2, an initial local positive disturbance in 
electron density is accompanied by a drop in electron temperature. 
This effect is the result of the electron temperature dependence on the 
energy exchange collision frequency vy? being stronger than its 
dependence on the momentum transfer collision frequency. Hence, 
electron energy transported by momentum out of a disturbance region is 
small compared to the electron energy loss due to collisions [2.4]. 

When negative ions are absent, ionization and recombination 
kinetics dominate electron production and loss (case A). For this 
particular case, Fig. 2.1 reveals that within a self-sustained plasma, 
k << ke, . Therefore, the electron density will decrease since 
recombination effects have a much stronger influence on electron 
density than due to ionization. Consequently, the electron-ion 
recombination process describes an effective local electron loss 
which exerts a stabilizing influence on the gas discharge. Sub- 
sequently, it can be concluded that the ionization mode, in the 
absence of negative ions, is stable. | 

When negative ions are present (case B), an ionization instab- 
ility will occur provided that K > k and that the electron 
attachment rate increase with electron temperature [2.1]. Since the 
ionization and attachment times are comparable in magnitude, the 
negative-ion density will increase. Naturally, recombination processes 
will become less significant in favour of detachment. The instability 
occurs because the electron density and temperature fluctuations are 
out of phase. In this circumstance, the electron production rate due 


to ionization and detachment, exceeds the loss of electrons by 


q sepegoo nobtent tape ATR 


ae ee ee 
srwteseaniod rordnete nt ore one bee | 
ait Wo soreabregeb: lad noses nts tes it 
2d? nerd NbunYte ‘pnbed . " . pa" 
sonal. tonebpen? norahikes ‘aon was TSIM: OR 

ef nora : ane WVYIZID § ™ Bi 7) va Shithon ip orpegchedt a 
PRS] ano t< tifes o2 sub aaat ovens ae 
no hhaotinpsey bie. hotresinat nerds s46 me 
orde yea (4 9250) 2207 OME ni tonbioNa ‘noah 
tiiemtg benisdeua~thae 6 niatin sedd Stesyey fe % “ah seen ng 
asavia~Seaetoeb [Fw \s renan La pea I] ae cratered : oo 
NOVYioSsia oO ad nevitn? vspnetye AOR: ‘Bb eve sestites 
not-~NoV joa ke: ont . idraypaario? “Wot ein oat say 


220) nontodie lnsar svidoat ta) meee 198 


“due pagventooth 2ap St) ne spat tith ‘pris tnahde +z 

ed at! abom «itserday eit tart bit 05; ode oe. 

Paty -Stste eh .ghot silt zsigem rs 

-dsdent worsesingt ne. .[d naa} sa ahaa aiyt enon ei | 

ouisosis sit shy ons ay i daria bet were woo tw a 
eit sant? .[f28) SUSE VRS: aordsate ‘Woe aagsnant over: Senetsn 

SHI . Shur tipom ni wee 246 at SANIT DEI bya: -notaestnol | 

Bt siont S41 Gt a wane 

walt Fdsaent ot gest i ‘oye wi suse Mhagte ‘ae 5 Atte 9 1" 

om amottaatoutt evs eroqad ‘bi iiteiet Hortons silt boul ews tas im ‘ 

eb ater noha oubarn nontoete ita sagonissari 1 ates’ ot wait Fo tuo: Bs 


an 
ed ‘gnondaefe to een ane sana ,Snantanteb a. 0} Ae boi ca 


rye (i 


28 


recomination and attachment during an electron temperature decrease. 

As a consequence, this mechanism will produce a local electron density 
increase which will create a positive feedback instability process as 
illustrated in Fig. 2.2. Nighan [2.4] has shown that the growth rate 

of this instability is strongly peaked in the direction of the applied 
electric field. This explains the striated and/or constricted manifest- 


ations so often seen within gas discharges. 


2e2 danke fects of External Preionization 


Kinetic processes which involve electron temperature dependent 
ionization and attachment rate coefficients, produce micro-instabilities 
which are independent of the discharge power density and gas pressure. 
One common technique to overcome this problem is to utilize an external 
ionization approach in the form of an electron-beam or photoionization 
source. With this method, charge production in the plasma is decoupled 
from the applied electric field. This decoupling can significantly 
reduce the effect of electron temperature variance on electron and 
negative-ion density fluctuations. In other words, the independently 
controlled ionization source allows the electron temperature to be 
adjusted to a value independent of electron production and losses. 
However, for this method to be successful, the contribution to the 
ionization rate due to direct low energy electron impact by drifting 
electrons, must be negligible compared with that provided by an 
auxiliary ionization source, that is, S >> nk, . Thus for optimal 
performance, the ratio #/n is lowered such that the discharge operates 


in the externally sustained region as shown in Fig. 2.1. 
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The electron production mechanism is highly dependent upon the 
dominant detachment process. Since the presence of detachment processes 
has a stabilizing influence, less external ionization is required. 
Consequently, the discharge can be sustained at a higher value of 
electron temperature. However, if detachment processes are less sig- 
nificant, the electron temperature must be maintained to the left of 
point C in Fig. 2.1 in order to avoid excessive negative-ion formation. 
In this manner, the plasma will be dominated by the stabilizing influence 
of recombination, provided the condition Kk. << (n,/n)k,, can be satisfied. 

The growth times for the ionization and negative-ion modes are 
Short compared to the characteristic times for other modes. Con- 
sequently, when they occur, the ionization and negative-ion modes 
dominate plasma behavior. The subsequent application of an external 
ionization source helps to control the plasma ionization mechanism and 


thereby reduces their destabilizing effects. 


2.3 TRANSPORT PROCESSES 


Neutral gas processes have somewhat longer characteristic times 
than the time associated with charged particle dynamics. Changes in gas 
temperature and density can cause fluctuations in neutral gas properties 
which can couple to the charged particle kinetics and electron energy 
transfer eres The coupling that does occur primarily depends upon 
which rate coefficients are dominant [2.1]. 

The presence of negative ions can trigger adverse effects with 
regards to discharge stability processes. Since detachment is strongly 


temperature dependent, a local gas temperature increase leads to an 
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increased electron production rate. Under these conditions, a large 
local increase in electron density causes additional heating of the 

gas by way of elastic collisions. This further stimulates the detach- 
ment process and thereby provides the positive feedback mechanism 
necessary for this process to be unstable. Thus, an initial disturb- 
ance in gas temperature will lead to an instability growth provided 
thermal conduction does not dissipate the resulting local concentration 


of heat. 


2.4 COLLISIONAL ENERGY TRANSFER PROCESSES 


Some of the other important collisional processes that occur 
are the thermal and the vibration relaxation modes. These modes are 
created by perturbations in local gas temperature and density which 
couple to the electronic and vibrational properties of the gas. An 
unstable thermal mode is developed when vibrational eae aeton causes 
local gas heating. The rise in temperature forces the V-7 relaxation 
time Top to decrease. This process prevents a rapid energy release 
within the gas and the stabilizing influence of thermal conduction will 
be unable to transport this energy away from its point of origin. The 
net effect is a runaway condition in which the gas temperature continues 
to rise. The gas density will then begin to decrease and the local 
electrical conductivity quickly rises to create a high concentration 
of current flow. The time constant for the evolution of this arclike 


current filament is t typically 10 “1004s, 
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This instability mode is likely to occur for all typical high 
power density discharges with some degree of ionization and a large 
electron temperature to gas temperature ratio. This mode also has a 


characteristic wavelength such that if » > Xr , the thermal mode 


critical 
is unstable. As the perturbation in temperature grows, the ability of 
thermal conduction to dissipate local hot spots decreases. Consequently, 
the critical wavelength decreases and current filaments become narrower, 
thus promoting the formation of an arc. 

The thermal modes are insensitive to external ionization sources 
Since they are relatively unaffected by charged particle kinetics. 
However, large transport processes can be used successfully if the 
critical wavelength is not too large. This necessitates the use of 
high flow velocities to provide heat removal and to convect locally 
unstable fluid elements from the discharge region. These effects must 
occur in a time less than that required for a glow-to-arc transition. 
The results of Nighan [2.6], have indicated that a convective laser 
having a power density of 10 W/cm? will require a gas residence time 


below 1 ms. Only in this manner will convective type CW gas lasers 


operate properly. 
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CHAPTER 3 


SYSTEM DESIGN 
3.1 INTRODUCTION 


The main features of the experimental apparatus used in the MGD 
stabilization concept are illustrated in Fig. 3.1. It consists of five 
major components: a vacuum system, a thyratron pulser, an appropriately 
designed electrode, several types of electromagnets and a 3-phase 
power supply. The overall system was designed with sufficient flex- 
ibility to accommodate a wide variety of gas discharge experiments. 
These included the testing of various electrode and magnet designs. 

The apparatus was not designed to operate as a laser, but rather as 
an experimental device used to test new ideas and components neces- 


sary for a more advanced laser. 


3.2 VACUUM SYSTEM 


Using available material, the gas discharge chamber was built 
from a PVC pipe 40 cm in diameter, 130 cm long, with a wall thickness 
of 2 cm. In order to allow easy assembly, operation, and modification 
of the vacuum system, demountable seals were employed. The open ends 
of the chamber were covered with clear plexi-glass plates to allow good 
visual observation of the gas discharge. In order to seal the system, 
Viton A elastomer O-rings were installed where necessary. Additional 
plexi-glass access and diagnostic parts were added at a later date to 


allow experimental gain and temperature measurements. 
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To evacuate the vacuum chamber, a Welch rotary mechanical vacuum 
pump, with a free air displacement of 500 1/min, was attached to the 
system. After several hours of pumping, the lowest attainable pressure, 
with this roughing pump, has only been 2 x 10 °torr. At this point, 
the chamber was then back-filled with the appropriate gas mixture. A 
mechanical diaphragm gauge was then used to measure the partial pressure 


of each gas species introduced into the vessel. 


3.3 ELECTRODE MATERIALS AND DESIGN 


To initiate an electric discharge, two metallic electrodes, 
each at a different potential, are required. Experience has shown 
that electrode materials and their configuration play an important 
role in determining the stability and general overall performance of 
an electric discharge. For a number of years, TEA lasers have relied 
upon the uniform electric field aspects of their Rogowski profiled 
electrodes. However, unlike TEA laser systems, the uniform field 
characteristics of a CW laser does not play a dominant role in discharge 
Stability. Of more importance are the dynamic parameters such as; 
flow uniformity, boundary layer effects and gas residence time be- 
tween the electrodes. Consequently, a variety of electrode shapes and 
materials have been evaluated for each CW laser system. It is common 
practice, however, that in most realistic devices care be taken to 
eliminate sharp edges and other prominent surface irregularities. 

Such imperfections can create high electric field concentrations and 
promote discharge instabilities. 

In an attempt to arrive at an efficient device design, several 


different electrode materials were examined; most notably aluminum, 
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graphite, steel and copper. In order to facilitate the choice of an 
optimum material, a survey of each was performed, with respect to its 
durability, work function (both photoelectric and thermionic) and 
behavior in a gas discharge. 

The first electrode design was constructed from aluminum. After 
a relatively short period of testing, it was observed that this alum- 
inum electrode performed quite poorly, particularly with respect to its 
durability when subjected to a glow-to-arc transition. More precisely, 
it was observed that a high current arc could easily damage the electrode 
surface by arc pitting. This pitted surface then promoted the forma- 
tion of hot spots. With further influx of energy these hot spots became 
potential sites for streamers and subsequently more arcs. The occur- 
rence of this phenomena was most likely due to the high electric fields 
and correspondingly excessive emission of electrons emanating from 
sharp edges along the pitted surface. A detailed inspection of the 
aluminum electrode surface revealed a myriad of tiny scratches, all 
running in the direction of the rotating plasma. It is most probable 
that these scratches are the result of a large number of streamers 
etching a pattern, as they rotated with the plasma. It was clearly 
apparent from the results of these experiments that an aluminum elec- 
trode could not withstand prolonged use in this type of electrical 
discharge. 

A carbon or graphite electrode was then tested because of its 
proported durability and slightly resistive qualities. It was antic- 
ipated that a resistive electrode would assist in spreading the dis- 


charge over the electrode surface, thereby resulting in a more uniformly 
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distributed glow. A slight spreading of the discharge was indeed 
observed but it was only confined to the electrode surface and conse- 


quently did little to improve the bulk plasma region. 


Several other qualities make graphite favourable as an electrode. 


For instance, graphite has the ability to withstand very high tempera- 
tures; it is available in a variety of composite forms; and is easily 
machinable. With a suitably shaped graphite electrode, several experi- 
ments were performed to determine its durability in an electric dis- 
charge. Unfortunately, it was also discovered that an arc or high 
intensity streamer could easily damage and pit the soft graphite 
surface. After continued operation, the chamber apparatus became 
covered with a fine layer of carbon dust, due to sputtered graphite 
particles. It is believed that the porosity and softness of the 
graphite material is at fault. Choosing a graphite-metallic com- 
posite would probably alleviate this problem. 

A steel electrode was then tested and observed to perform 
moderately well in the high current gas discharge. An examination of 
this metallic electrode surface revealed scratches similar to those 
found on the aluminum electrode, though not as deep. However, the 
use of a steel electrode was not advantageous since its presence could 
disrupt the magnetic field pattern necessary for discharge stabiliza- 
tion. In order to overcome this difficulty, a non-ferromagnetic 
Stainless steel could be substituted. 

Ultimately, copper was finally chosen as the material best 
suited for electrode fabrication. The copper electrodes were non- 


ferromagnetic and were less prone to arc damage than aluminum, steel, 
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or carbon. In addition, the thermal conductivity of copper is higher 
than the other materials, making it ideal for transporting heat from 
the electrode surface. Table 3.1 lists the various work functions of 
the electrode materials tested. Upon examination, it is seen that 
copper has the lowest thermionic work function in addition to a mod- 
erate photoelectric work function. These results indicate that copper 
is a good source and emitter of electrons which are essential to dis- 
charge formation. An additional benefit derived from using copper is 
the variety of shapes and forms that are available, thus permitting the 
design of several types of electrode structures. In order to confirm 
the usefulness of MGD as a viable stabilization technique, three types 
of electrode structures were constructed: transverse, co-axial, and 


fluid ballasted. 


Photoelectric thermionic. 


Element work function (ev) 
Al 2.98-4.43 - 
C 4.8] 4.39 
Fe 3.91-4.7 4.04-4.77 
Cu 4.07-4.8 3.85-4.38 


Table 3.1 Work Functions of Tested Electrode Materials [3.2]. 
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3.3.1 Transverse Electrode Structure 


Initial experiments in the transverse geometry employed two 
electrodes constructed from short sections of copper tubing. The tubes 
were arranged in parallel and connected to water inlet and outlet 
headers at either end. The 6 mm diameter copper tubes, selected for 
the electrodes, were equally spaced 6 mm apart to provide 50% transmit- 
tivity to a Flowing gas. During actual operation, cold tap water 
flowed through the electrode structure so as to remove excess heat 
generated by the gas discharge. 

Observation of the rotating plasma indicated excessive tur- 
bulence at the cathode. It appeared that the copper tubes were creating 
a large nonuniform boundary layer which severely inhibited plasma flow. 
This effect was highly noticeable from the ragged appearance and jit- 
tering motion of the cathode glow layer. Furthermore, it was probable 
that a number of stagnant gas cells were also present. These cells 
eventually coalesced into numerous hot spots as more energy was fed to 
them. As time progressed, these spots contributed to the formation of 
streamer instabilities. 

The tubular cathode was later replaced with a water cooled copper 
plate. The smooth plate eliminated turbulence and created an electrode 
Surface over which the cathode glow layer could move unimpeded under the 
direction of the Lorentz force. This new cathode was constructed from 
a 6 mm brass plate with a machined groove, 3 mm deep, meandering across 
the surface. A 3 mm copper plate was then soldered onto the machined 


brass surface, forming a leak tight water channel. Since the initial 
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installation, the cathode plate Structure has provided reliable, trouble- 
free performance with a minimum of maintenance and care. 

Inspection of the gas discharge near the tubular anode indicated 
Only minor turbulence. It would seem that a major difference in 
behavior exists between the tubular cathode and the tubular anode. There 
are two probable causes, responsible for such a wide difference. First, 
due to a reduced Lorentz force at the anode, a much lower gas velocity 
is produced and, correspondingly, a lower level of turbulence exists. 
Of more importance though, is the absence of a highly ionized region at 
the anode. The space-charge layers at the cathode are highly suscept- 
ible to a disturbance, more so than one occurring within the diffuse 
glow near the anode. The primary reason for this is the highly organ- 
ized nature of the cathode fall region and its importance in maintaining 
the overall structure of the glow discharge. 

To further assist the maintenance of a stable glow discharge, 
the copper tubes of the anode were bent into a Rogowski prone: It was 
later discovered that this profile also assisted the recirculation of 
gas through the anode. To minimize any disturbance to the transverse 
flow, the electrodes were separated by thin, adjustable ceramic posts. 
These and other details of the overall assembly of the transverse 


electrode structure can be seen in Figs. 3.2 and 3.3. 
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Fig. 3.2 Magnetically Stabilized Transverse Discharge Geometry. 
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3.3.2 Coaxial Electrode Structure 


Following the successful example of the transverse electrode 
geometry, a coaxial structure was then investigated. As in the previous 
geometry, the electrodes were made of copper. The cathode consisted of 
a 9 cm diameter pipe with a spiral groove machined into the outer sur- 
face. Another slightly larger diameter thin-walled pipe was placed over 
the inner pipe, thereby leaving a spiral water channel in the middle. 
The ends were soldered together and provided with Swagelock fittings 
to supply the cooling water. 

The anode was comprised simply of a 1.25 cm diameter copper pipe 
equipped with Swagelock fittings on either end. To vary the electrode 
gap, larger diameter pipes were used as the anode. Each anode was 
Supported in the center of the cathode by an insulating frame. The 
complete coaxial electrode structure was then mounted inside a solenoid 
as illustrated in Fig. 3.4. 

This electrode structure performed extremely well and had no 
major difficulties. Turbulence within the rapidly rotating discharge 
was not observed, due to the presence of smooth electrode walls. In 
addition, electrode edge instabilities did not present a problem, as 
the magnetic field kept the plasma confined to the central portion of 


the electrodes. 
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3.3.3 Multi-pin, Ballasted Anode 


The previous two electrode designs were only partially success- 
ful in their attempt to create an extended area uniform plasma distri- 
bution. It was determined, through the use of high speed photography, 
that only a rapidly rotating plasma column was created. The film 
revealed that the rotating plasma column was simply running across the 
continuously conducting anode surface. | 

As a result of the difficulty encountered with the continuous 
Surface anode, a new electrode structure, based on the multi-pin 
concept [3.3], was constructed. Now, instead of a continuous plate 
or tubular surface, the new anode was constructed such that the overall 
Surface was subdivided into an array of independent sub-electrodes, each 
with its own current limiting feature. This was accomplished by employ- 
ing an internal electrolytic solution to balance the current fed to 
each sub-electrode pin. In this manner, a small electric discharge 
was created at each pin, thereby spreading the applied current over 
the entire pin array. Early SEGARA by others, to create a uniform 
current distribution used discrete resistors as the ballasting com- 
ponents. However, these components were often prone to overheating 
and failure. 

The initial multi-pin electode concept is displayed in Fig. 3.5. 
This structure employs the use of an array of bent stainless steel pins 
mounted such that one end was submerged into a channel containing the 
recirculating electrolytic fluid. The sub-electrode pins were bent at 
right angles to allow cool gas to flow through the anode and replace 


gas lost by the centrifugal pumping action of the rotating gas 
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discharge. Furthermore, the long bent pins prevented damage to the 
electrode structure by the hot ionized gas. 

As discussed in [3.3], the electrolytic fluid best suited for 
prolonged use within a current ballasted electrode consists of KACO. 
dissolved in distilled water. As a solvent, distilled water is free 
of chemical and mineral contaminants often found in ordinary tap water. 
Potassium carbonate was advocated as a solute because it does not 
produce an excessive amount of gas, nor is it strongly corrosive towards 
metals. In addition, the solution resistivity is easily controlled 
and is nearly linear when conducting a large current. With the proper 
solute concentration, the electrolytic solution will create a uniform 
current distribution without the need for external ballast. Thus, the 
total power lost in the system can be minimized. Optimal performance 
was achieved with an electrolyte concentration of 0.5 g/1. At this 
concentration, the VI characteristics of the electrode is shown by 
Eigen oO. 

Due to the large amount of current that can be passed through 
the system, a significant temperature gradient can exist within the 
electrode structure. It has been determined that the resistivity of 
the electrolytic solution decreases almost linearly with increasing 
fluid temperature. This can be explained by the fact that the dis- 
solved electrolyte is more readily dissociated at a higher temperature, 
and thus creates a higher charge carrier density. To maintain a 
constant temperature throughout the electrode, a high speed pump is 
used to recirculate the fluid through a heat exchanger. This process 


removes excess heat and in addition forces gas bubbles, generated by 
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electrolysis, out of the system. As a result of this action, the 
electrolytic solution provides both cooling and current distribution 
to each element constituting the anode. 

As discussed earlier, the sub-electrode pin material was of 
primary importance in determining the performance of the gas discharge. 
Past experience with gas discharge electrodes has shown that copper 
and stainless steel worked best. However, with the multi-pin concept 
there are additional constraints to be considered. Since the pin 
material had to be in contact with the electrolytic solution, some cor- 
rosion was to be expected. It was discovered that corrosion by the 
electrolytic solution formed a prominent oxide layer on the copper 
pins. As the oxide layer formed, it began to inhibit the current flow 
to the pins. It was determined by Nam [3.3], that stainless steel pins 
had little or no oxide film on their surfaces exposed to the electrolyte 
solution. In addition, stainless steel had a very Tow corrosion rate. 
This could be determined By passing a specified current through the pin 
and measuring the weight lost over a period of time. These results sug- 
gested that stainless steel was ideal for a multi-pin electrode, which guar- 
anteed good performance and a longer lifetime over other metallic pins. 

The current carried, through the electrolytic solution, between 


the busbar and the sub-electrode pin can be expressed by Ohm's law 


Aw (3.3) 


where i, is the terminal voltage between the busbar and the pin, E 4 is 
the electrochemical decomposition potential of the electrolytic solution 


and o., is the fluid conductivity. From this equation, it can be 
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concluded that the current is also controlled by the water contacting 
Surface area of the pins, A as well as the gap spacing, d, between 
the busbar and the pins. 

With this knowledge, the multi-pin electrode can be resistively 
profiled by selectively adjusting the differential gap spacing, d, over 
parts of the pin array. This simple technique was used to increase 
the effective sub-electrode resistance near the electrode extremities. 
Consequently, this process was very effective in removing the tendency 
for a glow-to-arc transition to occur along the electrode's edges. 

Perhaps the best feature of the multi-pin ballasted anode design 
is its ability to suppress the formation of arcs and streamers. Should 
one electrode pin begin to conduct more current than its surrounding 
neighbors, additional gas bubbles, generated by this accelerated elec- 
trolysis action, momentarily surrounds the pin. These gas bubbles 
reduce the effective water contacting area, causing a decrease in fluid 
conductivity. This naturally reduces the current Anse carried to 
the pin, the arc is suppressed and equilibrium is restored. 

Another useful characteristic of this electrode design 4 its 
tendency for self-optimization. After final assembly, each pin 
experiences slightly different operating conditcine due to non- 
uniformities of construction and slight differences in pin length 
which can alter the values of d, and A, As a consequence, some pins 
initially conduct more current than others. This causes increased 
electrolysis and selectively accelerates pin corrosion. After a period 
of time, the individual electrode gap lengths are adjusted so that each 


pin carries an equal amount of current. 
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3.4 ELECTROMAGNET DESIGN 


An integral part of MGD discharge stabilization is the design 
of a proper magnetic field. Such a field is essential in producing 
the Lorentz driving force and the magnetic bottle used for containment. 
Various magnet designs were analyzed and tested by using a two 
dimensional iterative computer program called POTENT, which was devel- 
oped at Culham laboratory [3.4]. Inputs to the POTENT program included 
the specific overall dimensions of the magnets, iron core permeabilities, 
coil current densities and boundary conditions. The use of this pro- 
gram enabled a detailed evaluation of magnetic field and flux density 
components which were illustrated and plotted. Various parameters 


could easily be altered and the results quickly viewed on a CRT screen. 


3.4.1 Cylindrical Magnet 


To satisfy the necessary field requirements in the planar 
electrode geometry, a cylindrical re-entrant electromagnet was chosen 
as the most promising design. A cross sectional view showing the 
magnetic core and the resulting flux lines is reproduced in Fig. 3.7. 
Axial symmetry has been assumed so that the magnetic flux density is 
independent of 6. As illustrated in Fig. 3.8, the contours of magnetic 
field intensity reveal a significant non-linearity in field strength. 
This change in intensity creates the strong shear forces that are set 
‘up between the electrodes. 

The magnet core was constructed of mild steel, anodized to 


create an insulating layer on the surface. The magnetic field was 
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Fig. 3.7 Magnetic Flux Density of a Re-entrant Electromagnet. 
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Fig. 3.8 Magnetic Field Intensity Contours of a Re-entrant Electromagnet. 
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created by passing a current of 10 A through 3400 turns of number 18 
gauge magnet wire. However, in order to maintain this high current 
density within the winding, small gaps were built into the coil to 
facilitate the flow of water or oil. By water-cooling the magnet, the 
effective amp-turns could be increased and the excess heat generated 
could be easily removed. These considerations enabled the magnet to 
yield a flux density exceeding 0.2 T. The entire interior assembly was 
then sealed to allow the magnet to work within the discharge chamber. 
The computer program POTENT can also be used to determine the 
magnitude of the two magnetic field components B, and Boe. In this 
geometry, the current density vector points in the g-direction. The 
addition of a radial magnetic field B, creates a 6-directed I x B, 
force exerted on all charged particles. The axial seGnahe field B, is 
used in connection with the Hall current J to produce a radial force 
wich tends to force the gas discharge into an annular region. This 
Lorentz force sets up a pressure potential wel] that acts aS a con- 
finement bottle to prevent the discharge from blowing apart. For 
maximum effect, the electromagnet was placed directly behind the 
cathode plate to allow the strongest magnetic field to permeate the 


most active zones of the gas discharge. 
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3.4.2 Coaxial Magnet 


To test the radial discharge geometry, a short solenoidal magnet 
was chosen to provide the proper field configuration. As displayed in 
Fig. 3.9, the solenoid was short enough to produce both radial (B,) and 
axial (B.) field components. As expected, the contours of magnetic 
field intensity were symmetrical due to the imposed o-symmetry of the 
computer program. The results are clearly seen in Fig. 3.10. 

The coaxial magnetic field was created by passing a current of 
5 A through 6500 turns of number 20 gauge magnet wire. The magnetic 
field within the solenoid was enhanced approximately 45% by the intro- 
duction of a 1 cm thick steel shell. The additional steel intensified 
the magnetic flux, in the center, by redirecting flux lines towards a 
path of least reluctance. The bending of these flux lines can be seen 
in the computer plot of the solenoidal magnetic field profile (Fig. 3.9). 

In this discharge configuration, the axial magnetic field B. in 
combination with the radial current J, creates a Lorentz force oe x zB ; 
which causes plasma rotation in the negative 6 direction. At the 
magnet edges, the radial field component is most prominent due to the 
presence of fringing fields. This field component B, >» interacting 
with the Hall current cy » produces a Lorentz force ih x B » which 
serves to confine the discharge axially to the central portion of the 


electrodes. 
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Fig. 3.9 Magnetic Flux Density of a Solenoid Electromagnet. 
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Fig. 3.10 Magnetic Field Intensity Contours of a Solenoid Electromagnet. 
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3.5 THYRATRON PULSER 


It has been demonstrated [1.5-1.9] that gas discharge stability 
can be improved by any technique which externally controls the electron 
density. This procedure allows the ionization mechanism to be decoupled 
from the gas discharge current transport process. In an attempt to 
control and supply the electron density necessary for this experimental 
setup, a high power pulse generator was developed. While the pulser 
was in operation, a uniform preionized glow discharge was established 
containing an estimated electron density of 10®8 electrons/cm°. 

The primary element employed for the repetitive high voltage 
Switching, necessary in the pulser system, was a hydrogen filled ceramic 
HY-1802 thyratron. A thyratron was chosen over several other devices 
because of its reliable high speed switching action and lower production 
of electrical noise. In addition, hydrogen thyratrons have a very short 
recovery or deionization time. In this context, ta is the time interval 
from the end of the thyratron firing until the grid has regained control 
for the next trigger pulse. For this reason, hydrogen thyratrons are 
well suited.for high frequency applications. 

As a result of hydrogen's low molecular mass, the singly charged 
ion, accelerated within the thyratron, cannot gain sufficient kinetic 
energy to inflict significant damage to the cathode. Consequently, the 
hydrogen thyratron can be driven much harder and will therefore be able 
to handle higher voltages and currents as compared with other devices. 

In order to activate the thyratron, or start conduction, it is 


first necessary to draw current between the grid and the cathode. This 
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current is initiated by the application of a positive voltage trigger 
pulse to the grid. Following the trigger pulse, ions and electrons are 
produced within the grid-cathode region of the thyratron. As soon as 
a sufficient electron density is established, the anode field is able 
to produce ionization in the grid-anode region. This excessive ioniza- 
tion causes an avalanche breakdown to occur between the cathode and the 
anode. The breakdown coalesces into an arc and the thyratron is fired. 
The pulse generating circuitry, as displayed in Fig. 3.11, 
consists of a high efficiency resonant charging curcuit. When the 
thyratron is non-conducting (between trigger pulses) the impulse storage 
capacitor Cy is charged via the charging inductor L, and diode D. The 
holdoff diode is used to prevent any reverse current flow. After the 
thyratron fires, the energy stored in capacitor C, is transferred to the 
anode through the coupling capacitor Cy. 
The average current transferred to the anode by the pulser can 


be determined from the approximate expression given by [3.5, 3.6], 


T= Vf. 


In this equation, Y is the capacitor charging voltage, C, is the forall 
circuit capacitance and ne is the pulse repetition rate. For an applied 
voltage of 5 kV operating at 500 Hz, the average current is typically 
50 mA. 


w9pptyd apes toy witieeg fie) 


*6 ceortsels bad anot eee Yap 

ee TOG2 BA MONT RIV SS one te ane ite “BID 

elds at bist! shone sis bade tate an 

palnot syldesoxs etn . ‘iinet spotis-bhp ont of pers | 
ony “Die siddes 3 ont neawied beg a cwabinsid osmat eve ne me 


sap : 
a} ie? 
h 


re de 


bert? ef agresveds ofd bas ey NER. aan? escies ten, Poni 7 
TEE ar? 1 wT Ete 465 aitfopnis writ vianae nie 

ott erw wSPUINGD vier aa dnsnape eno aad 
wenje Balam at (aaehuged saping: neawted) ‘* 
ont 0 bol’ bhs. 1 vofaubat patients old at bape, #2 
ait (eTTA. aor? Tayo SevSyes. Lek snovert of beau ai 


* «s é 
| A ae ae, ac 


gid) od Det shana a LS nawh 25065 Ye ‘rote ene eo 


pee a an ‘ 
i ‘ 


| na bottawon as 
2 - in hee : 
nea seelud snd yd shone oft Of inerrio. sparevs oft 
ah a oy fh 
Lot (2,8) ye hewitp not datate sreinote oft 


ary 
oe i A 7 * 


iv bl ay f be N : ie ‘ se 
; + ihe ‘ah ea oe ake pth a 
° - - 


“J 


faded add 2t 9 sRttow wat tag oid of W. tial ine 


i, Aes 


bellaqs We wed Jae nets Thegee ee ‘ai as “bonds, songs oan 3 


: heats J 


A 
Uisstey a eri mayan . ee pe eee w we Sa i. | 
Ae tvs Li oa d 


| : 


= ry 0 - 
| si iv i 


59 


Ly D Cy Co Ds | 
pen 000007. Y 
DC Pulser DC 
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e = C3 
_ HY - 1802 > > 
Cathode 
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SinWre a 


HY - 1802, Hydrogen filled ceramic thyratron 
L,, Charging inductor = 13.2 H 

D, Charging diode, 20 kV P.I.V. 

C,, Impulse storage capacitor = 1 nF 

C5, Coupling capacitor = 0.01 uF 

L., Bypass recharge inductor = 375 pH 

L3, Isolation inductor = 460 wH 

C., Isolation capacitor = 0.1 uF 


Fig. 3.11 Hydrogen Thyratron Pulse Generating Circuit. 


hl a 


PS eer Leeeeees o0n2- ce aan 


a 
vi inguc . 
aeons } 
pai es 
¥ 
‘ 
Pe ad 


| & — 


abomise 


a a ss 


: — mages (3 Hoa ‘ . 
ys phy vim ban mn 


an hat age 


U 


“s1uSt79. pnt toreoaagst naraerNiT aByoNli «17-8 att 


*~ 
. 


60 


5.0. DC. POWER-SUPPLY 


The complete power supply circuitry is shown in Fig. 3.12. 

From this illustration, it is seen that power was drawn from the mains 
and fed into a 3 phase motor controlled variac. In this manner, the 
stepping motor would adjust the amount of power to be sent to a 3 phase 
50 kVA transformer. The 60 Hz signal was then rectified by a high 
current. diode bridge and filtered by a 375 uF capacitor to remove the 
remaining AC ripple. This setup produced a maximum current of 25 A 
that could be fed to the gas discharge. 

Since the gas discharge was well-behaved and operated in the 
normal glow region, less current limiting resistance was required. 
When operating the system, using the current ballasted anode, no 
external ballast resistance was necessary. As previously shown in 
Fig. 3.6, the total fluid resistance was only slightly more than 25 Q. 
This was fortunate since more energy could be deposited into the gas 
discharge resulting in a more efficient device. 

To isolate the pulser signal from the DC supply and voltmeter, 
a 460 uH inductor was placed in the main line. Any pulser signals 
penetrating the inductor were simply shunted to ground via a 0.1 uF 
high voltage capacitor placed near the experiment. High voltage co- 
axial cables, with the outer sheath conductor grounded, were used 
throughout the system. This was necessary in order to prevent the 
wires from radiating like antennae during the operation of the pulser. 

The IV characteristics of the gas discharge were continually 
monitored by two analog meters. In order to maintain a solid ground 


connection, a high current 2.5 ma shunt was utilized instead of passing . 
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the current directly through an ammeter. This procedure eliminated 
excessive line inductance which would impede the conduction of the 
pulser signal to ground. A voltmeter was connected directly across the 
incoming power line and ground. This resulted in a measurement of the 


voltage drop across the internal ballast fluid plus the drop across the 


gas discharge impedance. Thus, it was necessary to make a separate IV . 


measurement of the fluid ballasted anode. Subtraction of these two curves 


yielded the true IV characteristics of the gas discharge. 

An arc is characterized by a very low impedance and consequently 
draws a large amount of current. Thus, to provide protection for 
CUNCULC components — in the event of a glow-to-arc transition, a vacuum 
relay was installed, so that the operator could instantaneously break 
the DC circuit. In addition, to prevent damage to the power supply 
and diode bridge, several 25 A fuses were connected an series with each 
line. Utilizing these precautions, greatly reduced the potential 


damage that could be caused by an arc. 
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CHAPTER 4 
MGD MODEL OF A GLOW DISCHARGE 
4.1 INTRODUCTION 


In order to properly simulate the time dependent behavior of a 
rotating glow discharge, it was necessary to construct a consistent 
set of equations. The theoretical and numerical analysis used to pre- 
dict the expected flow pattern and gas discharge position was based 
upon the compressible MGD equations. However, a complete description 
of this laser discharge system must also include: an electrical con- 
ductivity which exhibits strong spatial dependence, spatially varying 
electric and magnetic fields, and the presence of significant radial 
and axial secondary flows. Furthermore, a typical glow discharge 
contains ionized and neutral particles, each following a complicated 
trajectory that is influenced by a variety of forces. To simplify the 
analysis of such a problem, a macroscopic fluid model is used to 
describe the average behavior of a collection of particles in lieu of 
an individual particle. 

A macroscopic theory can be applied to the analysis of ordinary 
gases, provided the mean free path is small compared with any natural 
length over which the gas properties may significantly change. One 
further criterion required for the validity of this theory is that the 
mean collision time be short compared with any relevant variation in 
the fluid. 

For a typical glow discharge, molecular collisions will prevent 


any molecule from straying too far from its neighbours. Thus, the 
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finite fluid element is provided with a collective identity in terms 
of the specific molecules it contains. In addition, molecular 
collisions produce a distribution that is locally Maxwellian and can 
therefore be described in terms of a local thermodynamic state. 

In this fluid approximation, the partially ionized gas is 
assumed to be composed of three interpenetrating fluids: a positively 
charged ion fluid, a negatively charged electron fluid, and a neutral 
particle fluid. The neutral fluid will only interact with the charged 
particles through collisions. The ion and electron fluids, however, 
can also interact as a result of self-generated electric and magnetic 
fields. 

To reduce the complexity of the problem, the three separate 
fluids will be combined into a single fluid. The MGD equations 
generated from this less detailed approach will be later implemented 


into the computer program MAGIC. 


4.2 SINGLE FLUID MGD EQUATIONS 


The behavior of a flowing ionized gas can be described by a 
system of conservative MGD equations written for each individual gas 
Species. These conservative equations reveal the manner in which the 
average properties of the flowing gas are conserved by the particles 
during a collision or transport process. The conservation of mass 
equation and the momentum conservation equation for each gas species can 
be obtained from the velocity moments of the Boltzmann equation. As shown 


in Tanenbaum [4.1], the mass conservation law may be expressed as 
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d 
ent ed iy (4.1) 


The global conservation of mass can be attained by summing 


Eqn. 4.1 over all species. The overall density becomes p = ) Po? and 
s 


the mass averaged velocity is ov = , Pag Thus, the global conserva- 
§ 


tion of mass or equation of continuity becomes 


=e Poy aan \e ea0ie: (4.2) 


0 
s 3 _ ‘ aoe 
vs = 0 since mass is conserved in collisions 
@ 


The collision term ) 
s 


involving chemical reactions. 

The term Ov, in the continuity equation, expresses the rate of 
flow (mass/unit area) of the ionized gas. If the divergence of this 
term is positive, the gas is expanding; and if the divergence is 
negative, the gas is being compressed. Since a glow discharge is 
primarily composed of neutral particles, little change in overall] 
mass density is expected. However, a significant change in the mass 
density profiles of the charged particles is anticipated. 

The conservation of momentum equation for charged particles can 
be written as 


> >> > * > >> 63 
— + Ve : + Vell = + +1) § 4.3 
3b? 5%s) : (00,4) z ie ae ne2,\# 5) ( ) 


Similarly, the conservation of momentum equation for neutral particles 
is given by 
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The term I, represents the partial pressure tensor of the s species. In 


these momentum equations, the collision term is expressed in the form 


ae) 
ae 2 eee tm. ‘vgi’s"’p (4.5) 


> 
where A. is the species momentum transfer between r and s particles. 
* 
The variable Ving |S known as the average collision frequency for 


momentum-exchange between r and s particles, and is written as 


‘ns mm m_D i (4.6) 


The variable oes represents the mutual diffusion coefficient between 
ry and s particles. Lastly, the variable i is defined as the diffusion 


or drift velocity with respect to the mass average velocity, that is, 


8 eee > 2 
ve nee eo (4.7) 


The last term in Eqn. 4.3 represents the driving force which sets the 


ionized particles into motion. The expression 


B, = 02, (#08) | (4.8) 


is more commonly referred to as the Lorentz force. 
The overall momentum equation is obtained by summation over al] 


species. The term ) a = 0, since momentum is conserved during each 
Ss 


collision. Thus, the sum of Eqns. 4.3 and 4.4 result in the formulation 
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The first term on the right hand side contains the net charge density, 


given by 


P, = ) Neo : 
s 


Since a normal glow discharge is primarily composed of an electrically 
neutral positive column, belo O. This assumption is not valid within 
the charged particle sheaths located in the cathode fall layer. The 
complete analysis of the cathode fall region requires the use of a 
particle code and is presently beyond the scope of this thesis. 
| The second term on the right hand side of Eqn. 4.9, contains 


the total electric current density which is written as 
> > 
J = L n,eZ v - (4.10) 


The total current density can be separated into a conduction current 


and a convection current by Eqn. 4.7. That is, 


oh 
iT 


) 77 pte See OFS j 
L neZV, +v L Micon Gt Up (4.11) 


As determined earlier, ais O within the positive column. Thus, the 
total current density is predominantly due to the conduction current. 
Using the expression for the mass averaged velocity, Eqn. 4.9 


can be simplified into the form 


200) + ¥e(ovv) + Women ns (4.12) 


To account for the transport of momentum and energy through the 
gas by the diffusion of particles, the total pressure tensor I is 
written in terms of a viscous stress tensor ¢ and the hydrostatic 


pressure p, 
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Il = x + pL es ag 
where I is the unit tensor. The viscous stress tensor is created as a 
result of internal friction produced when different fluid elements move 


past one another at different velocities. According to Landau and 


' Lifshitz [4.3], t may be represented in subscript notation as, 


Ov» oN ”, ov, 
tg ane: le a 3 (zu -)6,, dey, (4.14) 


where u and z are the shear and bulk viscosity coefficients. In vector 


notation this equation becomes 


aha T ; 
T = - 1 {95+(%5) } ts (Su -c) (F-3)z (4.15) 
where (¥)7 is the transpose of the dyadic Vo. 
To close the system of equations, the thermodynamic equation of 
state is applied. This equation relates the pressure p to the mass 


density o. Since a glow discharge can be treated as an ideal gas, the 


equation of state has the form 


D Cot (4.16) 
where C is a constant and y is the ratio of the specific heats, 
eee Canes ‘ 
To simplify the problem, each species within the gas discharge 
is assumed to have a constant temperature. In this isothermal case, 
y = 1, and the overall gas pressure is simply the sum of the individual 
partial gas pressures. Consequently, 


p=) n kr, 
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where k is the Boltzmann constant. Since the degree of ionization is 
approximately iOnae nys<che and the total hydrostatic pressure is 


primarily due to the neutral gas molecules. As a result, 
pn kr (417) 


where qr is the translational temperature of the neutral gas. 


4.3 CURRENT DENSITY DISTRIBUTION 


For a partially ionized gas, Sutton [4.2] derived an expression 
for the generalized Ohm's law in terms of the applied electric and 
magnetic fields and possible electron pressure gradients. This equa- 


tion can be written as, 


Se 2 * 
Vp >> an Te : 
= —— =. ——— - i + 
Af oii haere rsa (2 a) Vp, xB Bx (<8) |} (4.18) 
e e 2 4E 
where a s 
n € n € . 
oe) CUTS ve MALES T eS re Se . (4.19) 
mul ty. 3} e 
e en er 


is the scalar electrical conductivity parallel to the magnetic field 
vector. In addition, C is the mean collision time derived by the 
first approximation of the Chapman-Enskog method. However, due to 

its complexity the more accessible Lorentz model will be used through- 
out this thesis. It must be noted that both these models become 
equivalent when the particle distribution function is Maxwellian. The 
variable of is the electric field vector in a coordinate system moving 


with the mass average gas velocity, that is, 
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>*k > > > 
E 


MEE USOD tN (4.20) 


The symbol f represents the mass fraction of gas molecules which are 


not ionized. Mathematically, this can be stated by the expression 


For a partially ionized glow discharge n >> 1; 5 SO that =". 
Lastly, the symbol a represents the ratio of ion to neutral particle 
mass. Hence, a =m./m = 1. 
(2 Ids 

The first term on the right hand side of Eqn. 4.18 is the result 
of an applied electric field; the second term accounts for any electron 
pressure gradients, the third term represents the Hall effect, and the 
last term accounts for ion slip. 


For a partially ionized gas, influenced by a small magnetic 


field, the ion slip term in Eqn. 4.18 can be neglected if 


* * 
Rese ar °F Ube The < 
Ce e@e CoCr 


In this expression, Wee and wep are the electron and ion cyclotron 
frequencies, respectively. As will be shown in Section 4.5, this 
inequality will indeed hold true for the glow discharge under investiga- 
tion. If the electron pressure gradient is also neglected, Eqn. 4.18 


can be simplified to 


* 
te ok Yoleg + 3 
J = a = B 5.3 ° (421) 


When the vector terms in Eqn. 4.21 are expanded into component form 


[4.4], the scalar electrical conductivity is transformed into the 
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tensor o. This procedure produces the well known expression for 
Ohm's law 
F2G+R = G+ (BB) . (4.22) 

Thus, the applied driving force is the result of the combined electric 
and magnetic fields acting upon the charged particles. 

In the transverse geometry under investigation, the magnetic 
field consists of two components, namely BA? s8) and B(r.z). In 
order to obtain the electrical conductivity tensor, within the non- 
uniform magnetic field region, a similarity transformation, into 
cylindrical coordinates, must be performed. As derived in Appendix A, 


Gc is transformed into the matrix expression 


° 2) tf ®) - ° in ° x 
o, Sin?¢ +o, ,CcOos*¢ oSin ¢ Mees g, )sin $+ COS ¢ 
G = o7sin ¢ o, -o7 COS 6 |e (4.23) 
2 2 
-o )sin $° cos S 524 t+ j 
(o,, ? o° COs ¢ O7COS > 9, COS*H +o, SiNM¢ 


> 
where on, is the scalar electrical conductivity parallel to B, while 
o, and or are the Pederson and Hall conductivities, respectively. The 
argument ¢ represents the angle of rotation between the localized and 


general coordinate systems. This angle can be expressed by the relation 
ty -] 
¢ = tan (B_/B,,) 


In order to obtain an exact solution for the current density 
within the gas discharge, it is necessary to simultaneously solve for 
the electric field and the rotational velocities using the expression 


for co. Such a technique was employed by Wijnakker et. al. [1.21] who 
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used the simplifying assumptions of a constant electrical conductivity 
in a nonvarying magnetic field. However, for the problem under investi- 
gation, the electrical conductivity is very nonlinear, particularly in 
the cathode fall region. Consequently, the conservation of charge 
equation will be used to produce an approximate solution for the current 
distribution within the discharge. 

Assuming a steady-state condition exists within the glow dis- 


charge, the conservation of charge equation implies that 


ded 


=F pace 0 i: aay j 
Ved = el) reper ths (4.24) 
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Within a transverse electrode geometry, the electric field is primarily 


along the axial direction, so that E>? E Thus, it would be reason- 


Ror 
able to expect that the axial current density J, would dominate the 
radial current density J . Evaluation of Ohm's law (Eqn. G72) 
employing the conductivity tensor in cylindrical coordinates (Eqn. 4.23) 
allows a comparison to be made. As shown in [4.5], the relationship 


between J, and J, is given by 


uy) B2tan2 

eee e 
J. 1+(1+62)tan26 
J. ay hate tan2¢ 


* 
rae is the Hall parameter. Examination of the magnetic 


where 8 = w 
e 
field components B, and B. ,» uSing the computer program POTENT, has 
revealed that over the annular discharge region, B/B., << |; Con- 
sequently, for the geometry considered, ¢ 1s small. Estimates of the 
Hall parameter have indicated that B, ~ ]. Thus, under these conditions 


a << 2. and Eqn. 4.24 simplifies to 
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Since the electrical conductivity is a function of the unknown, 
Spatially dependent, electron density, it cannot be properly evaluated 
using the one fluid MGD approach. It is hoped that this problem can 
be resolved at a later date, through the use of a two fluid model. 

As a first approximation, a current density distribution which 
matches the discharge characteristics will be selected. It has 
frequently been observed that the gas discharge was forced to rotate 
in an annular region by a combination of radial Lorentz forces. This 
annular region was seen to encompass the entire glow discharge includ- 
ing the cathode fall region. Due to this nonuniform electron emission 
from the cathode, a Gaussian distribution function was chosen to 
represent the actual current density within the gas discharge. 

In order that the Gaussian distribution function accurately 
represent the gas discharge, it was important that this function be 
centered on the region of highest observed plasma intensity. Using 
this guideline, the distribution function can be formulated by the 


expression P 
(evn) 


Af (eS SS ee 
Pai ee (4.25) 


where b* is the variance and lhe is a constant. The width of the annular 
shaped discharge, d can be defined as the distance between the points 
on the distribution function curve where the current density falls to 


one-half its maximum value. These points are located at r = i alZs, 
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where r, is the point at which the peak current density occurs. Using 


0 
this standard, the variance of the curve can then be related to the 


annular discharge width in the following manner: 


Simplifying this relation and taking the natural logarithm of both 
Sides produces a simple expression in which the variance is directly 


proportional to the discharge width, 


Pa 
parade or (4.26) 


The constant J, can now be determined by integrating the total 


applied current over the discharge area, that is, 


(r-r,) 


Ori phe e -— 
nes | J_(v)dA = is eam West tentis (4.27) 
a 0 40 bv 2n 


Evaluation of this integral yields the expression 


AjiaSe lew, vent ac 
+ 7Pr er er 
| Bp : 


where erf(x) is the error function normally written as 


2 -¢2 
(x) = 2f e Ghee 
Ssult ‘tho jee! 


The results of Eqns. 4.26 and 4.28 can now be substituted back into 


(4.28) 


the Gaussian distribution function described by Eqn. 4.25. By using 
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data taken from high speed still photographs, the corresponding 
distribution function has been calculated and plotted as illustrated 


Mets 4 als 


4.4 TRANSPORT PROPERTIES OF IONIZED GASES 


The phenomena of viscosity, thermal conductivity and diffusion 
in nonuniform gases represent effects which tend to make the mass 
velocity, temperature and composition aniecen through collisional 
processes. Transport phenomena involves the convection of these gas 
related properties due to a gradient or disturbance within the medium. 
These transport processes are vital when computing the expected macro- 
scopic behavior of an ionized gas for a given set of boundary 
conditions. However, before any analysis is performed, the various 
transport coefficients must be identified for each gas component. 
Values for these transport coefficients can be estimated once a suit- 
able model has been chosen. 

In a weakly ionized gas discharge, the rigid sphere model 
assumes that charged particles are predominantly deflected by single 
collisions with neutral particles rather than by multiple scattering 
with other charged particles. This assumption is based upon a charged 
to neutral particle ratio of n/n, ~ 10’. With the application of 
an electric field, the electrons and ions now move in the direction 
of the field. As a result the distribution function is no longer 


perfectly Maxwellian, but slightly anisotropic. The ions, because of 


their greater mass, are assumed to move very slowly compared with the 
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Gaussian Shaped Current Density Distribution. 
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light, energetic electrons, The accelerated electrons, however, wil] 
eventually reach a finite drift velocity determined by the frequency 
of collisions with the neutral particles. 
In general, for particles with a Maxwellian distribution function, 
the average collision frequency between rv and s particles is given by 
1 
BSkE SKE 
r S 


fy 
TM 1m” 
L 


Via ri 


(4.29) 


S 
or in terms of the mean thermal speed 


<y >=n@Q (<e >+<ce >). 
HS: ‘Sy, ALS 1 s 


In these expressions, ane is the collision cross section between r and 


s particles [4.2]. 


4.4.1 Viscosity 


Viscosity is a process involving the transport of momentum 
through a moving medium due to a gradient in the velocity of that 
medium. The resulting momentum flux creates a shear stress in the 
direction of the velocity gradient. The coefficient of shear viscosity 
can be defined as 


abe 
Bde wouen’ hme: 
The coefficient of shear viscosity can then be defined as 


i 5 Eat oe (4.30) 


Substituting expressions for the mean speed and the mean free path 


produces the result 
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for collisions with the yr species. The viscosity for a gas mixture is 


simply the sum of contributions from the individual gases, that is 
= 2 eo edad (4.32) 
ue ah er - 
s 


where wah is the mutual viscosity of the sth and rth gases [4.6]. 

It can be readily seen that the coefficient of shear viscosity 
is independent of its density. Hence, the capacity of a gas for 
transmitting momentum is not changed when its density is altered. 
Consequently, an increase in pressure will not normally retard the 
motion of an ionized gas flowing between two electrodes. However, if 
the viscous gas is moving at a high velocity, ordinary laminar motion 
along irregular boundaries can become unstable, and turbulent motion 
will begin to dominate. In this instance, momentum is conducted away 
by turbulent eddies and not by ordinary viscosity. Therefore, the flow 
resistance of a high speed turbulent gas may depend upon its density. 

This problem was encountered during early experimentation with 
a water-cooled tubular cathode. The turbulence produced when the 
azimuthal ionized gas passed over this electrode created hot, 
stagnant gas cells which led to the formation of instability streamers. 
This problem was subsequently solved by the substitution of a smooth 
cathode plate. 

Another form of viscosity present in most gases is the bulk 
viscosity z. This second type of viscosity is concerned with the 


internal state of the gas molecules. As described by Hirschfelder 
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et al. [4.7] the coefficient of bulk viscosity is given by 


_ n*k2r 2) 
n= C2 ) ¢, ee 7 (4.33) 


The parameter t, is the characteristic relaxation time required for the 


Q 
transfer of energy from the translational to the internal degree of 
freedom for the 2th mode. The term o\®) is the contribution of the 
particular degree of freedom to the specific heat per molecule Cy ; 


The heat capacity at a constant volume can be written as 


C= 5nk (4.34) 


for an ideal gas. For dilute monatomic gases c¢ ~ 0, and for poly- 
atomic gases, zc is generally small and depends upon the relaxation 
times. These characteristic times are normally the same order of 


magnitude as the collision interval ae 


4.4.2 Electrical Conductivity 


For an electrically neutral glow discharge, the total current 
density carried by charged particles is primarily due to the conduction 
current. Thus, for electrons and positively charged ions, the current 
density is 


> > > 
J = ne(V.-V,) ; (4535) 


However, in a glow discharge, v << 7, » and the current due to the 
ion drift can be neglected. The drift or diffusion velocity of the 


electrons can be calculated from 


> > 
V. = <e >=— | ec f de 
e e n, eiew € 


(4.36). 
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=e . 
where e, 1s the random thermal velocity vector and Aes is the electron 
distribution function in the presence of a magnetic field. As derived 
by Sutton and Sherman [4.2], the current density can be written in the 
tensor form 


* 
ais = ) oe 
de ad 


(4.37) 


where On. is the conductivity tensor. During evaluation of Eqn. 4.37, 
it may be shown that the conductivity tensor can be written as the 


Single integral 


2 oe fomeiil de 
a ee? UJ e 0 € (4 38) 
Lda ok Wes ; 


where the collision frequency tensor 2; is defined for a magnetic 


field in the zg direction as 


es 5 ce 0 
ce es ce es 
“ee Yes 
ce es ce es 
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For a general magnetic field, De must be modified by a similarity 
transformation as described in Appendix A. 

The components of ay can be attained by integrating Eqn. 4.38 
with the approximation that ke is to be considered at some average 
value vee Since ere is independent of Css it may be taken outside 
the integral sign. Simple integration produces a first order approxi- 


mation to the conductivity tensor, that is, 
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She cage lt (4.40) 


Due to the form of the collision frequency tensor, only three distinct 


conductivity components remain: 


ne* PR le 12° Ww, 
yy Crepe yes Fe TS Cgettay 2) Spray, (w2 #<y >2) ° (4.41) 
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The curvature of an electron moving under the influence of a 
magnetic field, is responsible for the transverse or Hall current. 
Rewriting the transverse conductivity in terms of the mean collision 
time produces 


Cea 
ee es 


(4.42) 
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This equation can be regarded as the ratio of the Hall current to the 
current parallel to the electric field. The significance of this ratio 
is determined by the Hall parameter given as 8 =w tt. _ Thus for 

e ce es 
Bek 1, the Hall current can be of major importance to the operational 


characteristics of a glow discharge. 


The Hall parameter can also be written as 


(4.43) 


where Me is the electron mean free path and Pr, is the average Larmor 
radius. From this result, the Hall parameter can be viewed as the 


number of gyrations an electron makes between collisions. 
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4.5 DETERMINATION OF THE TRANSPORT COEFFICIENTS 


In order to analytically or numerically predict the behavior 
of the rotating glow discharge, an accurate estimate of the transport 
coefficients is required. However, under the influence of a magnetic 
field, these coefficients become tensor variables. Thus, in order to 
evaluate these coefficients, the magnetic field must be known through- 
out the discharge volume. Magnetic field values were computed by using 
the computer program POTENT. 

A direct calculation of the electron density is not possible 
through the use of a single fluid MGD treatment. For this reason, the 
total current density is used to estimate this parameter. For neg- 


ligible ion current, the steady-state averaged electron density is 


given by 
e§ 
a ui 
LL rears | (4.44) 
e 
where In is the total discharge current distributed over an area 4. 


The electron drift velocity within a composite gas mixture of He/N,/CO, 
can be determined through the analysis of Elliott et. al. [4.8]. 
However, their calculation does not include the effects of a magnetic 
field. As a result of the discussion in Chapter 7, the ratio F/n must 
now be modified to include the presence of a magnetic field. In order 
to satisfy this condition, the Hall parameter Be Aeon os must be known. 
The electron cyclotron frequency is easily determined, but the mean 


collision time is dependent upon the unknown electron temperature. 


Without a direct measurement of the glow discharge electron temperature, 
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this evaluation cannot be performed except through an iterative 
procedure. 

On the first step of the iteration, the ratio F/n = 2.64~x 
10° ®V/cm* is calculated for a 16/2/2 torr laser gas mixture with no 
magnetic field. Reference [4.8] then provides an estimate for the 
electron drift velocity and electron temperature. These values are 
used in turn to calculate the electron density (Eqn. 4.44) and the 
electron-neutral collision frequency (Eqn. 4.29). The Hall parameter 


is then computed so that a new E, /n ratio can be determined. After 


nd 


several iteration cycles are completed, Z,../n = 2.44x 10 !©V/cm? for 


ff 
a magnetic flux density of 0.04 T. This translates into an overal] 
E/n reduction of 8 percent throughout the glow discharge. Values of 
magnetic flux density, near the cathode, typically reach 0.1 T. In 
this region, the £/n ratio can be decreased by over 65 percent compared 
to its non-magnetized state. Consequently, a small change in magnetic 
field can have a substantial effect on the glow discharge. 

The gas discharge typically operates at a current denial of 
50 mA/cm? . Through the use of Eqn. 4.44, the average electron density 


£3 


is estimated to be Ont 1.2x 10¢!cm™? for an average drift velocity of 


V. = 2.6~ 10°cm/s. Thus, with a neutral gas density of 


n 4.54 10°°cm™ >, at a temperature of 425 K, the degree of ioniza- 


n 
tion is n/n, = 2.64% 10°. Subsequently, charged particle collisions 
can be neglected. 

During the iteration procedure the average electron temperature 


was also calculated. For the above discharge conditions <é> converged 


to 1.3 eV. 


83 


Li a, 
fee 


RR ae cot Bet — | 


et vo aya rae ne neo a 
ose feet ey aeoni Lor uecani icotts 
ott bet ce ga) iat’ ails otizuet a on deladiel 

se Seinen ae { gH edit 4 27 a Laaesanst ty : 
(ae. baneitera tet 4 won.a Mt 


oF VPN ip ee 2 : 


: ag oe 
ae wa ett oe Sly suo saan Phe 00 


Bi) A" 


TPievaye: sh 
8 gourTey 


fi] af i 


\urt eat dentate aprons oe er) ey 5 oi tee. > 
to Wo sgihey fab -apntavn rip wr tik tb ps ‘t: ty mae ot | mn aa 
*@ var ena Pah Bp ae - 


et vous — 
nmenot * sens awd, uN 2. A ‘rats » Bi, nal ‘or oi te : 
moieties otattiing baw | | a oa 


‘4 r 


| snyterecqins Aoto0! & pevone i eunsore notteaeat ait ena a 


bagrpinin <3 anos eae apruane _ Tauad 1% abies eta 


The average molecular mass of a gas mixture, can be determined 


by the expression 


pee 


(4.45) 


“nie ) N, 
where NV, is Avogadro's number, M, is the molecular weight, and f; is 
the fractional percentage of the 7th gas. For the particular gas mix 
under consideration, the average neutral or ion mass is 
M tm = 1-726 10° °kg/particle. 

To determine the momentum transfer collision frequency for 
neutral particles, ions, and electrons, it is first necessary to know 
the collisional cross section of each species. The electron-neutral 


collision cross sections for an ionized gas mixture with an average 


electron energy of 1.3 eV, are displayed in Table 4.1. 


10°27 °m*) 


ae 


Table 4.1] Electron Impact Cross Sections for Elastic and Inelastic 
Processes when <é> = 1.3 eV. 


As shown earlier, the collision cross section of two gas 
Species is 


2 
= Pay e 
Qs 4 n(x. rs) 


When considering a gas mixture, the average collision cross section 


for momentum transfer becomes the summation of terms 
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By using Eqn. 4.46 and the molecular radii in Table 4.1, the average 


cross section for the specified gas mixture is a A = 21.5x10 om. 


mtx 
Since ions have a low kinetic energy within the glow discharge, the 


collision cross sections for ion-neutral and ion-ion interactions can 


be assumed to be Qe ~ Qe © aan . Ina similar procedure, the average 


electron-neutral particle collision cross section is determined by the 


Summation 
oma ase (4.47) 
mis Jd 
Using the values listed in Table 4.1, @, = 7.264 teu ae 
mtx 


The average momentum transfer collision frequency can now be 
determined for each species in a partially ionized gas. Summing 


Eqn. 4.29 over all particle species results in the equation 


Spe Sut aut er, (4.48) 
jh 


Thus, the average collision frequency for neutral particles can be 


written as 


Cpe =Wo kl hy oo (4.49) 
nN NL ne nn 


However, since Nin, << nm, > collision frequencies involving charged 
particles are significantly less than those involved with neutral 
particles. Hence, <by Hpeecn os With similar modifications, 

Se Pie and poe ear Ue 


As stated earlier, the electron energy has been estimated to be 


kr, ~ 1.3 eV. Moreover, due to the large number of collisions in the 
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positive column, the ion temperature is relatively low. Thus, a good 
approximation to use is r. os zr where the measured neutral gas temper- 
ature is ite = 425 K. However, the assumption does not hold in the 
cathode glow layer where the ions can accumulate a significant amount 
of directed energy as they accelerate through the cathode fall. 

Now, using the calculated average molecular mass and the average 
collisional cross sections for the gas mixture, the species collision 
frequencies are evaluated by Eqn. 4.29 and 4.48. Substitution of these 


parameters yield 
Sicusae <yetsrs 1.285 fe108S7) (4.50) 


Deal erage) ura 


A 
Si 
Vv 

il 


For an average magnetic flux density of 0.04 T, the electron and ion 
cyclotron frequencies are nen 7.034 10? rad/s and 


ie On i.x 10°rad/s, respectively. Through these results, 


; 4 F : ; : 
Cen Ol re oe 8.2x 10 , which affirms the earlier assumption of 


neglecting ion slip. 
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4.6 ORTHOGONAL CURVILINEAR COORDINATES 


To maximize the versatility of the computer model, the MGD 
equations are expanded in orthogonal curvilinear coordinates. This 
procedure would allow the program to accept different electrode 
geometries using different coordinate systems. In order to utilize 
this powerful technique, the force fields and differential operators 
must be expressed in terms of a generalized coordinate system. Since 
the differential operators are not often found in this form, they will 
be listed in this section for convenience. 

In this particular system, a differential arc length may be 
described by 

Che asta el oe ede at pode aal, 


where the term h de. respresents the differential arc length in the 

e, direction. The quantity h, is known as the scale factor for the 

x ,-coordinate and is generally spatially dependent. Associated with the 
scale factors is a transformation matrix known as the metrix tensor G. 


The determinant of this matrix is given by 


2 
Gael 25 Meic aloles 


In a generalized coordinate system, special notation is required 
when dealing with the dot and cross products of the unit vectors e, , 


These operations are described in the following manner: 


rm le, “TOR = 7 
Cyne Se : (4.51) 
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> > 
Cake. 


< a 
ee eee (4.52) 


where c.., is the Levi-Civita permutation symbol. This symbol has been 


tk 
defined as 


+], for an even permutation of (zjgk) 
=4-]  , for an odd permutation of (zjk) 
O , otherwise. 


It now becomes necessary to describe expressions for various 
operations of scalar, vector, and tensor quantities in a generalized 


coordinate system. In this system, the gradient operator becomes 
> > 
e eo 
P) Pe Nee re 
apelt Smbaees ) fe ee (4.53) 


for zt = 1,2,3. By using this result, the gradient of a scalar ¢ is 


do2. 
Vo = ) a oe (4.54) 
ees L 


ey Lg | . (4.55) 


a .with-elebees (4.56) 


The expression for the divergence of a tensor Tf is given by 


ee) ee BONG) 
i baton Ds } aS (4.57) 


where re is the covariant derivative of 7. By using expressions 
i) 


developed by Morse and Feshbach [4.9], Eqn. 4.57 is transformed into 
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(F-%), = h a Facer (ygr"F) + J ati] (4.58) 


where 7! = T../n.n. is the contravariant form of the tensor 7... The 
ta Cog 7) 


expression eo is a Christoffel symbol with the following properties: 


: -h. oh. - 
4] = tied 3 Ge Nees on ee ° 
l. 4 Ti)? 3a, 5 \, . Oy" GOVar.7 seattedar erent. 
This notation, for the Christoffel symbols of the second kind, repre- 
sents the change in direction of the unit vectors as one moves in 
Space. For completeness, the first component of Ve? expands into the 


equation 


> » Ee ] aie) 3 3 
OS ie rg steht if ee 1) see ota 


3h, oh, oh oh. 
ZF hT > Ba, *P2t13 32, hT,> cE ak 3a, : (4.59) 


This equation contains terms which look like the normal vector 
divergence, plus four extra Coriolis terms that arise due to the change 
in the direction of the unit vectors with changes in position. 

The final relation required for the compressible MGD equations, 
is an expression for the gradient of a vector. As developed by [4.9], 


the resulting dyadic is given by 
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With these definitions it is now possible to write out the com- 
pressible MGD equations in component form. First, the generalized form 


of the conservation of mass, Eqn. 4,2, becomes 
re) re) re) ( 
e+ aly hov,) + mae iis! 2) *! Co hyh,ov,)} = 0. (4.61a) 


Next, the generalized form of the conservation of momentum, Eqn. 4.12 


is given as 
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The dyadic terms Ass are described by Eqn. 4.60. The two other components 


of the momentum equation can be obtained by interchanging all subscripts 
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CHAPTER 5 
COMPUTATIONAL TECHNIQUES 
5.1 INTRODUCTION 


The understanding of complex gas discharge systems has, thus 
far, been limited to very specialized applications. Although gas 
discharges can be analytically described, these studies generally 
represent an oversimplified view of the actual situation. Moreover, 
detailed realistic theoretical studies of the behavior of a gas 
discharge with electromagnetic fields (both static and time-varying) 
are indeed rare. The reasons for this lack of understanding is due 
to the extremely complex nature of such systems. A typical gas 
discharge may consist of electrons, fully and partially ionized ions, 
and neutrals all rapidly varying with time. Consequently, only 
computer simulation techniques are adequately suited to handle the 
complexity involved in these types of problems. 

A computer code has been developed to solve the nonlinear 
time-varying, two-dimensional coupled set of compressible MGD 
equations describing a partially ionized gas discharge. The complete 
description of such a plasma is impractical due to the enormous 
number of particles, their complex trajectories, and the resulting 
induced electromagnetic fields. This problem can be circumvented by 
the application of statistical mechanics to describe the plasma in 
terms of distribution functions. These averaged plasma properties 


are then used to compose a fluid description of the plasma. In this 
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manner, a macroscopic view of a rotating gas discharge under the 


influence of a nonlinear magnetic field can be developed. 


5.2 BASIC NUMERICAL CONCEPTS 


In the numerical solution of a system of differential equations, 
the continuous Space-time domain is replaced by a discrete space-time 
mesh. In addition, the differential equations are approximated, at 
each mesh point, by an algebraic finite difference equation which 


: to the known variables 


relates the dependent variables at time ae 
at time 27 If the finite difference equations are uncoupled in 
these unknowns between time levels, then the equations are classified 
as explicit. Therefore, in order to solve a set of m difference 
equations valid over a two-dimensional lattice containing WV x uM points, 
a set of mxNxM algebraic equations must be solved for mxN xu 
unknowns. 

Should the finite difference equations be coupled with respect 
to the unknown variables at the time level 4”, the resulting mxW xm 
implicit equations must be solved simultaneously. As stated in [5.1], 
the number of multiplications involved in solving Y simultaneous linear 
algebraic equations is approximately om, Therefore, the computing 
time per time step in a two-dimensional fully implicit calculation can 
be quite excessive. Thus, for multi-dimensional problems, this type 
of calculation is rarely used. 

In order to avoid this limitation, an alternating direction 
implicit (ADI) method may be employed. As derived by Peaceman and 


Rachford [5.2], this algorithm breaks the time advancement of the 
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equations into two steps. In the first half step, variables are 
advanced implicitly in one dimension and explicitly in the other 
dimension. During the second half step, the procedure is reversed. 
Consequently, in the ADI technique, there are alternately y uncoupled 
sets of mx WN coupled simultaneous equations combined with W uncoupled 
sets of mx M coupled simultaneous equations to be solved. This 
technique drastically reduces the computational costs in contrast to 
the fully implicit method. The mesh points involved in each of the 
explicit, fully implicit, and ADI techniques for second order operations 
in two dimensions are depicted in Fig. 5.1. 

The numerical solution of the multi-dimensional compressible 
Navier-Stokes equation requires a considerable amount of computer time, 
and consequently, efficient computational methods are necessary. A 
number of prior methods for solving this problem have been based upon 
explicit finite difference schemes. However, this procedure is 
Subject to stability restrictions on the time step sizes as dictated 
by the Courant-Friedrichs-Lewy (CFL) condition. This condition states 
that the difference solution, which can propagate information across 
the mesh at a maximum rate of one space step Ax during each time step 
At, must exceed all propagation speeds of the physical system. This 
includes such diverse effects as waves, fluid motion and diffusion [5.3]. 

Within a general viscous gas discharge, influenced by time- 
varying magnetic fields, an explicit difference scheme will encounter 
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Difference Points Involved At Time 
Technique ¢ int 


Explicit, 2D 
Uncoupled 


Implicit, 2D 
Coupled 


ADI, Odd: 
Number Steps 


ADI, Even 
Number Steps 


Fig. 5.1 Two-Dimensional Mesh Point Orientation for Several Finite 
Difference Techniques. 
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Through this relation, the CFL condition then ensures that any physical 
velocity, within the system, will be less than the lattice speed, 
Axc/At. In Eqn. 5.1, v is the advective velocity, u is the viscosity 


coefficient, and oe is the adiabatic sound speed given by 
<A ens (5.2) 
S 9 
The term v, is designated as the Alfven speed 


= « (5-3) 


In general, discharge regions penetrated by a high magnetic flux tend 
to have low densities as a result of a balance between kinetic and 
magnetic pressures [5.4]. This effect can create very high Alfven 
velocities which severely restrict the allowable explicit time step. 
Therefore, the significance of this stability restriction is that 
computational efficiency is lowered by imposing a smaller time step 
than would otherwise be desirable. 

In contrast to explicit schemes, the implicit finite difference 
technique is stable regardless of the time step size for linear 
differential equations with constant coefficients. In this case, the 
values w*** are fully coupled to the values 2 at the old time level. 
Thus, a disturbance at any grid point can propagate to any other point 
on the nesadining One time step. As a result, the lattice speed or 
velocity of information transmittal on the mesh is infinite and the 
CFL limitation is always Satisfied. 

For the implicit case under consideration where the differential 


equations are nonlinear, stability theory suggests that the time step 
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restriction is relaxed from the CFL condition [5.5]. This enhanced 
numerical stability should create a substantial increase in computa- 
tional efficiency since fewer time steps are required. 

The magnetic Reynolds number Re is a measure of the ease with 
which the magnetic field can slip through a fluid. More precisely, 
this fluid interaction number is the ratio of fluid flux to magnetic 


diffusivity or 


iE e e 
vL diffusion 
kR = Se, (5.4) 
m  Wugo Tproblem 


Here, v is the advective velocity, Z is the scale length or mesh 
increment, is is the free space permeability, and o is the fluid 
electrical conductivity. Within the slightly ionized gas discharge 
under investigation, i 2n x 10°’ and even moderate size currents 
will only produce very slight perturbations in the applied field. 
Therefore, the magnetic field equation can be neglected in the analysis 
of the flow pattern. With the absence of this HUTT! the critical 


velocity used in the CFL condition becomes the acoustic velocity. 


For an operational pressure of 30 torr and a neutral gas density 


of 0.0256 kg/m’, the acoustic speed, e, = 480 m/s. This translates to 
a maximum explicit Courant time step of approximately 10 us. Asa 
resultof this small time step size, it is hoped that the relaxed time 
step limitation of the ADI method will substantially reduce computa- 
tional cost. 

The choice of an implicit scheme has also been found to be 
beneficial for situations involving high Mach number flows. Briley 


and McDonald [5.6] have found that in these cases, the CFL condition 
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becomes increasingly restrictive when a small lattice spacing is 
desired. The subsequent use of an implicit technique reduces com- 


putation time while providing the necessary spatial resolution. 


5.3 COMPUTATIONAL MESH 


Since the mathematics dealt with in a computer is both discrete 
and finite, the values of a continuous function must be defined as 
discrete points in space over a finite mesh. This discretization 
technique also arises in many areas of physics and EM theory where 
the equations of force and motion are applied to small finite elements. 
In the limit, these discrete elements become vanishingly small, 
thereby producing a differential equation. 

To represent a continuous function over the domain (Xx, »X,)> Tt 
is first necessary to divide the domain into a set of W - 1 elements 
of width AG. . A vector set {a .} may then be constructed by 
defining the continuous variable x only at the points z. Thus, the 


Space lattice points are 


where be, = 2, - 2 The dependent function f(x) will now be 


t t-1 * 
approximated by its corresponding vector set {f, 3 defined on the 


independent variable mesh {x} such that f; = f(x.). The nature of 


L 
this approximation must be such that the scale length Ax» be small in 
comparison to any wavelength of the rapidly varying function f. Thus, 
the mesh element Ax essentially defines a cutoff wavelength, below 


which no phenomena can be described. 
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For a general initial value problem, the time dimension is also 
divided into small finite intervals which separates discrete levels or 
moments in time. This process can be represented by 

x . 
t= 40+ J vty 


ves 


ml Here again, at” 


where ¢2 is some initial time and At” = 4” - ¢ 
must be sufficiently small in order that significant time dependent 
information not be lost. 

The continuous function u(x,t) can now be expressed as the 
discrete function ue = u(ae, st”) defined over a temporal and spatial 
mesh. The computational mesh designated for the program MAGIC is 
classified as Eulerian since it has no velocity with respect to a 
stationary observer in the laboratory reference frame. 

The values of the discrete function ae) are specified to 
be cell-centered quantities. In this description, the value of bs 
is the average value of uz in the two-dimensional cell centered at 
aa The two-dimensional mesh is illustrated by Fig. 5.2 where 
Prewors represent the points at which the values of the state vector 
u are stored. Similarly, the crosses represent the interface points 


at which the fluxes and forces are to be evaluated. 


To adequately describe the computational mesh, the mesh elements 


ee 


n 


A.. and Ae, ane defined as the distance between points x and x 
ices z t+] 


where as {x sx} are members of the orthogonal curvilinear coordin- 


ates. That is, 


ee det] (5.5) 
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Fig. 5.2 Schematic of the Finite Difference Computational Mesh. 


i—1 i+1 


Fig. 5.3 Illustration of the Mesh Points Involved in a Mixed 
Derivative Calculation. 
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These variable mesh elements allow the specification of a nonuniform 
mesh. Therefore, the distance between cell interfaces, denoted by yd> 


is given by 


(570) 


The value of any variable at the cell interface is simply the 


average of two adjacent cells. For example, 
C7) 


Finally, to determine the flux or force difference across an interface, 

the 6 Operator produces the terms 

SpA Blas orl ny cy) 
1+ A tt1 , 


Lt 
i (5.8a) 


peel 
6, _(F) = BAe. Pea) : 


These operators represent the forwards and backwards differencing scheme 
along the a direction, which are evaluated at the cell interface. 


Similarly, for the x. direction 


(5.8b) 
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5.4 SPATIAL FINITE DIFFERENCING 


The conservative MGD equations (4.61) are composed of two basic 
types of terms. The first type arises from the divergence of flux 
transport terms as a result of advection and viscous processes. From 


the expanded MGD equations, all flux terms have the form —(aF) where 
. 1 


a iS normally a transport coefficient and F is any general factor 
including differential terms. The second type of term found in these 


equations have the form 8“ where 8 is some function of the state 


vector u. These force terms Originate from 5 x B and pressure gradient 
forces. The resulting flux and force terms are now used to denote the 
_ type of finite differencing that will be applied throughout the MGD 
equations. 

The application of a spatial finite difference technique to a 
set of differential equations, forms a system of algebraic equations 
which depict variations to the cell-centered variables. These changes 
are due to the transport of flux through each individual cell as a 
result of appried driving forces. Hence, flux and force terms are 
calculated at the cell interfaces. 

Since flux terms measure the flow through a computational cell, 
the finite difference form evaluated at the cell boundaries, is written 
as 

Orr iete 


3 
an) [ae ae are (5.9) 


1 


where a, and F, are simply the average of two adjacent cells as defined 
by Eqn. 5.7. If the variable F represents a derivative term, the 


finite difference form becomes 
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where the 6 operator is the difference across a cell interface. 
For the frequently occurring case of mixed derivatives, 


F = 2, it is first necessary to take the derivative in the x, 
2 


direction at the points 3,,2_ and average them to arrive at a value 
Pi at the cell interface. The procedure is repeated at the points 
A,,4_ to determine F_ . The differential with respect to x, is then 


- performed to give the expression 


it fe, S ie, 
+ ox ={92 
af mo 2 AT. 2 


P) aG . 

fun octet Pals Eb) or ay 
1 e l 
where 
GORE an erp 1 
Ear Zee, seat se iD ¥ De ee see 
j (5.12) 
1 : | 
‘ We fee g) ; De es Be Par ely wen} ; 


The points involved in this derivation are shown in Fig. 5.3. 

The force terms measure the overall effect that applied forces 
have on the edges of a computational cell. In this case, it is 
necessary to define these terms as the average of the forces applied 


on each cell face, that is, 


| | 
scan: paid co ace 5 all (5.13) 
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5.5 TEMPORAL FINITE DIFFERENCING 


The system of equations derived in the previous chapter can be 


represented by the vector equation 


yf > >--> >> >> > > 
st”) +: F(u) teG(u) + H(u) = 7 x B = 0 (5.14) 


where uw = Ue 2h sD is the state vector of the system, and 


>> > 
T(u) = {o¥,p} are the quantities advanced with time. The function 

>> 

F(u) contains all terms in Eqn. 4.61 that have leading derivatives with 
respect to wo. Several of the possible forms that these leading 


a: a aR CMe 07 do. ; 
derivatives can take include: Ae ; ae 14) and Be 2) foneatt) 7. 


In a similar manner, the functions G(u) and H(u) contain all terms with 
leading derivatives with respect to x, and x, respectively. However, 
due to the type of applied driving force, the rotating plasma is 
assumed to have 6 symmetry. As a result, all partial derivatives with 
respect to 6 are zero for the case under consideration. Therefore, 
the function G(u) = 0, 

) The ADI finite difference method, discussed by Peaceman and 
Rachford [5.2] allows for the solution of large matrix problems using 
a straightforward algorithm. This technique employs two finite dif- 
ference equations which are used in turn over successive time steps 
-each of duration At/2. Each part is treated implicitly on alternate 
time steps while the other part is treated explicitly. Thus, if u* is 
an intermediate value at the end of the first time step, then this 


time-centering ADI technique yields a two-step approach: 
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> > 
* Ss 
Pee Ds PG) + BG") - 7 x B= 0 (5. 15a) 


and 


Bee pees Oi (5.15b) 


These two equations describe the time advancement scheme employed 
in MAGIC. Initially, Eqn. 5.15a is used to solve for the intermediate 
values es which are then applied in Eqn. 5.15b to determine the 
solution uot at the end of the time step At. It can be shown that 
these equations are each accurate only to first order in At, but when 
the two equations are combined, a second order accurate equation is 
produced. Consequently, the overall ADI mroeeduve is second order 
accurate. 

The forms of the functions F(u) and H(u) are well defined 
except for mixed and product derivative terms. The application of 
finite difference techniques to these terms requires advanced time 
values for rows not yet computed by the tridiagonal ADI method. In 
order to accommodate these terms, banded matrix techniques are 
normally used. To avoid these complicated banded techniques and 
still use a tridiagonal ADI finite difference scheme, it becomes 
necessary to treat the mixed derivatives in a special way. Finan [5.5] 
Suggests several involved methods, but advocates treating the mixed 
derivative terms in an explicit manner. This procedure lowers the 
accuracy for these terms to order At and introduces a Courant time 
limit on the time step size. However, it was felt that the simplicity 
of this mixed derivative treatment would outweigh the occurrence of any 


Slight inaccuracies. The problem of product derivatives could be 
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treated in a similar manner. These terms are normally found in the 


energy conservation equation. 


5.6 LINEARIZATION OF THE FINITE DIFFERENCE EQUATIONS 


Applying the temporal finite difference techniques to the MGD 


equations, yields a set of nonlinear algebraic equations of the form 


> > 


Tue) eae y + = Ol Orig Ska fe ted Oat Fae SAF xB) = 0 (Se16) 


> > oe > er A: Sera 
where u = {vo} and 7(u) = {pv,0}. The term pv in the vector 7 is an 
example of a nonlinear term. Since Eqn. 5.16 is nonlinear, terms such 


as ov are present when both p and vw are being sought for at the next 


time step. Consequently, a straightforward solution of this system of © 


equations is not possible. Instead, an iterative technique will be 
developed to solve for the nonlinearly coupled unknowns in Eqn. 5.16 
at the new time step. 

The nonlinear set of vector equations given by Eqn. 5.16 can 


be written in the form 


>> > > > + > 
Wu" ur) = W(u") + W.(uu*) = 0, (5.17) 


> oy re voyn > 
Since W(u,u ) can be considered to have an explicit part Wu ) which 


a4) 


ce ey + oy ook 
only depends on uw; and an implicit part W (wu su) which depends on u 


> 
and 7". 


In order to linearize this set of equations, the function 


Wea must be continuous and differentiable so that it may be 


> > 
expanded in a Taylor series. In this manner, the function WO) 


be expressed in terms of the function and its derivatives at yi 
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as follows: 


ey) 2 


ee, (2 (2 | At + O(at)~ , (5.18) 
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where W=W and 2 is the iteration index. The linearization 


is performed by a two-step process of expansion about the known time 


>*k 
level. ¢” and Subsequent approximation of the quantity eat by 


alia Bos - this term arises from chain rule differentiation. If the 
vector ay is close to the actual root of Eqn. 5.18, then ay = 0. 
In addition, if. a is very close to ae then all terms of second 
(iy 


order or higher may be neglected, that is = Q. Applying 


these results, the Taylor series expansion simplifies to 


ead 
0 = H+ Ba ME) + ont)? . (5.19) 
ou 


The matrix (9i7/9u*) is a standard Jacobian whose elements are 


defined as 
> 
5 t Thewunec a(W7, 4H, 3°"7,.) 
ae Eras ed 
yn \V) 2V, 90490 


Since the Jacobian J only measures the ue dependence, only the implicit 
: ares eae ; : 
portion of W(u ,u ) need be differentiated. 
The expression in Eqn. 5.19 describes a generalized Newton- 
Raphson iteration method which provides second order convergence. 


Rearranging this equation reveals a simple iteration formula, namely, 


a +] - 
Bon Te ea ae oY (5.20) 


Ad 


where the first iteration iy takes on the value of wu” at the old level 


t”. As derived by Finan [5.5], this equation is representative of the 
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actual expression to be solved. The solution of the finite difference 
expression (Eqn. 5.16), at the grid point z, will involve state vector 
values at three lattice points, namely 7 + 1, ~, andz- 1. Asa result 
of the Newton-Raphson iterative procedure, a Jacobian must also be 
determined with respect to each of these lattice points. Thus, for an 
implicit expression in the zth direction, the finite difference scheme 


produces a formulation similar to Eqn. 5.20, that is, 


Qe Q+1> Qe Q+1> Qe Q+1> 
é Ned Pie Ger md ® (i ate 9 f 28 Z J 
Gt leg i Lad (bea Bee (AG | t~ 153 
27 {> QF > ae > Q 
= dle e = UR e aL oe e bd GRE e + e © 2 Use e - W ° BCA 
Ls 7 tn lig tod tod ee sie Lal oJ ( ) 


The function a also includes the explicit mixed derivative terms. 
The solution matrix can now be rearranged into the standard ADI 


form, namely 
a nse . a a 1 4 
sa Pere, aoe =p | (E22) 


iG = “t4ang td tug tg t-lsg tad 
where ne B, and C are the Jacobians differentiated with respect to the 
vector uw on the lattice points (7+1,7), (t,7), and (t-1,7) respectively. 
The vector D contains all the explicit terms and satisfies the 
equation 
Ree hei eee i Chane ome. §—7(5123)) 
UJ UJ tt1 59 LJ LJ oJ Lass 
Since the ADI vector equations (5.22, 5.23) are uncoupled with respect 
to 7, the two-dimensional problem can be expressed as two problems in 
which only one dimension is advanced implicitly per half time step. 
More simply, Eqn. 5.22 takes on a block tridiagonal matrix structure 


where each block is an m x m submatrix. Here, m represents the number 
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of state variables to be solved, that is, three velocity components 


and the mass density. 


written in the form 


By “C) 
-A, B, ao 0 
-A, B, ~, 
-A Be -C 
L £ Lt 
0 Taye By ee 
“Ay By 


Sy &Y S&F 


The resulting tridiagonal matrix can then be 


Sone) 


In this representation, the block elements are depicted by the Jacobian 


matrices A, B, and C. This convenient form is then readily solved by 


a Gaussian elimination method. 
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5.7 SOLUTION OF THE DIFFERENCE EQUATIONS 


The problem of solving the set of two-dimensional partial 
differential equations has now been reduced to solving sets of simul- 
taneous linear algebraic equations. The resulting difference equations, 


after the first half time step, have the form 


aA nw 


SA eMail CERNE peo Cle. SEDs (5.25a) 
tod Pach Ws) VJ tod tod t= lag to 
and on the following half time step 
bir >t 1 “* is nt 1 = “xk je nt] = Ds (5 25b) 
ta ee y (Se) SOT) ea POEL. De) : . 


For a set of m partial differential equations, A, B, and C represent 
m x m matrices while u and D represent m component vectors defined 
at the mesh point 7,7. Equations 5.25 are valid over the interior 
mesh points for which z = 2,3,...,W-1 and j = 2,3,...M-1. However, 
they do not form a complete set. The remaining equations are provided 
by the boundary conditions at z equal to 1 and W for all jg; and for gj 
equal to 1 and ™ for all 2 

The ADI finite difference scheme requires boundary conditions, 


on the first half time step, to have the form 


GE a Toe Ty eee (5.26a) 


and 
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and 


a nan 


+k +k 
©. e = e° = e ar e e 
FEN las i) anna a a Vs Is oy (5.276) 


for 1 <7z< WN. In the above equations, G and H are m x m matrices 

and # represents an m component vector. Each of these four boundary 
conditions satisfies one side of the two-dimensional mesh. On the 
second half time step, the boundary conditions will have the same 

form as Eqns. 5.26 and 5.27 but with the superscripts m and * replaced 
by * and + 1, respectively. 

The procedure for solving this system of equations, is a 
generalization of the method developed by Richtmyer and Morton [5.7] 
for a two-dimensional case. In this situation, their scalar equation 
becomes a vector equation with matrix coefficients. Since Eqns. 5.25a 
and 5.26 are uncoupled with respect to g, the state vector Us along 
lines of constant 7, can be determined independent of U along other 
lines of constant 7. This attractive feature of the ADI technique 
leads to the recursion relationship between the u values at adjacent 


mesh points, that is, 


(5.28) 


The values for the matrix Ey 


de 


calculated by first rearranging the boundary condition (Eqn. 5.26a) 


. and vector P, j at z = 1 can be 
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into the form 
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Comparing this result with Eqn. 5.28 yields 
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To obtain the recurrence relation for z j and ig j for 


1 <t<WM, it is necessary to substitute Eqn. 5.28 for z - 1 into 


Eqn. 5.25a. This produces the expression 
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Comparison of Eqn. 5.30 with 5.28 shows that 
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Valid? forec ="2,35...6n- 1. 
Lastly, the recursion relation of Eqn. 5.28 with 7 = WN - 1, 
can be substituted into the boundary condition (Eqn. 5.26b) to solve 


for ly al athatcis.. 


“wn 


* oY % >* > 
ia > ( ) 


= e + . 
Nig ON .4 4 ee one eaten Tn 


Rearranging this expression yields 
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The solution algorithm is now almost complete. During calculation, the 
recurrence matrix B and vector F are first calculated by making a 
forward sweep from z equal to 1 to W - 1 using Eqns. 5.29 and 5.31. 
Next UN ¢ is determined through Eqn. 5.32. Finally, the remaining 
State vectors uF are computed by backsubstitution using Eqn. 5.28. 

_ The method presented, thus far, enables the calculation of all 
u* values for all ¢ when 1 < j <M. The remaining values along 7 = 1 
and g = M for 1 < zt < NW can be calculated by rearranging the boundary 
condition (Eqn. 5.27) using previously determined values. Now, only 
the u* values located at the four corner POMMtS Gls l)so( 1M). (N.1).; 
and (V,”) remain to be calculated. Any symmetry condition along either 
boundary can be used to estimate these points. If no symmetry con- 
dition exists, then the values at the corner points can be determined 
by passing a polynomial through the adjacent mesh points, for which 
the u* values have already been calculated [5.8]. A second order 


accurate polynomial for the corner point (1,1) produces 


> ke oe Hk 
u =u +i) -Uu : (5.33) 
ial Ik Dranl! 292 


The remaining corner points are estimated in the same manner. 
On the second half time step, the difference equations are 


solved along rows of constant z, using the equation 
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The procedure followed when the zth component becomes the implicit 
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variable is identical to the previous half time step, but with the 
Subscript indices z and 7 interchanged. 

For a general set of m x M x N simultaneous equations, an 
m x M x NW matrix must be inverted. However, the method presented 
takes advantage of the block tridiagonal nature of the ADI finite 
difference technique. As a result, this procedure only involves the 
inversion of Mm x m matrices along one dimension and W m x m matrices 
along the other. The substantial reduction in the number of floating 
point operations in the ADI scheme greatly enhances computational 


efficiency and reduces roundoff error. 


5.8 BOUNDARY CONDITIONS 


To ensure a proper solution over the entire mesh, the boundary 
conditions must be accurately specified. In addition, the boundary 
conditions must be compatible with the block tridiagonal form of the 
ADI solution technique. Examination of the matrix expression in 
Eqn. 5.24, reveals that the boundary conditions are represented by 
the rows ¢ = 1 and z¢ =wW for the case when z is implicit. These rows 
have the general form 


Cores nis Ga (5235) 


where uw... is the boundary point and mi is the next interior point. 


BP INT 
All forms of Dirichlet, Neumann, and mixed boundary conditions can 


be depicted by the equation 
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where = is the derivative normal to the boundary. Applying the 
finite difference technique to Eqn. 5.36 transforms it into the 
Standard form of Eqn. 5.35. Consequently, all classifications of 
boundary conditions can be employed in the ADI scheme. 

In order to preserve second order finite difference accuracy 
throughout the program, ghost point boundary conditions were selected. 
In this method, the boundary walls are chosen to lie halfway between 
the first two mesh row points and the last two rows. Through this 
procedure, the derivatives along the boundary are second order 
accurate approximations. 

The simulation mesh covers the domain 0 < r < no anda Qu< 2° <-c. 


Due to azimuthal symmetry, the solution is valid for all values of 6. 


-The anode and cathode can be assumed to be perfectly conducting walls 


which lie along the rows g = 1+ 1/2 and 7 = M- 1/2. Thus, the ghost 


or exterior mesh points lie along the rows 7 = 1 and 7 =m. This 
particular arrangement of boundaries is illustrated by Fig. 5.4. 


Since the electrodes can be considered a solid wall, 


where n is the unit normal vector outward from the wall, and D is 
the velocity at the wall. Since B,, is parallel to n, this Dirichlet 
boundary condition can only be satisfied if ee 5 O. Therefore, at 


“the electrode walls 


pole,0)) =v. (25d) = 0 C53) 


a a 


The boundary conditions must also consider that the ionized gas cannot 
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Fig. 5.4 Arrangement of Ghost Point Boundaries on the Finite Difference 
Mesh. 
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Slip along these walls. Consequently, the tangential velocity By is 


fixed at the wall by the relation 


With this Dirichlet condition, 


v (250) = v,,(v 4) = 0 
and (5.38) 
v(x 50) = v, (rd) =O. 


In order for the solution to remain finite, an axis of symmetry 
exists along ry = 0. Furthermore, the following Dirichlet conditions 
must hold: 


v,(052) = v (0,2) = v, (0.2) =N0". (5.39) 


Along the outside radial boundary, viscous drag from the walls will 
start to reduce the azimuthal velocity. However, even though the 
driving force terms are now reduced to zero, strong viscous coupling 


of the fluid will keep oe high. Consequently, at r =R the 


Oey 
azimuthal velocity is subject to a Neumann boundary condition of the 
form 


Ah ve(R, »#)) =0. (5.40) 


The boundaries, for the secondary flow velocities, are sufficiently 
removed so that the solution should approach zero. Therefore, they 


satisfy the Dirichlet condition 


v R92) = v, (Ry 2) ale, (5.41) 
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Outside the domain of interest, the neutral mass density is 
assumed to be constant. Therefore, to satisfy this condition, the 
boundary condition 


—— =.(0 (5.42) 


is applied along all boundaries. 

Applying the finite difference operations of Section 5.4, to 
Eqns. 5.37-5.41, produces a set of equations which will represent the 
boundary conditions. Combining these difference equations into the 
matrix Eqn. 5.35, creates the submatrix elements used in the block 


tridiagonal form. Along the g = 1, and Jj = mM boundaries 


hi 0a ke 
‘ A Ons lsak creek 
G., = G. = (5.43a) 
PTE Rey aes Ton 
Ore One Ce) 
and 
-T 0 0-0 
2 2 Sie 0.0 : 
H. =H. = (5.43b) 
Di tee 0 see 0 
0 Oi? se] 
and 
7. =7.=0 (5.43c) 
De ates! te .43¢ 


where the subscripts UY and D refer to the ‘up' and ‘down' boundary 
walls. 


Similarly, for the t= 1, and z = W boundaries 
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] 0 0 0 
- 0 ] 0 , 0 > 0 
: = * R = N-1 (5.44a) 
0 ] 0 0 ] 0 
07a 10 ] 0 0 ] 
and 
-] 0 0 -] 0 
" 0 =] 0 e 0 5 0 
L F s R = N (5.44b) 
GLO 2-1] 0 0 -] 0 
Oe. Oca) ig 0 0 ] 
and 
J = J, = 0 (5.44c) 


where the subscripts Z and R refer to the 'left' and 'right' boundary 
Sides. 


5.9 ITERATIONS AND DYNAMIC TIME STEP CONTROLS 


The derived set of nonlinear finite difference equations 
employs a Newton-Raphson iterative procedure in order to accurately 
invert the block tridiagonal matrix. Since the ADI technique decouples 
the rows of the mesh, it becomes possible to vary the number of 
iterations per row depending upon the degree of change in a state 
variable along that row. Thus, rows with large changes will: iterate 
many times while rows that have little or no change will only iterate 
a few times. Each row is tested for convergence by comparing the 
maximum variation of the state vector between iterations versus the 
average value of the iterations. Specifically, the convergence 


criterion e is given by 
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Q+1> Y> Q+1> Q> 
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This condition is subject to the constraint that should Mang be less 
than a preset floor value, then the floor value is used for 

by in the denominator of Eqn. 5.45. The maximum value of e« along 
a row is then compared with an allowable convergence criterion een 


fe c¥s oh thet the row is declared converged; otherwise the iteration 
procedure is continued. The value of the convergence criterion depends 


upon the accuracy desired; typically e < 1036 for most operational 


max 
runs. 

In order to minimize computational costs, a dynamic time step- 
ping routine was also implemented. This algorithm permitted the 
program to operate over a range of allowable time steps. The extent 
over which the time step could operate was determined by monitoring 
the maximum change in the state vector between time steps. Should the 
state vector values be below a floor condition, no check was performed. 
Furthermore, by constantly revising the floor condition to 5 percent 
of the maximum mesh velocity, the program time step will not be limited 


by small variations on the mesh. 


The maximum variation over the mesh can be expressed as 


a nt 1 on -n 
Vmax = max [lp site) ay | (5.46) 


To implement this dynamic time stepping procedure, the maximum 
variation US is compared with a prescribed maximum allowable change 
e. For a typical computer run, the maximum fractional change, in any 


state variable, has been kept to less than 60 percent per time step, 
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that is, 2 =0.0.. J1f ea < 0.8¢, then the time step is increased by 
10 percent. For the case ee > 1.5e, a recycle is performed; other- 


wise if 0.8e < La < 1.5e, the time step is decreased by the relation 


At = 0.8a¢(c/V_..) } (5.47) 


This technique is very flexible in that it allows the program to check 
the variation of any variable; including velocity, momentum, or energy. 
As a result of these time step controls, the optimal efficiency of 


MAGIC, with respect to computational running costs, is maintained. 
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CHAPTER 6 


EXPERIMENTAL RESULTS AND DISCHARGE OBSERVATIONS 


6.1 NUMERICAL RESULTS OF THE TRANSVERSE GEOMETRY 
6.1.1 Velocity Profile 


As described in the previous chapter, a computer code was 
developed to model the time dependent behavior of a gas discharge. 
Its primary objective was to determine the three-dimensional gas 
flow pattern created by the Spatially dependent J x B Lorentz force. 
In order to accurately simulate this behavior it became necessary to 
include secondary flow terms in the MGD equations. Such terms are 
often assumed to be small and are neglected in order to simplify the 
analysis in gas discharge systems [1.17-1.21]. However, as will be 
subsequently described, this assumption can be an extreme over- 
simplification to a complex gas discharge problem. 

An analysis of the steady-state magnetically stabilized gas 
discharge was initially performed by neglecting secondary axial and 
radial mass flows [4.5]. The rotational gas discharge velocity 
obtained from the one-fluid Navier-Stokes equation is illustrated in 
Fig. 6.1. This plot indicates that a strongly sheared velocity 
profile exists across the discharge. Furthermore, a peak rotational 
velocity of 250 m/s can be noted. Subsequent experiments with the 
transverse geometry have revealed a much smaller rotational velocity 


in addition to the observance of a radial afterglow. This large 
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Fig. 6.1 Rotational Velocity Distribution for the Transverse Geometry 
Neglecting Secondary Flow. Discharge parameters include: 
of j=) 22, 0A /CMss i= 10.60%. 10-2 skG/(ines ),, 


i 


i78\2) 


3 


u 
ai 


“A 


<< 
z 


ond fi 
jet I 


ee a *) 


1 


rea 
, 


discrepancy between the initial theoretical model and the experimental 
Studies can only be explained by the neglect of the secondary flow 
terms. 

In-a gas discharge system, designed for a high power laser, it 
is preferable to excite the gas molecules into the upper vibrational 
and rotational levels rather than produce ionization. As a result, 
the applied #/n values are kept relatively low so that the degree of 
ionization is maintained in the range 10° to 10’. Since the number 
of energetic ions is orders of magnitude less than that used in 
highly ionized gas centrifuge devices, substantially less momentum 
will be transferred to the predominantly neutral gas via inelastic 
collisions. As a consequence, greatly reduced gas flow velocities 
are expected in a laser gas discharge. 

Upon application of a magnetic field, the charged particles 
are Subjected to a Lorentz force which alters their normal collisional 
behavior. The derivation in Chapter 7 will show that the usual 
direction of the charged particle drift velocity will be changed such 
that the electrons and ions will have a longer inter-electrode transit 
time. During this transit period, the charged particles undergo 


additional inelastic collisions with the neutral particles. This 


process allows more momentum to be transferred to the neutral particles 


than would otherwise be possible. The additional momentum transferred 


to the neutral gas molecules is preferentially applied in the direction 


of the driving force, that is, the azimuthal direction. However, due 
to the existence of pressure gradient forces, ion momentum can be 


transferred to other directions as well. This process is responsible 
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for the creation of secondary flows. Within spatially dependent 
magnetic fields, additional Lorentz forces are produced which can 
further stimulate these flows. 

In this section, the results of the ADI program MAGIC will be 
presented along with an evaluation of its performance. As stated 
earlier, the rotational velocity produced by an applied Lorentz force 
1s expected to be significantly reduced due to the presence of 
secondary flows in this highly viscous, partially ionized gas. The 
plot in Fig. 6.2 is a representation of the azimuthal velocity Vp > 
4.0 ms after the Lorentz force has been applied. The parameters used 
in this calculation include a 2:2:16 torr gas mixture of C0,/N,/He 
with a neutral gas temperature of 425°K. mecdirions the total 
discharge current was 5 A and the viscosity coefficient, up was set 
to 6.6 x 10° kg/m-s. With these discharge parameters, the maximum 
rotational velocity was 35 m/s. This value is substantially below 
that predicted by the original model which neglected secondary flows. 
The peak rotational velocity occurs near the cathode surface since 
the driving forces are largest in this region. It is also evident 
that the rotational flow is strongly sheared throughout the discharge 
volume as a result of the spatially dependent Lorentz forces. The 
conservation of momentum equation (Eqn. 4.12) indicates that the 
rotational velocity will be proportional to the Lorentz force and 
hence to the discharge current. Thus, it appears feasible that 
larger flow velocities can be produced through the use of higher 


current densities. 
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Figure 6.3 1S a computer plot of the secondary flow pattern 
generated by the MGD model employed in MAGIC. The arrows show the 
direction and magnitude of the secondary flow at a time, 6.0 ms after 
the driving force has been applied. The longest arrows in this 
diagram represent velocities of approximately 32 m/s. This figure 
clearly shows the existence of a convective cell where the secondary 
flow is forced to recirculate upon itself. In order to save com- 
putational costs, the outer radial finite difference mesh boundary 
was moved closer to the discharge. The actual outer boundary wall is 
somewhat further away, however, little change in the recirulating 
flow pattern is expected. 

In this MGD system, the normally neglected secondary flows have 
been found to be approximately equal in magnitude to the principle 
rotational flow Ve: The existence of such large axial and radial 
gas flows must therefore be highly significant in determining the 
overall behavior of the gas discharge. It is obvious that such 
characteristics play an important role in the exceptionally effective 
Stabilization mechanism found in this system. It is suspected that 
any localized nonuniformities in density, temperature, or current 
are rapidly distributed over the entire discharge volume in a time 
very much less than the normal instability formation time. 

It may also be possible to exploit the large radial gas flows 
in order to cool the hot gas produced within the plasma. It is 
anticipated that with further adjustments, a large recirculating 
convective cell may be established. This would allow the hot neutral 


gas to be transported radially outwards towards a heat exchanger where 
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it can then be cooled and returned back to the discharge. In this 
radial flow system, the hot gas need only travel across one-half the 
electrode length to escape the excitation processes. Consequently, 
such a system will only have one-half the normal gas residence time 
found in conventional transverse devices. 

It must be noted that the calculated flow velocity values 
represent the mass averaged flow of the gas discharge system as 
determined by a one-fluid MGD model. Since the degree of ionization 
ts. onky-2.6.x Oa a: the plotted values are essentially those of the 
anne neutral gas and do not represent the flow pattern of the 
charged species. This additional information can only be attained 
through the use of a two-fluid model which could account for the 
motion of charged and neutral particles and their interaction through 
collisions. 

To ensure that the computer code was functioning properly, an 
elaborate series of checks were performed as outlined in Section 5.9. 
An additional accuracy test was also implemented whereby the state 
vector solution was backsubstituted into the original set of momentum 
conservation equations. Through this process it was determined that 
each equation was satisfied to within a relative error normally less 
than 3%. A detailed analysis of the program has shown that this 
iterative ADI code can normally be operated at twice the CFL time step 
limit without creating any significant errors. This result makes the 
ADI scheme very attractive since most explicit codes of this complexity 
can only run at a fractionof the Courant limit (usually 20%) just to 


ensure stability. The net result is a substantial savings in 
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computational costs. Further improvements in the accuracy and time 
step size are currently under investigation. 

It has been found that ADI direction splitting methods can 
lead to difficulties when large magnetic fields and low electrical 
conductivities are involved. Craxton and McCrory [6.1] have found 
that the direction splitting technique produced smooth solution pro- 
files in the last implicit direction scanned, but noisy profiles in 
the other direction scanned explicitly. This noise was also observed 
in MAGIC but could be eliminated by reducing the time step. 

At present the isothermal model cannot predict the expansion 
of the ionized gas due to a temperature rise. In the next phase of 
this project, MAGIC will be expanded to A an aneate the conservation 
of energy equation. The inclusion of this equation will be especially 
useful for modelling the heat transfer effect established by the 
rapidly recirculating flow. This will then complete the MGD model and 
Should produce an accurate representation of a magnetically stablized 


gas discharge. 


6.1.2 Radial Pressure Profile 


A thorough analysis of the rotating gas discharge system must 
consider the density distribution of the charged and neutral gas 
molecules. In order to determine these pressure profiles, a two- 
fluid model must now be employed. 

The momentum conservation equations for the charged and neutral 


fluids may be written as 
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aS 
dv 
ec@_t > > > Sor 
ee eras 2 a i a 
CAE TxB Te, p,vlv, v) UV20 a) 
ee Y - (3-0) + p V2 (6.2) 
en dt Py an VG ee : 


where d/dt is the convective derivative for each species, B, is the 


charged fluid velocity, 2 the neutral fluid velocity, Dees Da tip. 


is the sum of the electron and ion partial pressure, P, the neutral 

gas pressure, and ue and u, are the viscosity coefficients for the 

charged and neutral fluids. The parameter vw is the collision frequency 

for the transfer of momentum between the charged and the neutral fluids. 
The degree of coupling between the two fluids can be obtained by 


considering the steady-state 6-component of Eqn. 6.2, which yields 


pnsttys eltsauah a: (6.3) 
The Laplacian in Eqn. 6.3 has been approximated as -1/6%. An estimated 
value for 6 may be obtained from a second-order difference approxima- 
tion for the Laplacian evaluated at the midpoint of the discharge. 
This result is given by 6 =~ Z/4 where £Z is the anode to cathode 
distance (Z=0.0535 m). For the discharge under examination, 
u, = 6.6 x 10°, 6, = 2.07 x 10°, and the jon-neutral collision 
frequency is v = 1.285 x 1o8 The expected neutral velocity v6 
is therefore approximately 0.88 Der » which indicates a very strong 
coupling between the two fluids. It must be noted that the coupling 
that does exist is primarily between ions and neutrals since the heavy 


ions dominate the charged particle fluid equation. 
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A pressure profile for the neutral fluid can be determined by 
examining the steady-state r component of Eqn. 6.2 with bem svg) ) 
An ordering of terms will show that the viscous drag term is small. 
Thus, a neutral gas pressure profile will exist through a balance 


between centrifugal and radial pressure terms, that is, 


ey dp 
ns = 
eae per OR (6.4) 


This simple expression indicates that a strong centrifugal force term 
will transport neutral particles radially out of the discharge region. 
This result has indeed been confirmed by visual and experimental 
observations. 

Results from the computer simulations show only slight varia- 
tions in neutral mass density; although a small accumulation of mass 
was noticed near the axis of symmetry of the glow discharge. This 
effect can be explained by the small azimuthal velocity term found in 
the region near r = 0. Since the overall variation over the discharge 
volume is small, the neutral gas can be assumed incompressible. 

The pressure profile for the charged particles is more complex 
because of Lorentz forces tending to force the plasma into an annular 
region. A detailed analysis of these forces will be given in Section 
7.2. Since the neutral gas body is strongly coupled to the motion of 
the charged particles, On aganoea n° The radial component of the 
charged particle fluid equation, in steady-state, can now be approx- 


imated by the expression 
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where the largest convective term has been included. Through an 
ordering of terms, all remaining terms have been found to be small. 
The term JB. may be evaluated through the use of Ohm's law 
(Eqn. 4.22). The conductivity tensor c, has been expressed as a 
matrix given by Eqn. 4.23. From these equations, the current density 


in the azimuthal direction can be written as 


— he + = + - 
a 0150 B, o,,(v,B,, v Be) aI E(B vB) (6.6) 


where o; . refers to the tjth matrix element in the conductivity tensor. 


This equation may be simplified by neglecting terms of order oe 10; 2. 
a 
Therefore, to a good approximation 
Sp Gk tet Ps (6.7) 


Upon substitution of Eqn. 6.7 into Eqn. 6.5, the charged particle 


pressure distribution becomes 


2, 
"Pe Be eee 


=e a 0, COS($) EB, (6.8) 


where ¢ is the angle between the magnetic field lines and the general 


coordinate system. 


dp 
The plot of ree in Fig.6.4 was obtained by using the magnetic 


field profile for the re-entrant electromagnet shown in Fig. 3.7 and 
the predicted velocity profile obtained from Fig. 6.2. This figure 
illustrates the creation of a radial pressure well since B. changes 
direction and therefore sign. The pressure well depth varies directly 
with the magnetic field strength and is therefore greatest near the 


cathode surface. 
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The profile of the radial pressure gradient in Fig. 6.4 reveals 

that the Lorentz force is responsible for confining the gas discharge 
in one of two possible locations. For very low discharge current 
values, the gas discharge is located in the central portion of the 
pressure well, that is, region I. In this mode, the discharge forma- 
tion resembles a typical plasma column. As the current is increased 
further, the plasma cannot be contained within the inner well because 
of rising centrifugal forces and consequently jumps into region II. 
In this second mode, a rapidly rotating annular discharge is formed 
which is centered in the pressure well. Since this region contains 
a much larger volume, the plasma is capable of supporting a greater 
current density. 

In the tubular anode geometry, it was observed that if the 
current density was too high, or the magnet field too low, the centri- 


2 
) 
cé 


p 
fugal force | = | could exceed the magnetic radial confinement force 


(FB). At this point, the discharge could not be contained in 
region II and would subsequently blow out into region III, forming 

a glow-to-arc transition along some point on the electrode extremity. 
This process could be controlled by increasing the magnetic field 
strength thereby creating a deeper pressure well. Unfortunately this 
procedure caused the discharge to narrow. An alternative method for 
confining the discharge is to redesign the magnetic geometry to create 


a "saucer" shaped pressure profile which would maximize the discharge 


volume. However, when the anode was redesigned into the fluid ballasted 


geometry, this problem was essentially eliminated since the ballasted 
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anode was very effective in distributing the plasma over the discharge 
volume. 

A comparison of Fig. 6.4 with the photographs shown in the next 
few sections reveals that the theoretical predictions for the radial 
position of the discharge is in good agreement with the experimental 


results. 


6.2 ANALYTICAL MODEL AND RESULTS OF THE COAXIAL GEOMETRY 


The behavior of a glow discharge in a coaxial geometry can also 
be determined by a MGD analysis. In the structure shown in Fig. 3.4, 
the axial magnetic field B. in combination with the radial current 
J, » provides a Lorentz force which causes a rotation of the ionized 
particles in the negative 6 direction. The motion of the charged and 
neutral particles is coupled through collisions and results in a bulk 
- gas rotation in the same direction. The radial magnetic field B., in 
combination with the Hall current J. produces an axial Lorentz force 
which serves to confine the discharge to the central region of the 
electrode. 

The Hartmann number given by the expression 

H = BL rab 
n 

gives the ratio of magnetic forces to viscous forces. For the 
experimental parameters. appropriate to the magnetically stabilized 
gas discharge, Hs 1. This implies that lv, << E/B, and con- 


sequently, the Lorentz forces in the azimuthal and radial directions 
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can be approximated as 


(6.9) 


Pee -o7C0S($) °F B ; 


The force ig is negative for z > 0 and positive for z < O and there- 
fore serves to confine the plasma about z = 0 as indicated by Fig. 6.5. 
To simplify the analysis of this problem, Fig. 6.5 is employed 
to model the current emission from the rotating disk shaped plasma. 
Here it is assumed that 7, << J, so that uniform cathode emission 
occurs over the region -D/2 < z < D/2, and zero elsewhere. The current 


density may then be approximated in the following manner. 


Ly 
Deep? LOH Dieta ay SD 2 
AE = 
5 0 , otherwise (6.10) 
gon GB 
2 


where I, is the total discharge current, and D is the width of the 
discharge, usually 2.0 cm. 

The solution for this coaxial electrode geometry was initially 
carried out for a single fluid model as described in [6.2]. In that 
analysis, a steady-state solution, neglecting secondary flows, was 
obtained by considering the 6 component of the momentum conservation 
equation, that is, 


o> _ 
1s u (Vv oe oa Oe, (6.11) 


By using the approximation for the current density (Eqn. 6.10) and 


assuming that the azimuthal velocity was independent of 2, Eqn. 6.11 
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Fig. 6.5 Current Density Distribution Model for the Coaxial Geometry. 
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~ could be written as 


(6212) 


The one dimensional analytic solution to the partial differential 


equation in Eqn. 6.12 was then expressed as 


| C, 
v(x) = rilertetr in(e)] (6.13) 
where 
I,B 
a5 7 Anu D 
and 
Bo Ha 


1 2 
eid be ni: ain 
1 BS Re l Ry 2 


AiR, Ry 
one (R2-R2) in | A 

The parameters Be and Hare the inner and outer radial boundaries of 
the discharge. | 

A two-dimensional steady-state solution, which neglected 
secondary flows has also been carried out by numerical methods in 
[6.3]. Figure 6.6 illustrates the computed rotational velocity profile 
for a 4:16:20 torr gas mixture containing C0,./N./He with a total 
discharge curent of 7 A and a viscosity of 1.43 x 10% kg/(mes). A 
comparison of the peak rotational velocities determined from the one- 
and two-dimensional models is displayed in Fig. 6.7. It is clearly 
seen that both velocity profiles are similar, although the velocities 


obtained for the two-dimensional model is 40% less. This result can 
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Fig. 6.6 Rotational Velocity Distribution for the Coaxial Geometry 
Neglecting Secondary Flow. Discharge parameters include: 
J = 155 mA/cm2, »p = 1.43 x 107° kg/(mes). 
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Fig. 6.7 Comparison of the Peak Rotational Velocity Profile Between 
a 1D and 2D Model, Neglecting Secondary Flow. 
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be attributed to the presence of a larger active discharge volume with 
Bee ecoonieaiy larger viscous drag associated with the regions outside 
the uniform plasma emission. 

These initial models reveal that extremely large amplitude, 
highly sheared flow velocities can be expected to exist in the co- 
axial electrode structure. In the near future, the computer code 
MAGIC will be modified to determine the secondary flows generated by 
this geometry. Again, it is expected that the actual rotational 
velocities will be much lower than suggested by the initial models. 

In both electrode geometries, secondary flows are created as 
a result of pressure gradient forces which convect momentum into 
regions of high viscosity and low azimuthal drseiey The radial and 
axial secondary flows are eventually restricted by boundary walls 
where momentumis continually lost. The flows are then forced to 
return to the Hiscranae thereby forming a small convective cell. In 
the process of returning, the secondary flows intercept the region of 
high speed rotational flow causing the azimuthal velocity to be 
loaded down. Thus, the Lorentz driving force must perform additional 
work in order to maintain the azimuthal momentum. It is this dynamic 
process which is responsible for lowering the overall rotational 
velocity. From these results it is obvious that a modified con- 
vective cell is needed to create a less restrictive flow pattern which 


would substantially increase the rotational flow velocity. 
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6.3 OBSERVED PHYSICAL PHENOMENA 
6.3.1 Preliminary Observations 


During the initial phase of this project the effects of a 
magnetic field on the behavior of a self-sustained and a pulser 
Sustained glow discharge was studied. It was observed that in the 
absence of a magnetic field, the gas discharge would undergo a glow- 
to-arc transition at an electrode edge, when the applied voltage 
approached the gas breakdown voltage. To prevent this sudden tran- 
sition, a uniform preionization was established by means of a high 
voltage thyratron pulser. Now, as the gas breakdown voltage was 
approached, a pulser sustained, normal glow discharge column began to 
form between the electrodes. With a further increase in discharge 
current, the gas discharge entered the abnormal glow regime and 
eventually deteriorated into a glow-to-arce transiton. This problem 
in conventional laser systems is normally solved by employing. high 
speed blowers to spread the constricting plasma column and to provide 
convective cooling. However, it was proposed that magnetic fields 
could accomplish the same task but with improved performance and 
increased power loading. 

The general behavior of a magnetized glow discharge was 
determined for three different electrode geometries, through the use 
of high speed photography. The resulting photographs revealed the 
kinematics of motion present when a glow discharge is subject to a 
Lorentz force. In addition, an estimate of the rotational discharge 


velocity was obtained for botha transverse and coaxial electrode design. 
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A. Tubular Anode Geometry 


The first electrode geometry examined was a transverse design 
with a flow through tubular anode (Fig. 3.3). It was observed that the 
behavior of the glow discharge column was radically changed with the 
application of a magnetic field. Under the influence of a small 
magnetic field, the plasma column began to slowly move over the tubular 
anode in a circular motion. The direction of this motion was in accord- 
ance with the resulting Lorentz force. As the discharge current or 
magnetic field current was increased, the plasma column began to spread 
as it rotated with increasing speed. At this stage, the glow discharge 
resembled a uniformly distributed annular ring and gave little indica- 
tion that it was actually rotating. 

The cathode glow layer was also observed to be annular in shape. 
This particular characteristic shape is the result of radial Lorentz 
forces (F,B,) interacting with the epahed Fluid. This process can 
create a radial pressure gradient for the charged fluid as Meviously 
explained in Section 6.12. A detailed explanation of the physical 
mechanism for this discharge shape will be discussed in Section 7.2. 

The photograph in Plate 1.a, taken at 1/60 of a second, reveals 
some of the internal structure in the rotating discharge. However, 
it can be clearly seen that the discharge 1S still  oceupying “its full 
annular region. At this point, it was still unknown whether the 
discharge was simply a single plasma column rotating at a very high 
speed or whether the discharge was actually spread out over the entire 
electrode surface. Consequently, higher speed photographs were needed 


in order to resolve this dilemma. The photograph in Plate 1.b, taken 
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1. High Speed Photographs of a Tubular Anode 
Geometry. 


(a) 


Plate 1. Gas discharge parameters include: C0.,/N,/He laser gas mixture 


Qt orc oLorn. Tne = 640 V, Ino 


Photographs are of (a) 16.7 ms and (b) 1.0 ms duration. 
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at 1/1000 of a second, still shows definite rotational structure. How- 
ever, the glow discharge is now incomplete and appears to be just a 
smeared plasma column. This premise was corroborated by other photo- 
graphs showing an ionization trail as the plasma column rotated. 

It is well known that the plasma column is a region of intense 
current density and subsequently high conductivity. It is suspected 
that the uniform electric field profile that exists between the 
electrodes becomes distorted in the vicinity of the plasma column. 

As the gas discharge rotates, the distortion in the electric field 

_ follows the motion of the column. As shown in Fig. 6.8, the ions 
emitted from the anode surface are now subject to a 8-component in 
the electric field and are forced to follow the distorted field lines. 
The collective motion of numerous charged particles produces the 
corresponding ionization trail. 

While examining the behavior of arc discharges, it was observed 
that the application of a strong magnetic field was sufficient to 
extinguish an arc and form an annular shaped glow discharge. During 
this process, the arc looked as if it was slowly being pulled into 
the center of the electrode by the magnetic field. It is reasonable 
to speculate that the magnetic field slowly draws charged particles 
from the arc into the radial pressure well. This procedure would 
reduce the overall current density since charged particles would now 
be free to rotate in a large discharge volume under the governing 
influence of the Lorentz forces. The process of extinguishing the 
arc was most noticeable at the cathode surface where the magnetic 


field is strongest. It is therefore suspected that if the magnetic 
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field can stabilize the cathode glow layer, the remaining arc discharge 


must also follow suit. 


B. Coaxial Electrode Geometry 


The second electrode geometry examined, consisted of a coaxial 
design operating in a magnetic field produced by a solenoid. In this 
arrangement, a radial electric field combines with an axial magnetic 
field to produce a Lorentz force in the 6 direction as shown in 
Fig. 3.4. As in the previous geometry, a pulser was employed to pro- 
vide a uniform preionization over the discharge volume. The high 
voltage pulser created numerous streamers which radiated from the 
anode and formed small 1 cm diameter spots on the cathode surface. 

As the DC voltage increased, the spots became brighter and larger. 
Eventually these spots coalesced into a uniform cathode glow layer 
approximately 5 cm wide. 

As the gas breakdown voltage was exceeded, a discharge column 
was formed. With the application of a magnetic field, the radial 
discharge column began to rotate about the anode and rapidly accel- 
erated into a thin rotating disk. As the discharge current was 
increased further, the width of the rotating plasma body expanded 
longitudinally into the inter-electrode volume. The Hall current 
elie coupled with the radial component of the magnetic field, B 
produces a Lorentz force which confines the discharge axially to the 
central portion of the eiectrode Structure. This magnetic contain- 


ment feature is especially useful, since it becomes unnecessary to 
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profile the strong fringing electric fields created by the coaxial 
electrodes. 
The photographs in Plate 2 were each taken at 1/1000 of a 


second with a gas composition of He and N. (10:10 torr) and 20 torr of 


2 
C0, » respectively. Both photographs indicate that the disk is 
rotating at very high speeds since no internal structure could be 
detected. Thus, for the 7 cm diameter coaxial electrode, the estim- 
ated minimum rotational velocity is 200 m/s. This very high velocity 
is due to a large driving force acting upon a small active discharge 
volume. In addition, both coaxial electrodes are smooth and con- 
tinuous and consequently do not subject the rotating discharge to 
excess turbulence. 

A wide variety of gases and mixtures were tested in the coaxial 
structure including Ar, He, No » and C0. - All gas combinations 
produced excellent results for pressures exceeding 100: torr. Stable 
high power discharges were also obtained in electronegative gases such 
as 0, and SFE ; } 
particularly interesting since observations revealed an absence of 


The electric discharge containing 5 torr of SF was 


a Faraday dark space. This effect was due to the extremely high 
electron affinity that exists for this electronegative gas. Conse- 
quently, this rapid recombination process to form negative ions, 
prevents the formation of any large concentrations of electrons. 
These tests have demonstrated that the magnetic stabilization 
technique is largely insensitive to gas composition and pressure. 
Furthermore, in order to uniformly distribute the rotating disk 
discharge over a larger active volume, it may be necessary to 


construct a fluid ballasted electrode. 
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2. High Speed Rotating Coaxial Discharge. 


Plate 2. (a) Gas discharge consists of a N,/He mixture at 10:10 torr; 
Ve 480 V, lip OG4IN aes 400 G. (b) Gas discharge 
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Both photographs are of 1.0 ms duration. 
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C. Multi-Pin Ballasted Anode Geometry 


The previous experiments with the transverse tubular anode 
geometry indicated that the gas discharge operated as a localized 
plasma column. The resulting small active discharge volume made very 
inefficient use of the potential volume between the electrodes. In 
an attempt to solve this problem, a multi-pin ballasted anode 
(Fig. 3.5) was constructed to replace the tubular device. The bal- 
lasted anode was designed to create a uniformly distributed plasma 
through an array of sub-electrode pins. Since the other components 
of the transverse structure remained unchanged, a thorough comparison 
of the two electrodes could be performed. 

To initiate the glow discharge, a high voltage pulser was first 
applied to the anode. Due to the performance of the individually 
bal jaeved sub-electrode pins, a homogeneous preionized discharge 
volume was established. Upon initiation of the gas breakdown voltage, 
a widely dispersed: glow discharge began to form in the center of the 
electrodes. With the application of a magnetic field, the glow dis- 
charge rapidly accelerated into a uniformly distributed plasma, 
covering the entire electrode surface. The resulting gas discharge 
was only slightly annular in shape, being dispersed by the individ- 
ually ballasted sub-electrode pins. Thus, it can be surmised that the 
radial Lorentz forces are still playing a factor in determining the 
gas discharge characteristics. 

The photographs in Plate 3, taken at 1/250 and 1/1000 of a 
second, revealed an evenly distributed plasma with little internal 


structure. As a result, the rotational discharge velocity could not 
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3. High Speed Photographs of a Discharge in a Multi- 
Pin Ballasted Anode Geometry. 


Plate 3. Gas discharge parameters include: CO,/N,/He laser gas mixture 
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Photographs are of (a) 4.0 ms and, (b) 1.0 ms duration. 
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be accurately inferred from these photographs. Subsequent determina- 
tion of the parameter will require more sophisticated diagnostic 
techniques. 

Since the entire discharge volume was now more effectively 
utilized, the local discharge current density and luminosity was 
Significantly reduced compared with the tubular anode device. This 
uniform distribution of current density enabled the electrode to 
handle a large increase in overall discharge current. Typically, 
the maximum value of current density achieved was 50 mA/cm2. Through 
a variety of experiments, it was also determined that the electrode 
system was capable of stable operation at pressures exceeding 60 torr. 
These results were repeated numerous times with an assortment of gas 
mixtures. Once again the magnetic stabilization approach has proved 
highly successful and thus far, has only been limited by power supply 


restrictions. 


6.3.2 High Speed Photographic Interpretation 


The electric discharges produced by the different electrode 
geometries were all initially analyzed, under various conditions, 
with a 35 mm single lens reflex camera. The resulting high speed 
color photographs revealed that a uniform glow discharge had been 
established in the coaxial and ballasted anode geometries. However, 
it was observed that the supposedly uniform glow discharge found in 
the tubular anode geometry was actually a plasma column rotating at 
high speeds between the electrodes. In order to study this phenomenon 


further, a high speed framing camera was employed. The gas discharges 
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were examined with a 16 mm DBM-54 Milliken motion picture camera 
capable of frame rates in excess of 400 frames per second. The 
Shutter used in this special camera consisted of a rapidly rotating 
plate with a small adjustable slot. By varying the slot width, a 
shutter speed as low as 1/4000 of a second could be achieved. 

Previous experiments have indicated that the cathode glow 
layer was also rotating under the influence of the Lorentz force. 
High speed photographs revealed that the cathode glow layer only 
existed as a half-moon shaped region on the cathode surface. As 
stated by Nasser [3.2], a proper glow discharge cannot function with- 
out a cathode glow layer. Now, since the plasma column is strongly 
influenced by the magnetic field, the positive column must move in 
unison with the cathode glow layer. This effect was observed by 
viewing the film taken by the high speed framing camera. The sequence 
of photographs shown in Plate 4 show the cathode glow layer rotating 
in conjunction with the positive column. Furthermore, each frame in 
this series confirms that the rotating discharge is actually a single 
column rotating between the electrodes. 

From the film sequence, the positive column was estimated to be 
2-3 cm in diameter. Also observed was the emergence of a radial after- 
glow from the top of the positive column near the cathode. Due to the 
loss of these particles, this region became diffuse and less luminous 
than the bulk of the positive column. By examining the film, frame by 
frame, the rotating discharge column was observed to jump along the 
tubes of the anode as it made its way around the electrode. It is 


believed that this process reduced the velocity of the rotating plasma 
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column since a new ionization channel had to be established with the 
cathode glow layer. 

The cathode glow layer was estimated to cover an area of 50 cm2 
on the cathode surface. Since this discharge region is more extensive 
than the positive column, it was hoped that further increases in the 
applied driving force would be sufficient to enlarge the cathode glow 
layer and thus extend the positive column. Unfortunately, the magnetic 
field intensity was limited by saturation effects. In addition, the 
power supply was unable to provide the required amount of current to 
test this speculation. 

By examining the film sequence in Plate 4 it becomes possible 
to estimate the velocity of the rotating column. This can be accom- 
plished by simply counting the number of frames that the discharge 
takes to complete one revolution around the anode. Since the framing 
speed of the motion picture camera was 400 frames/s, each photograph 
was taken at 2.5 ms intervais. Viewing this 16 mm film on a Film 
editing machine, indicates that for the conditions specified, there 
are 5 frames per revolution. Subsequently, the rotational velocity 


of the plasma column is given by 


» = md: framing speed 
# of frames/rev. 


where d is the diameter of the dishcarge. This formula indicates a 
rotational velocity of v = 27.6 m/s which confirms the results of 


the photodiode experiment discussed in Section 6.4.1. 
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6.3.3 Afterglow Effect 


As the plasma rapidly accelerates to a steady-state value, an 
increasing radial pressure profile begins to develop across the dis- 
charge. The resulting centrifugal force will in turn drive secondary 
flows. Since the applied driving forces are stronger near the cathode, 
larger rotational velocities are grea ad in this region. Consequently, 
the outward flow of gases will also be predominantly emanating from 
the cathode. In the transverse geometry, this flow will be radially 
outward along the cathode surface. To complete the convective cell, 
the inward flow is expected to return upwards through the tubular or 
multi-pin flow through anode. In the coaxial geometry, the outward 

flow will be driven axially along the outer electrode and return 
axially along the inner electrode. 

This outward flow of gas was visually confirmed in all electrode 
geometries, by the light emitted from metastable particles. The light 
in the form of an afterglow could be seen spewing outwards from the 
Faraday dark space and cathode regions. Further examination by the 
high speed film sequence, indicated that the radial afterglow actually 
originated from the top-half of the positive column as it rotated over 
the tubular anode surface. 

The metastable particles are normally the result of a recom- 
bination process involving the capture of an electron by a positive 


ion according to the reaction 
+ 
X +teoXx +h. 


From this reaction, a metastable particle is produced and a photon is 


emitted. These metastable particles can only be created when the 
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captured electron cannot return to a lower energy level without 
violating certain selection rules. These long-lived metastables can 
only be restored back to their ground state through an energy absorb- 
ing process. This will allow the electron to move to a higher energy 
level from which it can then return to its ground state. The energy 
absorbing nnuceee is usually the result of inelastic collisions. 

The afterglows were visually observed as a result of the photons 
emitted during the metastable production process. A bright orange 
afterglow was observed for large concentrations of No Coez0etorr) 
while a light green afterglow was produced in a C0, gas discharge. 
Since the electric discharge dissociates C0., into CO and 0, eet 1s 
believed that the C0. afterglow is associated with the presence of 0, : 

The afterglow emanating from the rotating ionized gas is known 
to result from collisional processes occurring within the discharge. 
Since the metastables are uncharged, they are not bound to the charged 
particle pressure well. Consequently, these uncharged particles are 
convected radially outwards, in the transverse geometry, by the cen- 
trifugal pumping aside of the rotating gas discharge. In the coaxial 
geometry, this flow was generated axially. The afterglow velocity was 
experimentally determined by the Pitot tube measurements discussed in 
Section 6.4.2. These results indicate that the afterglow velocity was 
- proportional to the magnetic field and the current density, that.is, 
to the applied driving force. The significant outward flow generated 
by these forces can be particularly important in the design of a laser 
system since the effects can be exploited to recirculate and cool the 


gas. 
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6.3.4 Standing Wave Striations 


One of the more prominent characteristics exhibited by the 
rotating gas discharge concerns the occurrence of fluctuations in 
light intensity. In the transverse geometry, these fluctuations 
appeared as horizontal bands in the plasma column. These fluctuations 
were viSually observed upon careful inspection of the gas discharge. 
Through the use of high speed photography, these bands were clearly 
evident since the amount of light entering the camera had been reduced 
(Plate ig The luminous bands were also recorded by a high speed 
framing camera. An examination of a sequence of frames revealed that 
the light intensity fluctuations were stationary even as the plasma 
column rotated. It therefore became evident that these stationary 
luminous bands were the result of standing wave striations in the 

positive column. 

A detailed study of the coaxial electrode geometry again 
revealed the occurrence of fluctuations in light intensity in the 

form of concentric rings. These luminous rings, located in the 
rotating positive column were visually observed and photographed. The 
pictures displayed in Plate 5, were taken of gas discharges containing 
10: 10° tomnaog He/N., and 20 torr of CO, respectively. The number of 
rings varied according to the discharge current density and applied 
magnetic field, that is, to the ina) driving force. As the Lorentz 
force increased the number of rings grew until they merged and became 
indistinguishable. This effect can be illustrated by Plate 2 where 


only one ring remains. This finding is consistent with those of 
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5. Standing Wave Striations Observed in a Coaxial 
Geometry. 


Plate 5. (a) Gas discharge consists of a N,/He mixture at 10:10 torr; 


no = 425 V, Ina 


duration. (b) Gas discharge consists of 20 torr C0.; 


“DE B.», ~ 300 G. Photograph is of 8 ms 


Vine = 660 V, Inn a7. 0 A 5 Be. = 300 G. Photograph is of 


16.7 ms duration. 


- 
if Ve 
1 
<, 
° 
r 
AN 
! 
1 
» as 
; " 
j 
Lo 7 
abe 
) 
. ee 
’ 
y 
vi 
' 
' ee 
; 
A ‘ 
' 
‘@ 
.@ 
ar, 
ey | 
Mt 
a 


i - 
Rie Yne : ; de 


i? ‘cdiiial C Loh pare 7 ) 
ma 8 POE Hootie Pat. | 
oF ee ea gist 


a he veenpedady ig 


| rn 2 hie ad if vie 


v ‘ 


160 


| Crandall and Cooper [6.4] who observed a single standing striation 
usually located near the cathode. 

Since the gas discharges under examination are highly col- 
lisional and only partially ionized (n/n =2x 10°’), it is highly 
probable that the standing wave striations are due to acoustic waves 
propagating in the positive column. Within an ionized gas, an acoustic 
wave can be generated as a result of an external disturbance. If this 
disturbance is due to a local change in the electric field, an increase 
in the dissociation of particles may result. The ensuing local 
pressure increase can then initiate an acoustic wave propagating in 


the neutral gas at the wavespeed 


The passage of the acoustic disturbance is recognized by a local 
increase in light intensity due to a rarefraction wave, and a local. 
decrease in luminosity for a compression wave [6.4]. The light 
intensity decrease with compression is attributed to a decrease in 
the electron temperature resulting from an increased collision fre- 
quency. The less energetic electrons are now unable to fully excite 
the neutral particles to produce spontaneous emission of light. 

The analysis of Pekarek and Krejci [6.5] has shown that the 
formation of a striation is caused by a disturbance in the local 
electric field. The resulting perturbation in the electric field 
causes a rise in the rate of production of ion pairs, which in turn, 
produces a positive space-charge. The electrons are easily separated 
from the ions, due to their large diffusion coefficient, and eventually 


create a separation of charge carriers. As the positive ion density 
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accumulates, an electric field produced by the space-charge begins to 
oppose the original perturbed field. Consequently, the perturbed 
electric field is reduced and the positive ion density begins to 
decrease in this discharge region. Further along the discharge axis, 
the perturbed field now combines with the space-charge field and a 
new striation is formed. The analytic solution of the resulting 
perturbed electric field and positive ion density, after an initial 
aperiodic disturbance, can be described by an oscillating Bessel 
function. Thus, it can be shown that moving striations will result 
from the superposition of many repeated disturbances in the electric 
field 16.5); 

It is suspected that the initial disturbance in the electric 
field is created by the periodic application of the high voltage 
pulser. Unfortunately this recently known cause was not verified at 
the time of the experiment. It is also believed that the occurrence 
of standing waves rather than moving striations can be attributed to 
the use of parallel electrodes in both the tubular anode and he 
coaxial systems. This was especially true in the coaxial geometry 
with its smooth concentric electrodes. These solid electrode walls 
were perfectly reflecting and thus ideal for the production of 
acoustic standing waves. While operating the multi-pin ballasted 
anode system, no standing wave striations were observed, since the 


pin anode is a poor reflecting surface for acoustic waves. 
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6.3.5 Dark Space Width 


Another interesting feature exhibited by the magnetized glow 
discharge pertains to the dark space width and its dependence upon 
the applied magnetic field. It was observed that for the three 
electrode geometries, the physical width of the cathode dark space 
could be easily varied by an order of magnitude ce adjustment of 
the magnetic field. For very large field values, the dark space 
could be reduced to within 2 mm of the cathode surface. 

The variation in dark space width is attributed to the effective 
reduction of the ratio #/n due to the presence of a magnetic field. As 
reported earlier in Section 4.5, F/n was reduced by over 65% in the 
cathode fall region compared to its non-magnetized state. This sub- 
Stantial change in electric field will significantly reduce the 
electron temperature and correspondingly its energy as it accelerates 
through the cathode fall. The maximum energy attained by these 
accelerating electrons occurs at the end of the cathode dark space. 
Here, they participate in inelastic collisions with the negative glow. 
Normally, Since the electron energies are too high, their inelastic 
collision cross sections are small. Consequently, the excitation and 
jonization rates only increases after the electrons have been slowed 
down. However, with the reduction in cathode fall voltage, due to 
E/n, the less energetic electrons are able to produce inelastic col- 
lisions sooner. This effect has been substantiated by Brewer and 
Westhaver [6.6] who have correlated the width of the negative glow 


with the electrons which have gained energy in the Crookes dark space. 
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‘Their results show that as the electron energy is reduced, there is 
a corresponding decrease in the width of the negative glow. 

From the analysis presented, a reduction in #/n leads to an 
enhanced electron collision frequency. Therefore, the addition of a 
magnetic field is equivalent to an increase in pressure. Under these 
conditions, the negative glow and the two dark spaces that surround 
it, will shrink towards the cathode, and the positive column will 
extend to fill the void. Subsequently, the positive column is estab- 
lished closer to the cathode and more of the active discharge volume 
can participate in laser excitation processes. A more theoretical 
approach to describe this change in gas discharge behavior is des- 


cribed in Chapter 7. 


6.4 GAS DISCHARGE DIAGNOSTIC TECHNIQUES 
6.4.1 Photodiode Measurements 


The results of the high speed photographs revealed that the 
rotating gas discharge was not continuous, but merely a single plasma 
column moving rapidly over the tubular anode surface. Though high 
speed motion pictures gave a good indication of the azimuthal velocity, 
a more convenient and accurate approach was needed. 

It was proposed that a sensitive photodiode could be used to 
detect the light emission from the plasma column. As the column 
rotated, the resulting periodic signal would be used to determine the 
discharge velocity. The experimental photodiode apparatus is illust- 


rated in Fig. 6.9. In this schematic, two shielded photodiodes were 
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employed to detect the passage of the plasma column as it rotated 
around the tubular anode surface. The initial deployment of two 
photodiodes were used to ensure that the rotational velocity was 
indeed constant. Each photodiode was mounted against a 3 mm diameter, 
15 cm long brass collimating tube. The long narrow tubes were nec- 
essary in order to collect the incident light and prevent stray light 
from reaching the photodiodes. The tubes were mounted in parallel 

6 mm apart so as not to overlap the incoming spectral signal. 

During operation of the discharge, a pulsed beam of light was 
received by the photodiode with each revolution of the plasma column. 
This signal was then amplified by an emitter follower circuit. The 
resulting voltage peaks were viewed on an oscilloscope and recorded. 

A typical voltage trace can be seen in Fig. 6.10. Note the occurrence 
of a primary and secondary voltage pulse. As the plasma column passes 
point A, it is detected by the first photodiode. A Short time later, 
At, the second photodiode observes the passage of the column past 

point B. As the plasma column continues around the electrode it passes 
point C where it is seen by detector 2 and detector 1, as it passes 
point D, At seconds later. Note that the detected amplitude of the 
discharge column light intensity, as it passes the far side, is much 
smaller and broader. This effect is due to the increase in the depth 
of field seen by the photodiode. This experiment was repeated, but 
with the magnetic field polarity reversed. This action changed the 
orientation of the Lorentz force and accordingly reversed the direction 


of discharge rotation. 
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Fig. 6.10 Oscilloscope Recording of Photodiode Voltage Measurements. 
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The rotational velocity of the plasma column was initially 
ealeulated by measuring the time delay, At, between detected pulses 
using both photodiodes. These initial measurements provided con- 
clusive evidence that the rotational speed was constant. Since the 
detected light pulses were also periodic, only one photodiode was now 
‘ required. The rotational discharge speed could now be calculated by 


the formula 


ae 


the period of revolution determined by the oscilloscope time base. 


where nd is the peripheral path taken by the discharge, and aa 


As shown in Fig. 6.11, an increase in magnetic field and/or 
applied current is accompanied by a rise in rotational speed. The 
dotted lines indicate the expected initial velocity relationship. 

After a proper plasma column has been established, it is clearly seen 
that this relationship becomes linear. These eepoerederetiits cor- 
respond to a linear increase in the applied I xB driving force. 
However, as indicated in Fig. 6.11, the magnetic flux density becomes 
slightly nonlinear as the magnetic field current is increased. This 
effect is due to saturation of the magnet core. Since the magnetic 
flux density is of nonuniform strength throughout the discharge volume, 
it was decided that the applied magnetic field current is best suited 
to represent magnetic field effects. This parameter is easily measured 
and is strictly linear with respect to the magnetic field; neglecting 
iron core saturation effects. However, as a rough guide, each ampere 
of magnetic field current corresponds to an average magnetic flux 


density of approximately 200 Gauss. 
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Fig. 6.11 Dependence of Plasma Column Rotational Velocity on Applied 
Discharge Current. Gas discharge contains a 2:8:10 torr 
mixture of C0,/N,/He. 
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The width of the plasma column was estimated by measuring the 
half maximum width Tee of the light intensity pulse received by the 
photodiode. Since the rotational period is known, the discharge column 


width d,, >» 1S calculated through the ratio 


As Fig. 6.12 indicates, for a given magnetic flux density, the plasma 
column expands in response to the applied discharge current. This 
expansion is the outgrowth of two possible competing effects. First, 
it is probable that the increased current density stimulates further 
ionization which eatvanutee to the expansion of the plasma column. 
Secondly, an increase in the Lorentz driving force would spread the 

discharge column due to acceleration of the rotating charged particles 
| in the Yo direction. 

It was also noted that the discharge column width decreases with 

applied magnetic field. It is apparent that the constricting plasma 
column is caused by the charged particle pressure well which narrows 
in response to an increase in magnetic field. The largest diameter 
of the plasma column occurred in the absence of a magnetic field, 
however, this was an unstable situation and a glow-to-arc transition 
could easily develop. At higher pressures, the plasma column simply 
constricted to the center of the electrode. At this point, the radial 
magnetic field component, responsible for creating the rotational 
driving force, was quite small and consequently no rotation was 


observed. 
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Fig. 6.12 Variation in Plasma Column Width. 
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One further observation noted in this experiment, concerns the 
intensity variation of the plasma column. It was ascertained that the 
discharge luminosity fell off rapidly as the photodiode axially scanned 
from the anode to the cathode. This investigation confirmed results 
_ from the high speed photographs which revealed a substantial expansion 
of the positive column near the cathode. This observed effect is 
attributed to the presence of a stronger Lorentz force, near the cathode, 
which assists in spreading the discharge column. 

The photodiode apparatus was not implemented to determine the 
discharge velocity achieved by the coaxial elecends design. It was 
felt that since both the transverse and coaxial rotating plasma 
experiments were based upon the same principles, concentration on 
both projects would deviate from the central theme of this thesis. 
Nevertheless, the coaxial discharge geometry proved interesting and 
worthy of continued study. 

As a result of the undesirable qualities manifested by the 
tubular anode; a fluid ballasted multi-pin electrode was constructed 
to provide a uniformly distributed glow discharge. Since the new 
anode evenly spreads the glow discharge over the entire electrode 
surface, the rotating plasma column phenomena was no longer apparent. 
The photodiode apparatus was used to confirm this assumption. The 
resulting investigation only revealed a uniform background light 
intensity level present along any horizontal plane. The light 
intensity level did become more diffuse and eventually disappeared 
as the detector scanned the positive glow and the Faraday dark space, 


as in the previous case with the tubular anode. These observed facts 
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suggests that under the influence of a magnetic field, the new electrode 
design produced a uniformly rotating plasma with no discernible 


Structure. 


6.4.2 Pitot Tube Measurements 


The rotating gas discharge was suspected of also producing a 
radial gas flow that might originate due to the presence of powerful 
centrifugal forces. This suspicion was confirmed through the obser- 
vation of a radial afterglow being ejected out of the gas discharge. 
An experiment was then devised to measure this flow and determine its 
characteristics. As illustrated in Fig. 6.9, a Pitot tube was mounted 
to measure the static and stagnation pressures. A capacitance man- 
ometer was then employed to sense the difference in these pressures 
and produce an electronic output signal. From these results, a radial 
flow velocity can be determined. 

From any standard gas dynamic text, the ratio of static to 


Stagnation pressure is given by the expression 


_ = ; - este (6.14) 


where y = Grey, is the ratio of specific heat capacities, and /¥ is the 


Mach number. Rearranging this equation yields 
ie 


ae ray 1- a u (6.15) 


y-1 ie 


where M = vie. ites vyRP is the sonic speed and FR is the universal 
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gas constant. With these substitutions Eqn. 6.15 becomes 


y=) )2 
2yRT P, Y 
aia et 1- whe . (6.16) 
0 
phe ty, 18 Se 
It is also known that [> = = so 
yt 
P v 
ve 2c, \1 - > § (6! 17) 
0 


it then becomes convenient to simplify Eqn. 6.4 through the use of the 
binomial theorem, that is, 


Gena 


ee 


0 


For a laser gas mixture, the overall specific heat capacity per 


unit mass, at a constant pressure, is given by [6.7] 


vais 7 5 
k dela |S gn akhy 15h oon! fiz fue) 
Cy = Saerery = poe . (6.19) 
phi To00|4f Het ohn mae co,| 


In this expression, k is Boltzmann's constant, N, is Avogadro's number, 
f; is the fractional concentration of the zth gas species, and v. is 


the number of degrees of freedom for the zth species. 
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The expression for the radial gas flow velocity can now be 


written as 
1, 
i Bel een) le 
a= frog of1-YfreZ}}r (6.20) 
ee ; 
merey = oR ° and ec, 1s defined from Eqn. 6.19. 


The differential pressure, as determined by the capacitance 
manometer, was correlated with the appropriate gas mixture to yield 
the velocity curves displayed in Figs. 6.13-6.14. These results were 
obtained using a 1:7:12 torr gas mixture consisting of CO,/N,/He. The 
curves of Fig. 6.13 show the change in radial gas flow as the Pitot 
tube detected pressure changes just outside the active discharge 
volume. As expected, from visual observation of the afterglow, the 
radial gas velocity measured near the cathode region, was substantially 
higher than the measured values at the center or near the anode. These 
measurements indicate that the centrifugal forces are definitely higher 
near the cathode. An analysis of the magnetic field profile within 
the discharge volume, confirms this conjecture. It was found that the 
magnetic flux density was significantly larger in this region and 
thereby contributed to a larger Lorentz driving force. 

In the absence of a magnetic field, the Pitot tube apparatus 
detected no pressure differential. With the initiation of a magnetic 
field, the plasma column began to rotate. However, no radial gas flow 
was detected until the rotating gas discharge was fully established 
and an afterglow observed. This effect is speculated by the dashed 
lines in Fig. 6.14. In addition, it was noticed that a stronger 


magnetic field was needed to fully establish a rotating discharge for 
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Fig. 6.13 Effect of Applied Discharge Current on the Radial Gas 
Velocity Profile. Gas discharge contains a 1:7:12 torr 
mixture of C0,/N,/He with Be =# £00 G. 
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Fig. 6.14 Dependence of Radial Gas Velocity on Magnetic Field. Gas 
discharge contains 1:7:12 torr mixture of C0,/N,/He. 
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_ larger values of applied current. It is believed that this effect is 
due to a stronger Lorentz force required to overcome the high prob- 
ability of a glow-to-arc transition for discharges with large current 
densities. 

With further increase in magnetic field strength, the afterglow 
noticeably changed its shape by expanding radially. Reducing the 
magnetic flux density, immediately caused the afterglow to shrink in 
size and lose some of its luminosity. By maintaining a constant 
discharge current, an increase in magnetic field produced a linear 
rise in the detected radial gas flow. Figure 6.14 illustrates this 
dependence for several values of discharge current. 

For a constant magnetic field, Fig. 6.13 also reveals that the 
radial gas velocity is an increasing function with respect to the 
applied discharge current. This result was also expected since the 
driving force is proportional to the current density. Thus, for large 
values of ¢ x B, a strong centrifugal force is created. The effect 
of this force, on the rotating plasma column, would tend to throw Of 
neutral particles since they are strongly coupled to the motion of the 
charged particles by viscosity. However, the neutral particles are 
not bound by the applied electric and vuanatic fields as are the 
charged particles. Therefore, they are free to establish a convective 
cell in which momentum can be transported out of the discharge volume. 
The resulting centrifugal pumping action is thus able to quickly 
remove hot spots from the discharge and replace it with fresh, cool 


gas. 
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6.5 VI CHARACTERISTICS OF A MAGNETIZED GLOW DISCHARGE 
6.5.1 Glow-To-Arc Transition 


It has long been known that most high power CW laser gas dis- 
charges operate in the abnormal glow region [1.2]. This region is 
characterized by a rapidly increasing plasma impedance as shown by 
the lower curve in Fig. 6.15. To maintain the gas discharge in the 
abnormal glow region it is necessary to continually increase the 
_external ballast resistance. This procedure will shift the load line 
to a more stable point on the VI characteristic curve, at the expense 
of system efficiency. Failure to follow this guideline, will ulti- 
mately lead to a glow-to-arc transition, as speculated by the dashed 
curve in Fig. 6.15. 

Under the influence of a magnetic field, the same gas discharge 
now has the characteristics of a normal glow, as illustrated by the 
upper curve in Fig. 6.15. In order to understand this marked dif- 
ference in gas discharge behavior, it is instructive to follow the 
transition from a normal glow to an arc discharge. While operating 
in the normal glow regime, the discharge only partially covers the 
cathode surface. As the discharge current is increased, the current 
density and voltage remain constant. Consequently, the area of the 
discharge must increase proportionately. At the termination of the 
normal glow region, the cathode is completely covered by the glow, and 
its entire surface serves as an electron emitter. In order to further 
intensify the discharge current, it becomes necessary to increase the 


number of electrons created by ionization in the cathode fall. This 
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Fig. 6.15 Comparison of VI Characteristics for a Magnetized and Non- 
Magnetized Glow Discharge. Gas Discharge contains 1:7:12 torr 
mixture of CO,/N,/He. 
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process can be accomplished by increasing the voltage, and hence the 
electric field, in order to stimulate the rate of ionization by electron 
impact. The characteristic rise in voltage signals the onset of the 
abnormal glow. Since the cathode now receives a larger number of 
positive ions, each with a higher energy, its temperature rises. This 
process continues until the thermionic emission from the cathode 
Surface becomes comparable to the electron emission by positive-ion 
bombardment. Since the thermionic emission current is a strongly 
increasing function of temperature, an additional rise in discharge 
current leads to greatly increased electron emission. The excess 
electrons increase the discharge conductivity by initiating an inter- 
electrode ionization channel or streamer. Upon further influx of 
energy, the streamer builds in intensity to form an arc discharge, 


characterized by a high current and a low voltage. 


6.5.2 Behavior as a Normal Glow 


A comparison of the two VI curves in Fig. 6.15 reveals that the 
magnetized glow discharge has an extended normal glow region. It is 
suspected that as the cathode glow layer expands to cover the cathode 
surface, an abnormal glow will eventually develop. Thus far, this 
transition has been far beyond the range of the existing power supplies. 
It is therefore believed that the magnetic field is creating an effect 
“which prevents or delays the occurrence of an abnormal glow and sub- 
sequently a glow-to-arc transition. This section will explain the 
physical mechanisms involved in maintaining an otherwise unstable 


gas discharge as a normal glow. 
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As reported in Section 6.1.2, the annular shaped discharge is 
confined to a radial pressure well as a result of the radial Lorentz 
pressure Z, x B - In addition, this force interacts with the space- 
charge sheaths, creating an annular shaped glow layer which only 
partially covers the cathode. Consequently, it would appear that the 
Lorentz forces can inhibit the expansion of the electron space-charge 
Sheath and the positive column, thereby preventing the transition to 
an abnormal glow. 

It has been determined from the theoretical analysis in 
Section 6.1.2, that the discharge can only rotate in regions where 
the forces created by electric and magnetic fields are balanced with 
radial pressure forces. Furthermore, the magnetic field profile 
controls and continuously adjusts the emitting cathode area to supply 
the desired discharge current. Therefore, the cathode becomes a more 
efficient emitter of electrons by producing more current through the 
cathode glow layer. At present, there are two known effects which 
can account for this process; involving the use of ETecercns and/or 
positive ions. 

In the presence of a magnetic field, electrons and ions are 
required to gyrate about magnetic lines of force. It has been 
speculated that near the cathode, some electrons will acquire a suffi- 
ciently small Larmor radius such that they will be forced to return 
to the cathode, causing further secondary emission. This process is 
possible provided the electrons have not undergone any collisions. 
This can be expressed mathematically as 


en, Sie (6.21) 
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where hn is the mean free path, and Yr, is the Larmor radius defined by 


aur <e,> fos, - Equation 6.21 can be rearranged to yield an expression 


for the required magnetic field necessary for electron back-diffusion, 


thats ti's?, m4 
; pe 2 cs 10 yi 
——_—— 2 SS 
B > ax (Snm kL) 1, , (6222) 


This equation indicates that for typical gas discharge conditions, 
the electrons emitted from the cathode will be too energetic to be 
affected by the moderate strength of the magnetic field. Only electrons 
with energies less than 1 eV will have a chance of returning to the 
cathode. However, these electrons have insufficient energy to cause 
secondary emission and will thus be simply recaptured by the cathode. 
Ion kinetic behavior can also be strongly influenced by a 
magnetic field. In particular, the applied Lorentz force creates a 
Significant azimuthal] velocity component. However, due to the non- 
linear nature of the magnetic field, the rotational driving forces 
- are strongly spatially dependent. Thus, the highly effective rota- 
tional shear will cause the ions to bombard the cathode surface at 
angles far from normal incidence. It has been wel] documented that 
the secondary electron emission yield can be effectively increased by 
over 100% for ionic angles of incidence greater than 45° [6.8]. The 
reason for this increased emission yield is due to secondary electrons 
being produced closer to the cathode surface for angles of oblique 
incidence. Thus, these newly liberated electrons have a higher 
probability of reaching the surface with correspondingly higher energy. 
There is some experimental evidence to support this latter 


hypothesis. In particular, after running the gas discharge for an 
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extended period, the cathode surface, upon examination, revealed a 
Surface texture resembling a myriad of fine, closely spaced, hair- 
ine Scratches. These scratches, up to 1 cm in length, were all 
parallel and uniformly distributed in the directon of discharge 
rotation. This effect was particularly noticeable for an aluminum 
cathode. The resulting scratches are believed to be due to sput- 
tering action of the cathode's thin oxide coating by bombarding ions 


near grazing incidence [6.9]. 


6.5.3 Voltage-Current Curves 


As shown in the previous section, the VI characteristics of a 
magnetically stabilized discharge, indicate that the plasma operates 
at a near constant E/n within the normal glow regime. Since the gas 
discharge was extremely stable, even for large current densities, only 
a small amount of external ballast was required for operation, 
typically about 15%. Further analysis was necessary to determine 
whether these beneficial characteristics could be affected by other 
parameters. This section will examine the operating conditions of 
the gas discharge and its dependence upon geometry, gas mixture, and 
pressure. 

Through numerous experiments, the MGD stabilization approach 
was found to work effectively with each geometry considered. In each 
case, a Stabilized gas discharge was observed to operate as a normal 
glow. However, it was quickly determined that the magnetic field was 
not of sufficient strength or configuration to establish a uniform 


plasma throughout the discharge volume. During operation of the 
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tubular anode geometry, the plasma preferred to remain as a discharge 
column which rotated between the electrodes. In addition, the high 
speed rotating plasma, found in the coaxial geometry, preferred to 
exist as a thin disk. Thus, it can be surmised that the primary 
configuration of an electrical discharge is one that will occupy the 
least amount of volume in order to remain self-sustained. 

Continued study of the gas discharge systems, revealed that an 
externally applied force has no ability to create additional current 
sources on the electrode surface. It can only distribute the existing 
current source in the direction of the applied force. As a conse- 
quence, the establishment of a uniform glow discharge is only partially 
promoted. 

This problem was overcome through construction of a fluid 
ballasted, multi-pin anode. In this geometry, each pin electrode 
creates a small uniform glow discharge between it and the cathode. 
The resulting array of pins establishes a uniform glow discharge 
which fills the entire active volume. Now, the application of a 
magnetic field can impart convective motion to each individual dis- 
charge column, thereby forming a uniformly distributed rotating 
glow discharge. 

Further testing revealed that the magnetic stabilization 
technique was relatively insensitive to the type and mixture of 
“gases utilized. |All were found to behave as a normal glow in 
the presence of a magnetic field. This included the electronegative 
gases 05 and SF. » which could be stabilized when subjected to large 
magnetic fields. Accordingly, the magnetic stabilization technique 


may find numerous applications for a variety of laser gas mixtures. 
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These tests were repeated to test the effectiveness of this 
technique on varying gas pressures and magnetic fields. Figure 6.16 
illustrates the VI characteristics of the magnetically stabilized 
coaxial discharge for a 2:8:10 torr mixture of CO,/No/He. beats 
clearly seen that the magnetic field can appreciably change the 
operating discharge voltage. This variation in voltage is approx- 
imately linear with the applied magnetic field. However, as stated. 
previously, saturation of the electromagnet core tends to reduce the 
field strength for larger values of field current. In Fig. 6.17, 
the partial pressures of each gas have now been doubled. The result- 
ing curves are similar in shape to the previous case, but show higher 
operating voltages. This expected result is due to the larger break- 
down voltages required for gases at higher pressures. By comparing 
Figs. 6.16 and 6.17, it can be stated that the influence of a magnetic 
field is equivalent to an increase in pressure of the gas discharge 
system. 

It was also noticed that as the magnetic field or pressure 
increased, the gas discharge requires more current to establish a 
normal glow (constant voltage region). This effect can be explained 


through the use of Townsend's equation, that is, 


= Ae BPE A (6.23) 


 |2 


His equation indicates that a, the number of electrons produced per 
unit length in the field direction, is reduced upon an increase in 
pressure. If a constant electric field is maintained, an increase in 


pressure will result in a larger number of low energy collisions, 
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Fig. 6.16 VI Characteristics of a Magnetically Stabilized Coaxial 
Discharge. 


Gas discharge contains 2:8:10 torr mixture 
of co, /N/He. 
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Fig. 6.17 VI Characteristics of a Magnetically Stabilized Coaxial 
Discharge. 


Gas discharge contains 4:16:20 torr mixture 
of CO,/N./He. 
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thereby reducing the number of collisions involving ionization. The 
presence of a magnetic field also promotes less ionizing collisions 
Since the charged particles must now gyrate about the magnetic field 
lines. During this process, the electrons and ions suffer periods of 
deceleration as they travel against the electric field. Thus, with 
less kinetic energy, the existing charged particles are less pro- 
ficient in stimulating additional ionization. The resulting loss of 
electrons changes the VI characteristics by reducing the current 
density. Consequently, more applied discharge current is required to 
maintain a normal glow behavior. 

This behavioral pattern was repeated for all gas mixtures 
analyzed including the 30:30:30 torr mixture of CO,/N./He shown in 
Fig. 6.18. Since the supply voltage was limited, the anode to cathode 
distance was reduced in order to achieve a breakdown voltage. Even 
at these very high gas concentrations, the magnetic field is highly 
effective in producing a normal glow discharge from an inherently 
unstable one. | 

Similar VI curves were made for the fluid ballasted, multi-pin 
anode structure and are displayed in Fig. 6.19. Examination of these 
curves also show a voltage variation with magnetic field. In addition, 
the VI characteristics of this electrode design show a glow-to-arc 
transition when the magnetic field is absent. This brief transition 
is speculated by the dashed curve. From these curves, it is also 
quite apparent that the constant voltage discharge is operating in 
the normal glow regime. This behavior is consistent with the smal] 


annular shaped cathode glow layer. Additional data were taken from 
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Fig. 6.18 VI Characteristics of a Magnetically Stabilized Coaxial 


Discharge.. Gas discharge contains 30:30:30 torr mixture 
of C0,,/N./He. 
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Fig. 6.19 VI Characteristics of a Magnetically Stabilized Fluid 


Ballasted Discharge. 
mixture of CO,/N,/He. 


Gas discharge contains a 2:8:10 torr 
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a similar gas composition but with a total pressure of 40 torr. 
Similar results were obtained but at a higher operating voltage. 

The upper experimental points in these figures represent a 
Specific discharge power loading of about 10 kW/1. This value is 
particularly encouraging in that it is substantially more than that 
which could be achieved under comparable conditions with a conven- 
tionally ballasted electrode structure. To date, an upper instability 
threshold has still not yet been observed, suggesting that even higher 


power loading is possible. 


6.5.4 Variation of the VI Characteristics with Magnetic Field 


It is well known that numerous changes in a gas discharge's 
characteristics can occur when it has been subjected to a magnetic 
field. Some of these erreces have already been observed in the form 
of a discharge rotation due to Lorentz forces, oblique angle bombard- 
ment of the cathode surface and a variation of the length of the 
cathode dark space. As demonstrated throughout this thesis, the 
applied Lorentz force has usually been found to be responsible for 
the plasma's distinct behavior. This section will now investigate 
the variation of the discharge current and voltage in the presence of 
a magnetic field. 

It has been previously determined that the normal glow behavior 
exhibited by the magnetically stabilized discharge, was due to the 
creation of a more efficient cathode emitter. Furthermore, the motion 
of charged particles is directly influenced by the combination of 


electric and magnetic fields. Thus, under the action of a magnetic 
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field that is transverse to the applied electric field, a significant 
number of electrons and ions will be deflected away from the discharge 
axis. In the transverse geometry, the radial magnetic field forces 
particles to travel in the radial and azimuthal directions. Similarly, 
in the coaxial geometry, an axial magnetic field forced the particles 
to flow radially and azimuthally. 

Through the use of Boltzmann's transport equation, it can be 
Shown that the magnetic field does not affect the component of charged 
particle drift velocity parallel to B. However, the drift velocity in 
the direction of #, perpendicular to B, decreases and a velocity com- 
ponent perpendicular to both # and B appears. The newly created drift 
velocity components are thus responsible for the dramatic change in 
gas discharge behavior. A more detailed analysis will be presented 
in Chapter 7. 

It is well known that a particle's drift velocity is related 
to its mobility. Thus, it becomes possible to associate the effect 
of a transverse magnetic field with the transport properties of the 
gas. In this respect, the mobility becomes a tensor and may be 
written in the form 


ee? 
nie (6.24) 


e +d 


Utes 
Ps) 

where Qs is the collision frequency tensor defined by Eqn. 4.39. 
Examination of the mobility component perpendicular to B, indicates 


that charged particle transport is now reduced in this direction since 
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In other words, the transverse magnetic field reduces the mobility of 
the charged particles in the direction of the applied electric field. 
Thus, the effective electric field normal to the electrodes, is 
decreased. This result confirms that the ratio Z/p is reduced upon 
application of a magnetic field. 

Due to the transverse magnetic field, deflected charged par- 
ticles can undergo several revolutions of the discharge active volume 
during an inter-electrode transit. Thus as a result of this process, 
a new electric field perpendicular to the original field will be 


created. The total effective electric field will then take the 


form 
ye 
F2 = f2 + KF ye (6.25) 
eff iy 4 a 2 


as derived by Golant [6.10]. In this expression, ries is the new 
component of electric field parallel to B and the second term is 
the reduced electric field normal to the electrodes. The effective 
increase in the overall electric field could then be detected as a 
rise in the measured discharge voltage. 

The VI characteristic curves in Figs. 6.16-6.18 illustrate that 
a Significant variation in discharge voltage does accompany a change 
in magnetic field. This association is depicted in Fig. 6.20 where 
the initial operating voltage is observed to rise linearly with 
applied magnetic field. As expected, Fig. 6.20 also shows that the 
discharge voltage is dependent upon pressure. 


Additional investigation has revealed that since the current 


density is a function of the electron drift velocity, it too will 
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Fig. 6.20 Variation of Discharge Voltage with Magnetic Field. 
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decline with applied magnetic field. This relationship has been con- 
firmed by Fig. 6.21 which shows the variation of discharge current for 
different values of B with different initial discharge currents. Over 
the region of interest, this dependence appears to be linear. However, 
due to saturation effects of the magnet core, the magnetic field values 
at high field currents tend to be slightly inaccurate. The plot in 
Fig. 6.21 is typical of those measured for each discharge geometry. 
In addition, the linear variation of discharge current was observed 
to be independent of gas species and mixture. 

Similar results were obtained for a small gas discharge tube 
by Sen and Gupta [6.11] where variations in discharge current were 
observed with respect to changes in pressure. However, their pressure- 
current relationship decreased as the pressure exceeded 150 mtorr. 
Consequently, this second variation was not observed in the present 
investigation since operating pressures normally exceeded 10 torr. 

ena tinther parameter of interest concerns the dependence of 
the plasma WEES » on the magnetic field strength. This 
relationship is displayed in Fig. 6.22, and is found to be initially 
linear for low magnetic field values. Overall, the graph indicates 
that the plasma impedance rises as the applied magnetic field 
increases. An examination of the conductivity tensor reveals that 
the perpendicular component o, is a decreasing function of magnetic 
field strength. Thus, current flow is reduced in this direction as 
previously seen. Consequently, this process will not favour the for- 
mation of highly conductive ionization channels or arc-like streamers 


along the discharge axis. 
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Fig. 6.21 Variation of Discharge Current with Magnetic Field. 
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CHAPTER 7 
DYNAMICS OF A MAGNETIZED GAS DISCHARGE 
7.1 INTRODUCTION 


Throughout this thesis, the influence of a magnetic field has 
been shown to have a substantial effect on the behavior of an ionized 
gas discharge. It has already been determined that the magnetic field 
can alter discharge characteristics through variations in the charged 
particle drift velocity, and collision frequency. Further analysis will 
show that the magnetic field is capable of modifying the expressions for 
the Townsend coefficients as well. All of the effects may be derived 
from an analysis of the electron distribution function. 

It is well known that a stabilized gas discharge is the first 
requirement for the construction of a high power laser. Thus far, 
it has been observed that the magnetic field is extremely effective 
in maintaining a stable, normal glow discharge. This chapter will 
elaborate upon the MGD stabilization scheme in both the space-charge 
sheaths and the plasma bulk regions. Furthermore, it will be 

explained how the reduction in the effective value of #/p leads to 
a less self-sustained discharge. 

The magnetic field has been found to be responsible for 
increasing the total number of inelastic electron-neutral collisions. 
This effect has the potential of increasing the excitation rate for 
vibrational pumping and thereby improving laser excitation efficiency. 
Verification of this theory will be revealed through a detailed gain 


analysis. 
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7.2 DISCHARGE STABILIZATION MECHANISM 


The initial study and preliminary analysis performance on the 
magnetic discharge structure of Fig. 3.2, suggested that the mechanism 
responsible for the stabilization of the glow discharge involved a 
highly sheared rotational velocity. The resulting velocity, driven by 
Spatially dependent Lorentz forces, developed an effective mixing of 
the ionized gas within the space-charge region near the cathode sur- 
face and in the positive column. 

It is well known that within a glow discharge, the largest 
electric field strengths are found within the cathode fall region. 
Here, field strengths can be several orders of magnitude larger than 
in the plasma bulk. This implies that the largest rotational mixing 
processes will occur in the space-charge sheaths directly below the 
cathode surface. Since the electric field strength within the 
eathoda fall region drops linearly with distance from the electrode 
surface [7.1], a strong rotational velocity shear will be developed 
in this region. This particular aspect has already been speculated 
by the oblique angle bombardment of positive ions upon the cathode 
surface. 

The analysis of the stabilization mechanism can be accomplished 
by examining the Lorentz forces acting upon the glow discharge. As 
indicated by Fig. 7.1, the discharge can be approximated as consisting 
of three radial and three axial regions. The radial regions are 
defined as zones of different magnetic field strengths and direction, 
that is, region I (B, smi) |B, | > B)> region II (|B | Se ), and 


region IIT (B, < 0, |B, > B Similarly, the three axial regions 
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Fig. 7.1 Lorentz Forces Acting Upon a Glow Discharge in a Transverse 
Electrode Geometry. 
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are defined by (1) an ion space-charge sheath found in the cathode 
dark space and negative glow, (2) an electron space-charge sheath 
located in the Faraday dark space, and (3) the neutral plasma bulk 
found in the positive column. 

The total current in the positive and negative sheath layers 
must include a conduction current J and a convection current e: result- 
ing from the azimuthal rotation of the space-charge layer. An 
examination of force terms in the © direction reveals that the J x B, 
force will cause a rotation of both the charged and the neutral fluids 
in the negative 6 direction. It is assumed that both fluids are 
strongly coupled together within this highly collisional gas discharge. 
The radial forces (7, xB.) acting on the charged fluid in regions I 
and III will tend to force the discharge current into an annulus 
(region II). This process has been confirmed by experimental 
observations. 

The convection current J, in the ion space-charge sheath, will 
create an outward radial force (Wee) which will tend to spread the 
ion sheath along the cathode surface. It is expected that this force 
will be minimal since the ion sheath velocity is low due to viscous 
drag. The actual location of the discharge within region II will be 
dependent upon a balance between centrifugal, radial pressure, and 
Lorentz force terms as illustrated by Eqn. 6.5. 

The strong spatial variation in the discharge current and the 
direction and magnitude of the imposed magentic field, implies that 
the rotational driving force (7, xB) will also exhibit a strong 


spatial dependence. This result confirms that the rotational velocity 
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within the charged sheath regions and in the positive column will be 
highly sheared. As previously discussed, the sheared flow profile, 
within the ion and electron sheaths, creates a uniform and controlled 
secondary emission from the cathode. The plasma bulk is therefore 
uniformly supplied with charge carriers from this electron emitting 
region. The strongly sheared flow within the remainder of the dis- 
charge volume will rapidly distribute any local nonuniformities in 
density, temperature, or current throughout the active volume. It 
is this convective process that is instrumental in preventing the 
development of self-regenerative electrothermal bulk instabilities 
which would normally lead to the formation of a constricted arc 


discharge [2.6]. 


7.3 ELECTRON DISTRIBUTION FUNCTION IN A MAGNETIC FIELD 


In order to thoroughly explain the variation of gas discharge 
parameters in the presence of a magnetic field, an analysis of the 
electron distribution function is necessary. The technique employed 


to solve the kinetic or Boltzmann transport equation given by 


(nf) + ¥-V (nf) + 2 (hte) -V (nf) = =-(nf) 9 C7a3) 


is based upon the assumed low anisotropy of the electron distribution 
function in a uniform electric field [6.10]. When the electric field 
strength is not excessive, the electron velocity distribution function 
depends only upon the velocity » and the angle 9 between the directions 


of velocity > and the field #. This assumption is quite valid for 
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values of # normally utilized in gas laser discharges. The dependence 
of the angle 6, which is due to the anisotropy of the distribution 
function, must also be small. In this limit, the distribution function 


can be expanded in orthogonal Legendre polynomials, that is, 
re foe) 
f(v) = ) f, (v)P,_ (cos 8) 


‘where the distribution function f(v) only depends upon the velocity. 
If the perturbation in the distribution function is small, it is 
possible to restrict the expansion of #(d) in terms of spherical 


harmonics to the first two terms in the series. Consequently 


flv) = f,(v) + cos(e)F, (v) (7.3) 


tt 


where the function f,(v) is the isotropic Maxwellian distribution 
and f(r) is the directed component of the distribution function. 


In the presence of electric and magnetic fields, the solution 


for the electron distribution function is sought in the general form 


URAC Oa eta) |= Deion Ce) 
p03) = f,(v) +2 — + b+ 


ase - (7.4) 


which depends upon all velocity components. Substituting the represent- 
ation for the altered distribution function, Eqn. 7.4, into the kinetic 


equation will produce a set of equations for the functions 7 ; i : 
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ene Ye and ar . The first equation is obtained by averaging the 
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kinetic equation over all the velocity directions by integrating over 
the solid angle dw. The equations for the directed components are 
obtained by averaging after multiplying by the direction cosines of 
the angles between the velocity vector and the coordinate axes [7.1]. 


This produces the equation 


+2 (nF) - 3s LF) - ee 
aE "alg? ae 3m v2 aye Ef, ) cee (7-5) 
and 
en af en 
5 ~ e> 1-0 Z Sa Ve ee 
ee) MC Rais) ay eed 1) oo va 


The collision terms 35 and S. are given by the terms 


tal ty 6(1,f) 
a| ae (7.7) 
and 
3 6(mf) 
Bs = Tale SE cos 6,,dw (7.8) 


where dw = sin edé6d¢ and ey is the angle between the & axis and the 
velocity vector. 
If the distribution function is independent of time and 


coordinates, Eqns. 7.5 and 7.6 may be simplified to 


-n 
ae dp? i f) 5 a (7.9) 
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The collision integral eh, can be obtained by substituting the elastic 
and inelastic collision integrals for the term Se]. The electron- 
neutral momentum transfer collision frequency is large compared with 
the electron-electron and electron-ion collision frequencies and thus 


Ss may be resolved into the form [6.10] 


> > 
Ce a vee! | (7 lily} 
Substituting this expression into Eqn. 7.10 produces the equation 
df 
eno 6 eos). 4 
m, ip * in? XE) ea Oe (e7eeliz:) 


If the magnetic field is directed along the z-axis, the pro- 


jections of Eqn. 7.12 on the coordinate axes take the form 


eBf 
Hin ang ie 
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where 9 is the collision frequency tensor described by Eqn. 4.39. 


These results indicate that the magnetic field can introduce anisotropy 


throughout the ionized gas in the form of a perturbation to the electron 
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distribution function. In the absence of a magnetic field, the 
collision frequency tensor becomes the unit tensor and Eqn. 7.13 
simply indicates the perturbation of the distribution function due 


to the imposed electric field, that is, 


Equation 7.9 can be conveniently expressed in the form of the 


divergence of a spherically symmetric flux in velocity space, that is, 


= = 2 
me Dune © oe ay? r) 
where 
2 p aes 
Poe ne Toe n 2 z Gaze 
E 3m, 1 3mé dv 


The term Z- (2°B) can be simplified by introducing an effective 


electric field, defined as 


Be ep = EF. (9-F) 


Vv e 
es 
If no electric field component exists in the y direction, then 


Eqn. 7.14 simplifies to 


Z 


Vv 

E2 = f2 + #2 | 
+2 

ef yr z CW Ne 


as derived by Golant et. al. [6.10]. By defining an effective electric 


field in this manner, Eqn. 7.9 can be written in a form common for both 


the presence and absence of a magnetic field. The dependence of 


ie on B appears only in terms of E, . Therefore, to determine the 


ies 


(7.14) 
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effect of a magnetic field on the distribution function, it is only 
required to replace £ by E, 


pica 


7.4 DRIFT VELOCITY IN A MAGNETIC FIELD 


As was earlier surmised, the magnetic field can dramatically 
change the drift velocity pattern for charged particles. This section 
will closely examine their altered behavior and derive an expression 
for the electron drift velocity in the presence of a magnetic field. 

The average electron drift velocity can be obtained by applying 


Eqn. 4.36 which states 


<vy> = | vf (Bao 4 


By substituting the approximation given by Eqn. 7.4 for f(v) into 
the above integral, the following result is obtained for the drift 
velocity 


1 UL 


. | 
ioe | fol) +] Le (v)|v2dv sin ededs. (7.15) 


cae e Q 


By inserting the direction cosines and integrating over the angles, 
the kth component of the drift velocity becomes 

“x ce k 
since all remaining terms contain odd functions which vanish upon 


integration. The substitution of Eqn. 7.13 into Eqn. 7.16 yields 


Vag z | vif, (v)dv (7.16) 
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Provided the electron collision frequency, ve) is a constant, the 


integration can be done by parts to give 
Ane 


ee ces 3 
tes 3n_m 2a? fo (v) 
ee 


- 3 0%, (v ao} 
0 0 


where the first term approaches zero at both limits when the Maxwellian 
distribution function is used for f, - Integration of the second term 
produces the result, . 


PEN e ..1'7) 


Equation 7.17 indicates that the electron drift velocity is also a 
function of the collision frequency tensor and accordingly the magnetic 
field. 

This result shows that the velocity component parallel to the 
magnetic field is unaffected. In addition the effect of the magnetic 
field on the perpendicular velocity components is substantial when 
Des ta kN Poh: In this limit the electrons are able to make numerous 
gyrations during an Tmaeecconineional time. As the electrons gyrate, 
they undergo an Bx B drift which increases the inter-electrode transit 
time and correspondingly the distance traveled. This process also 
increases the number of inelastic collisions, and as indicated by 
Eqn. 7.17, reduces the electron drift velocity in the Z, direction. 
Consequently, the rate of energy acquired by the electrons from the 
E, field component, is significantly decreased. 

A similar analysis for the ions can be performed using the 


perturbed ion distribution function in electric and magnetic fields. 
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The results will be identical in form, however, the collision frequency 
tensor will now contain the ion cyclotron and ion-neutral collision 


_ frequencies respectively. 


7.5 IONIZATION AND COLLISIONAL RATE COEFFICIENTS IN A MAGNETIC FIELD 


As previously seen, the application of a transverse magnetic 
field, complicates the analysis of an ionized gas. This complication 
extends to the determination of the ionization and collisional rate 
coefficients. In this section, the effect of a magnetic field on these 
coefficients will be examined to account for the noticeable change in 
gas discharge behavior. 

It was first theoretically determined by Wehrli [7.2], that 
the presence of a magnetic field was equivalent to an increase in the 
pressure of the .gas discharge. It was later established by Blevin and 
Haydon [7.3], that the discharge reacted as though the pressure had 


2.4, 
increased by a factor {1+ (w ) }*, when a magnetic field was 


ce es 
present. This meant that the charged particles behaved as though 
only an electric field was present and the pressure was increased 


fROM p< top! ,. thatitsi, 


el 
oe 


abe ae ae a (7.18) 


Here ia is the mean free time between electron and gas molecule 
collisions. 

The expression in Eqn. 7.18 has been experimentally ver- 
ified over the range 20 < #/p < 150 by Bernstein [7.4]. His work 


indicates that this formula is also applicable over a wide range of 
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Hall parameters, that is 0 < w s 25. Further extensions to 


ce es 
these limits may be possible once reliable collision cross sections 
are obtained. Nevertheless, these values indicate that the increased 
pressure concept is a valid approach to the study of most types of 
laser gas discharges under the influence of a magnetic field. 

In Section 7.3 it was shown that the presence of a magnetic 
field could be accounted for by defining an effective electric field. 
The expression for the new electric field includes the collision 
frequency tensor Q. The incorporation of this term is equivalent to 
an increase in the electron collision frequency in all directions 
perpendicular to the magnetic field. Hence, the corresponding 
increase in the number of electron-neutral collisions constitutes 
an effective rise in gas pressure. The results of this derivation 
are in agreement with the results of Blevin and Haydon [7.3]. 

The expression for the first Townsend ionization coefficient, 
a must now be modified to include this effective pressure increase. 


Therefore, the expression 


P 
becomes 
a(x) _ 2445 _|B 2545 
a Aqt Sy ate) } exp ney folate Oo te) } (7.19) 


where 4 and B are ionization constants. From the above equation, it 
can be shown that in the presence of a magnetic field, a(x) is 
reduced. This effect can be explained by considering the electron 
kinetics as they leave the cathode surface. Under the action of 


crossed electric and magnetic fields, a force is exerted on the 
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electrons which diverts their path from the anode. Thus, the electrons 
Spend more time between the electrodes. During the long inter-electrode 
transit time, electrons undergo more elastic and inelastic collisions. 
“Tf the applied Lorentz force is in the proper direction, the electrons 
will] gain less energy from the electric field and thus cause fewer 
ionizations per unit distance in the electric field direction. The 
reduction of a(x)/p with increasing magnetic field has been exper- 
imentally verified by Bernstein [7.4]. 

In the experiments undertaken, the crossed electric and magnetic 
fields create an azimuthal velocity, Vg - As the magnetic field is 
increased, more force is applied to the electrons and ions. This 
produces more rotation and the charged particles are able to spend 
more time between the electrodes. Furthermore, each electron or ion 
is involved in more collisions which transfer excitation energy to 
the neutral particles. Consequently, the charged particle temperatures 
P, and r. must also be lowered. This process is effective in reducing 
the charged particle temperature dependence on the production of 
instability modes. 

In the previous section, it was shown that the form of the 
electron distribution is altered in the presence of electric and 
magnetic fields. Through this variation, it was determined that the 
electric field component and the electron drift velocity were reduced 
in the direction perpendicular to Be Consequently, the magnetic 
field is also able to reduce the average electron energy, <&> and 


thereby its temperature, since 
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coo = | 6? e(ean)ae . (7.20) 
0 


From the previous discussions, the electron velocity is known to be 
vaécerident upon the electron energy and the ratio Z/n. Therefore, the 
electron distribution function can be written as a function of these 
variables. In this expression, the ratio Z/n is shown to be uniquely 
related to the average electron energy through the electron distribu- 
tion function. In addition, it has been shown by Nighan [2.2], that 
the electron distribution function can be extremely sensitive to 
variations in the parameter £/n. 

From the arguments presented, the ratio #/p will be decreased 


by the ratio 


E BE ] 
ELE Oayaly 
ie ania ase at sly 


Ww 45 
ce es 
for a gas discharge under the influence of a magnetic field. Should 
this change occur, Fig. 2.1 indicates that the gas discharge will 
become less self-sustained. Furthermore, as Z/n is lowered, the rate 
coefficients Ks are also reduced as shown in Fig. 2.1. This result 
can be mathematically derived since the collisional rate coefficients, 
responsible for instability growth, were functions of the ratio E/n 
through the relation 
1 
wid (22) 5 f 
a - (22 | Ef(EsE/n) + @, (Ede . 


J ne Ue 0 


where the electron energy distribution function now takes on the form 
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Consequently, the presence of a magnetic field helps to reduce the 


ability of the rate processes to initiate a micro-instability. 


7.6 GAIN CHARACTERISTICS OF A MAGNETIZED GAS DISCHARGE 
7.6.1 Temperature and Absorption Effects 


The values of vibrational population density and their distri- 
bution over the rotational levels are important quantities basic to 
the understanding of the mechanism responsible for the production of 
a population inversion in molecular gas lasers. Ina C0,/N./He gas 
system, the population inversion density and the vibrational tempera- 
tures can be obtained by analyzing the distribution in gain over the 
(00°1) + (10°0) vibrational transition which is known to exhibit 
stimulated emission. In this investigation, it is also important to 
determine the effect of a magnetic field upon these parameters and 
accordingly examine this behavior in terms of laser gain. 

Gain measurements were attempted with the oscillation-amplifier 
arrangement shown in Fig. 7.2. The oscillator consisted of a water 


cooled C0, laser with a 50 cm active discharge length. The amplifying 


2 
media was provided by the magnetically stabilized glow discharge. To 
increase the detectable signal, the effective gain length was tripled 
by folding the path of the probe laser beam through the gain media. 
This was accomplished by placing two gold coated mirrors within the 


vacuum chamber. The resulting gain length path was thereby increased 


(22) 
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to 60 cm. Salt windows were used to allow optical access into the 
chamber. 

The optical signal output was detected by a gold doped germanium 
crystal mounted in a vacuum vessel. In order to increase signal sen- 
Sitivity, the chamber was evacuated and cryogenically cooled with 


liquid N This procedure increased the Ge crystal impedance to 


2° 
approximately 1 M2 and with the proper bias circuitry, a voltage 
proportional to the detected optical signal could then be viewed on 

an oscilloscope. To avoid detection of the wide optical spectrum 
emitted by the glow discharge, a notch filter centered at 10.6 um was 
placed in front of the detector. This ensured that only stimulated 
emission would be detected. An increase in the signal to noise ratio 
was achieved by focusing the emitted radiation onto the detector with 
an anti-reflection coated Ge lens. 

The power gain was measured as the ratio of the detector out- 
put power with the discharge on, to that with the discharge off. How- 
ever, in all cases, the output signal was dramatically reduced upon 
operation of the gas discharge. It can therefore be concluded that 


the CO.,. laser beam is undergoing absorption as it passes through the 


C 
ionized gas. 

The absorption coefficient a, can be estimated from Beer's law 
which states 


I= ean (7.23) 


where a is the initial laser intensity and 2 is the gain path length. 
Under normal operating conditions, the measured absorption coefficient 


was a = 0.023 cm, for a discharge current of 4.0 A and a voltage of 
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660 V. The fact that only absorption was noted is attributed to the 


presence of a hot CO, gas. The value obtained for the absorption 


2 
coefficient compares favorably with those of Strilchuk and 

Offenberger [7.5], who measured similar values at an elevated gas 
temperature of 875 K. 

In order to determine the neutral gas temperature, an iron- 
constantan thermocouple was inserted into the glow discharge. The 
thermocouple was contained in a thermally conductive but electrically 
insulating ceramic tube in order to shield it from electrical and 
chemical reactions caused by the plasma. Figure 7.3 depicts the radial 
gas temperature profile of the fluid ballasted design operating at a 
current of 4.0 A. Examination of Fig. 7.3, reveals the significant 
effect created by the application of a magnetic field. For low values 
of magnetic field, the hot plasma was essentially confined to the 
center of the electrode. Near the electrode edge, the gas discharge 
was minimal and the temperature began to drop off rapidly. As the 
magnetic field increased, the recorded temperature at the electrode 
edge began to rise while the temperature in the central region 
dramatically decreased. This effect can be explained by considering 
that as B increases, the discharge moves radially outward from the 
center to form an annular shaped glow discharge. 

The results of Fig. 7.4 reveal the effect of the magnetic field 
upon three points located at the center (r = 0 cm), middle (r = 4.0 cm), 
and outside edge (r = 8.0 cm) of the glow discharge. The dramatic 
temperature change in the center of the discharge illustrates the 


rapid transition to an annular shaped discharge. The initially hot 
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Fig. 7.3 Radial Gas Temperature Profile. Discharge parameters include: 
J = 18 mA/cm2, 2:8:10 torr gas mixture of CO,/N_/He. 
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Fig. 7.4 Variation of Discharge Temperature with Magnetic Field. 
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gas in the middle of discharge is also observed to cool upon applica- 
tion of the magnetic field. Furthermore, the substantial rise in gas 
temperature along the electrode edge can be attributed to the convec- 
tive radial gas flow produced by the Lorentz force. From these 
results, it is obvious that the magnetic field is reducing the 
temperature of the ionized gas by convective motion. Through this 
technique, the ability of the discharge to generate local hot spots 
and thermal instabilities is being curtailed. 

Attempts were made to cool the gas temperature through the use 
of a heat exchanger and small electric fan. However the high power 
loading permitted by the magnetically stabilized glow discharge, 
overwhelmed the 500 W capacity of the existing heat exchanger. Since 
the input power loading of the gas discharge could easily exceed 
several kilowatts, the gas temperature continued to rise until an 
equilibrium was achieved. At steady-state, the average gas tempera- 
ture, as determined from Fig. 7.2, was approximately 550-600 C. 
Furthermore, it was not possible to install additional heat aehanders 
due to the small size of the vacuum chamber. In order to maintain a 
low gas temperature, the discharge current was limited to less than 
1 A. At this low current value, the discharge could not be uniformly 
excited nor the current density evenly distributed. 

It was initially believed that the natural radial flow of the 
gas discharge formed a closed loop convective cell and thus did not 
permit efficient use of the heat exchanger. This hypothesis was 
later confirmed by the computer simulation of the secondary flows as 


displayed in Fig. 6.3. In the next phase of this investigation, an 
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improved design will be constructed to allow the hot gas to recircu- 
late through a larger heat exchanger and up past the transparent anode 
with the assistance of a pair of tangential blowers. This new design 


is illustrated by Fig. 7.5. 


7.6.2 Excitation Processes in a Magnetic Field 


In order to ensure that the MGD stabilization technique is 
suitable for C0. laser applications, it is necessary to examine the 
effect of a magnetic field upon the gain of an excited media. However, 
due to the problems encountered with the removal of excess heat, no 
definite conclusions could be made. As an alternative, this section 
will speculate upon the possible beneficial aspects created by the 
presence of a magnetic field. 

The optical gain coefficient for the P and # branch rotational 
transitions between vibrational levels J and J + 1 ina C0., gas 


discharge, can be written as 


Moca! %s (7.24) 


In this expression, g(v) is the lineshaping function dependent on the 
dominant broadening process, and m is the index of refraction of the 
plasma. Upon examination of Eqn. 7.24, it is clearly seen that y(v) 


is primarily a function of ¢ and the population densities. It 


spont 
has already been shown that the magnetic field is responsible for 
lowering the temperature of the gas discharge through rapid convec- 


tion motion. It is hoped that this process will prevent the growth 
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of the lower laser level population density, thereby ensuring a large 
population inversion. This efficient technique of rapid gas convec- 
tion is widely established in the high power laser industry. However, 
their methods have relied upon the deployment of large high speed 
blowers to cool the gas. In this MGD technique, large velocities can 
be created by the Lorentz force which provide a non-mechanical means 
of cooling the gas discharge and hopefully improving the gain 
characteristics. 

The lifetime of the upper laser level t, 1s determine by mol- 
ecular collisions and is inversely proportional to the pressure [7.6]. 
Thus, as t, decreases with a rise in pressure, the pumping power per 
molecule will increase since it can be recycled more rapidly. It 
has been previously stated that the influence of a magnetic field on 
a gas discharge is equivalent to a rise in pressure due to the increased 
number of molecular collisions. In this manner, the magnetic field can 
effectively reduce the lifetime of the upper laser level and cor- 
respondingly increase the saturation intensity. | 

The analysis of Nighan [2.2] has determined that for an average 
electron energy of less than 1 eV, the electron-C0,, molecule energy 
exchange is clearly dominated by vibrational excitation, while above 
2 eV, electronic excitation dominates. More specifically, Boness 
and Schulz [7.7] have determined that the maximum excitation of the 
(00°1) asymmetric vibration mode occurs for an electron energy of 
0.9 eV. At this point, the collision cross sections of the (10°0) 
symmetric and (01°0) bending modes are an order of magnitude lower, 


and are thereby negligible. Thus, the population inversion by electron 
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impact and therefore stimulated emission at 10.6 um is enhanced for 
this low energy region. For electron energies greater than 3 eV, 
laser action is likely to be suppressed since the population inversion 
will be reduced by the diminished cross section of the (00°1) level, 
and by the enhanced cross section for the lower laser level (10°0). 
Consequently, by choosing the proper Z/p ratio, corresponding to an 
electron energy of approximately 1 eV, it is possible to optimize the 
lasing process. 

Throughout this chapter it has been shown that the presence of 
a magnetic field can alter the electron distribution function, the 
electron energy, and consequently the ratio Z/p. Thus, by proper 
adjustment of the magnetic field strength, #/p can be sufficiently 
Tray edi such that the energy transferred to the upper laser level can 
be maximized. With the peer ce ih of a larger vacuum chamber, a 
complete verification of the anticipated effects of introducing a 


magnetic field will be investigated. 
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CHAPTER 8 


CONCLUDING REMARKS 


8.1 SUMMARY OF DISCHARGE STABILIZATION PROCESSES 


The presence of a transverse magnetic field has dramatically 
demonstrated its stabilizing influence on a glow discharge. It has 
been experimentally and theoretically shown that a magnetic field can 
effectively suppress both plasma micro-instabilities and bulk thermal 
instabilities. This highly effective stabilization technique has, thus 
far, achieved a specific power loading of over 10 kW/1 without exper- 
iencing a glow-to-arc transition and without the benefit of external 
gas cooling. These results are especially encouraging since an upper 
instability threshold has not yet been observed, suggesting that even 
higher power loading is possible. 

One of the reasons for this exceptional degree of stability can 
be attributed to the trapping of charged particles into Larmor orbits. 
As a result of this orbital motion, the charged particles spend more 
time travelling between the electrodes, and are therefore involved in 
more inelastic collisions. The increase in the electron-neutral and 
jon-neutral collision frequencies essentially acts as though an 
increase in gas discharge pressure had occurred. Through this 
mechanism, the application of a magnetic field lowers the parameter 
E/p such that the glow discharge becomes less self-sustained. Since 
E/p is a measure of the average energy gained between collisions, the 


charged particle's drift velocity and its energy are correspondingly 
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reduced. As a result of a decrease in electron temperature or Z/p, the 
magnitudes of the dominating ionization and attachment rate coefficients 
are also reduced. Since the rate of these collisional processes is 
decreased in the presence of a magnetic field, the probability of an 
attachment induced ionization instability is dramatically lowered. 
Consequently, charge production and loss in the plasma are effectively 
decoupled from the applied electric field. 

It has been determined that the Lorentz forces, acting upon the 
Space-charge sheaths, are responsible for maintaining the gas discharge 
as a normal glow. Furthermore, the magnetic field is suspected of 
controlling the electron emission from the cathode by increasing second- 
ary emission yield through oblique ion bombardment of the cathode. 
Together, these processes dynamically inhibit the growth of electrode 
induced instabilities. 

The Lorentz forces are also responsible for stabilizing the 
bulk thermal instabilities through a rapid rotation of the charged and 
neutral fluids. The high rotational velocity exhibited by the glow 
discharge can rapidly convect nonuniformities in electron density, 
temperature, or current and distribute them throughout the discharge 
volume in a time less than the normal instability formation time. 

Since the applied magnetic fields are spatially dependent, large 
Shear forces exist which are able to break up streamers and arc-like 
filaments. 

Due to the imposed magnetic field, pressure gradient forces are 
also created which drive secondary radial and axial flows. A one-fluid 


MGD computer model has predicted the existence of a convective cell 
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formed by these secondary flows. The discovery and prediction of large 
secondary flows is an important advance in resolving the behavior of 

a glow discharge subject to a transverse magnetic field. Furthermore, 
the secondary flows play an important role in the stabilization process 
by allowing the convection of hot gas out of the active discharge 
volume. This radial flow system is highly efficient since the hot 
neutral gas has only one-half the discharge residence time of con- 


ventional transverse flow designs. 


8.2 SUMMARY OF COMPUTATIONAL WORK 


The development of the computer program MAGIC has produced a 
detailed simulation of a time dependent, partially ionized gas discharge 
hes an imposed magnetic field. Through the use of a one-fluid MGD 
model, this program can accurately predict the magnitude of the highly 
sheared rotational flow and the existence of strong secondary flows. 
These rapid flows provide an efficient non-mechanical means of 
reducing the gas residence time, thereby convecting electrothermal 
instabilities from the active discharge volume. 

The use of a large time step, iterative ADI finite difference 
technique has enabled MAGIC to operate with time steps larger than the 
CEenlimice. suensa technique can be regarded as state of the art in 
terms of advanced computer codes. The utilization of variable time 
Step size controls has proven to be an excellent method of maintaining 
maximum computational efficiency while avoiding the limitation in 
time step size that is imposed by explicit techniques. This effort, 


though indeed a major one, must be considered an initial one, 
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Suggesting many new ideas and areas for future consideration. Thus 
far, MAGIC has proven to be a valuable research tool and hopefully 
upon further refinement, a variety of discharge devices and structures 
may be entirely designed and simulated using this powerful computa- 


tional aid. 
8.3 FUTURE CONSIDERATIONS 


In the near future, the computer code will be extended to 
include the conservation of energy equation. This additional feature 
will allow simulation of heat transport throughout the gas discharge. 

A more ambitious project is currently planned which will include the 
development of a two-fluid MGD model designed to determine the complex 
interactions formed between the charged and the neutral fluids. 

It is also planned to extend MAGIC's capabilities to handle a 
variety of discharge configurations. This would initially include the 
coaxial discharge geometry discussed earlier and a new axial discharge 
device. In this geometry, the Lorentz driving force is supplied by an 
axial electric field and a rapidly rotating radial magnetic field. The 
' time dependent nature of the code will also be utilized to ascertain if 
the magnetic stabilization technique can be adapted to pulsed TEA laser 
systems where the duration of the discharge is very brief. 

To further substantiate the observed and numerical results, 
additional diagnostic measurements will be undertaken. Of primary 
importance is the determination of the electron density. This could 
be accomplished by the deployment of an electrostatic probe at various 
locations within the plasma. Alternatively, an externally generated 


microwave source canbe used as a non-perturbing probe. In this 
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case, the electron density is simply determined by varying the signal 
source frequency until the phenomenon of cutoff is reached. The 
rotational velocity of the ions and neutrals can be determined by 
measuring the Doppler shift of the spectral lines emitted by the 
plasma particles. The ion and neutral gas temperatures can be cal- 
culated from the Doppler width of the emitted spectral lines. 

Due to the broad scope of the investigations described in this 


thesis, a thorough study of all aspects of this project was difficult 


to accomplish. However, a significant amount of data and understanding 


now exists which confirm some basic theories. It has now been con- 
firmed tnat magnetic fields can shape, confine, and stabilize glow 
discharges suitable for use in high power gas lasers. With the 
promising results shown thus far, the MGD stabilization technique 
should have a major impact on the design of future gas lasers. It 

is anticipated that such a laser will be built and tested in the near 


future. 


228 


oe Re 

‘ a nt : i 
i. — in 

“i ie 7 


apie stile gnbere’ vi meieny rave NES , 
aT Bat ita e: ‘siaitag: 0 reno ot ; th | | 

vd banterietah od ea ‘eh waduse! one anor ‘ott , wat | 
eto batt tie dant raced ete itd. v0! pnmarts 

“fs » on ve aoiutoveynis? 2h Teena nas not oat iat 
s2entt Tevdzege boat tina ‘ene te Asetw s9Taqod i: oe 

pry nt heder ray ondtteg traavat wie 4g sgooe: badd att ‘as: ieee 
‘ " sey bts oew 2o0% on eins 26 espayes (ié ts bua fipuerods Bs 
grinned mesons bie dah to tome sonar tinea 5 reveal: ve 
Od Hewd we RENE sivod atewd aoe rN: Honig 


wote ahh ides > ‘tar ean? OME TAS weet 
eas ahh aeget pete yh. at ow oe x y 40 
Sit Anaad Hi antes Miiiade ae ott ee aus oe 

$1 .syaeet> p swt % aeewy wi ne sane 
yean of? “i bictzant ‘tina: shies ad irtw “t 


ye iW 
a 
2 a 
: : fl 
ee ft 
i ih 
x 1a te a 
‘ is Ae 
\ on ie Te - 
a on ; \- / 7 


a By. pe i : > | i 
bn an ie othe 
-_ : aw W 
; tae to o] ; oo ’ 
iy eke : vt ’ ‘ ‘i _ 


. aid 7 ; — An ue = ip Cts 
srs cece ch aij onal ne if 
| Ee i a 
| 1. ee ay Heures - ane 
; i — . r ig iy a mah * ye! 


REFERENCES 
CHAPTER 1 


1.1 W.d. Wittermann, "High Output and Long Lifetimes of Sealed-Off C0, 
Casers= <7 Appr apnyse Lect... VOU 1 Vs. Now. 10: psi3375 1967. 


1.2 A.J. Demaria, "Review of CW High-Power C0. Lasers"; Proc. TEBE. 
Vol. Gis INO. 26:5) po 735. 1973. 


1.3. W.L. Nighan, "Progress in High Pressure Electric Lasers", 
presented at the 11th International Conference on Phenomena 
in Ionized Gases, Prague, Czechoslovakia, Sept. 1971. 


1.4 A.C. Eckbreth and J.W. Davis, "RF Augmentation in CO, Closed- 
Cycle DC Electric Discharge Convection Lasers", Appl. Phys. 
ettaeVolaaciey Nowsiigsp.. 125, 92: 


1.5 C.A. Fenstremacher, M.J. Nutter, W.T. Leyland, and K. Boyer, 
"Electron-Beam Controlled Electric Discharges as a Method 
of Pumping Large Volumes of CO, Laser Media at High Pressure", 
Ropisw Rhys Sheela VO Wer cUsaNOn 2.00.00. 4 ole 


1.6 H.d.d. Seguin, J. Tulip, "Photoinitiated and Photosustained Laser", 
Appl. Phys. Lett. Vobk. 21; No. 9, p. 41451972. 


1.7. A. Jowan and J.S. Levine, "The Feasibility of Producing Laser 
Plasmas Via Photoionization", IEEE J. of Quantum Electronics, 
Voill.; OB=8,, No. BM, ip 6827, 1972. 


1.8 Q.P. Judd, “An Efficient Electrical C05 Laser Using Preionization 
by Ultraviolet Radiation’, App): :Phys. Lett...,Vol. 22, No. 3, 
DieI56@oro: 


1.9 0.P. Judd and J.Y. Wada, "Investigations of a UV Preionized 
Electric Discharge and C09 Laser", IEEE J. Quantum 
Electronics, Vol. QE=10, No. 15 p. 12, 1974. 


1.10 A.E. Hill, “Uniform Electrical Excitation of Large Volume High- 
Pressure Near Sonic C09-N9-He Flowstream", Appl. Phys. Lett. 
Vols 18s NO St pa IA e197: 1:. 


1.11 G.A. Garosi, G. Bekefi, M. Schulz, "Response of a Weakly Ionized 
Plasma to Turbulent Gas Flow", Phys. of Fluids, Vol. 13, 
No Shi 0927950 97 02 


1.12 W.C. Roman and T.W. Myers, "Experimental Investigations of an 


Electric Arc in Transverse Aerodynamic and Magnetic Fields", 
ATAA Journal, Vol. 5, No. 11, ps 2011,. 1967. 


ae) 


ents Obes 7s Ye 2atat 3 slike pred ba ine 
PBT SEE. cE ON ET LOY eating a8 


1341 .oond , Yeheedd o anes nae WO ¥9 une . 


a ol af aie 


| “razed sthteeti spies: thot ate we ig! | 
pnamoiet? noe Sone metro Fenot tented) qtr oi rh ; ste pio 
ee ae «ty of sortase Supa mae ‘bastgat at a ita 


 ailgise@ TD 302 at nats scammpa we at ah so 36 ters fies a re ue 7 
aye 108 5 . pan Lit aah govAe) aorarnd be 1339 5a 7 i 
| ayet eee fh af os 


a weeat hea: ,Deletged Tal ieatihel om er ‘tens 
horlist < 28 caguedoalt vile peflfevtnc2 | ‘sata 

, aruegs fG iptH 45 eT! tM ee te eter so gatgeet 8 
OT 08 a oft. OS FOF 5.2594 


,‘esael te i Teds weosti one sien ater ¥ ae 
A 2 chet a ot «AB Thiele. as AT in 


Avant ani: abord Pe: ig tteaiaa ott". 
stag tae wierd ng hid als ty Beat wot ty 


¥ 
ow y 
es 


nottas iniot ary patel aed ae bape 
ef On, e .foW .. 30a)". i 


ata ined ah isa oben 
bins nent) a tea) oO) 
vate ST sa iS 08 Oh 


: eit a ol %0 tat 3 x3 ions rss 
iv Oe. i ho rel) ith rahe eae 
py rue «8, al 81 id 


4 freies a, pa Pe ai 7 
ioae teeer se rete 7 a ie oe y 
eT eae 


‘bate fansl 6 10 aenagaen 
“et ¥ oY cebbolt 20 OUT 4 


at3angam brs. atiysnybor 


: os te enotsaptteeuns inal : 
, eb: oe ‘Bi 
ae Wet, Pe 


230 


1.13 Y.Y. Winograd and J.F. Klein, “Electric Arc Stabilization in 
Crossed Convective and Magnetic Fields", AIAA Journal, 
Vol S72, Nowr9, p.. 1699, 1969. 


1.14 C.J. Buczek, R.J. Freiberg, P.P. Chenausky, and R.J. Wayne, 
“Magnetic Stabilization of the Plasma Column in Flowing 
Molecular tasers,” 9\Proc., TEER, Vol. 59, No. 4, p. 659, 1971. 


Tous: BA. Tozer, © Rotating Plasma’ seProce IEE,» Vole VIZ ,°Now 1, 
De285. 1965. 


1.16 B. Lehnert, “Rotating Plasmas", Nuclear Fusion, Vol. 11, p. 485, 
1971. 


1.17 B.W. James and S.W. Simpson, "Isotope Separation in the Plasma 
Centrifuge", Plasma Physics, Vol. 18, p. 289, 1976. 


1.18 B.W. James and S.W. Simpson, "The Viscous Dynamics of a Rotating 
Phasmary,, (PlasmayPhysicss Wolly 20% ipie/59.,. 1978 


1.19 H.E. Wilhelm and S.H. Hong, "Plasma Rotation by Electric and 
Magnetic Fields in a Discharge Cylinder", J. Appl. Phys., 
Voslex: 485n.NO.. t2sac pe h6 ler cliO7 7. 


1.20 S.H. Hong and H.E. Wilhelm, "Boundary Value Problem for a 
Counter-rotating Electrical Discharge in an Axial Magnetic 
Fqe@ld. A) Appl aPnyS VOL. 49.: NOs Is Di £46, 21978, 


1.21 M.M.B. Wijnakker, E.H.A. Granneman, and J. Kistemaker, "A Study 
of a Weakly Ionized Rotating Plasma", Z. Naturforsch, 
Voln sdasep 672701079: 


CHAPTER 2 


2.1 R.A. Haas, "Plasma Stability of Electric Discharges in Molecular 
Gases", Phys. Review A, Vol. 8, No. 2, p. 1017, 1973. 


2.2 W.L. Nighan, "Electron Energy Distributions and Collision Rates 
in Electrically Excited No, CO, and C02", Phys. Review A, ~ 
Vols 22, NG. Oa Dango 9s) 11,9702 


2.3 W.J. Wiegand and W.L. Nighan, "Plasma Chemistry of CO-No-He 
Discharges  ssApOweEnyS. | Lette. .vob. ccs No. i; 
Den OS we 197 34 


2.4 W.L. Nighan and W.J. Wiegand, "Influence of Negative-Ion 
Processes on Steady-State Properties and Striations in 
Molecular Gas Discharges", Phys. Review A, Vol. 10, No. 3, 
px 922, O44. 


ne notiast Teded? iy a}vaae Li a 7 wig 
et tsrrnt BALA a “er ot tans 


, Sega b A bees Jateinal 4.4. re 
entwott ay peered TD mets ant TO not: 
NGF 288 2g ob a Tot Tait .30%8 


Gh ST! toy ‘aa. 304 Pomests anitate 


at i 


Et att POM Naren somal Teamen!9 eatera® } 


- 


Sh Tied" 

ames ’S oth ab cotsenhad panes sidan’ My? od : 
ARE SBS eT  Fow seh sme T9 ary “i ‘a 
) F: “= oe Pa ee. 
ond pagnn a ty eat mist eogsety eer” jhwanat? Wid. ng’ 28 8! 
aa "PVG? Ret uy 28 ae eaoheye anes t? ee 
ang oF ysoetd ad pebted to em £4" sunol oat . f ina low 
eo ends. . Tag ct Paitin agredoetd tnt zhistto 
_ RE 188 qed omigt v 
6 “or mot de A autee ena 1s" sales He Be gna iil IS hy 4 
2 doncem tates ca ob aptitgere sgh ty oatistor ‘G. aire. 
BOE RE gol a CAs ll gate 4 Nga ale ein ' 1 a 

| en} ae ri 


ybor? A" oRRNOENT,. ity Debt: yoamanans hut Cea 
shoe iaren Ss Pasi ent Ls 


ysiuoshom at zegnarnh) shag 
BOE | LIDS ae * ay 


29t6A cohen ited bas: gnottud hates 0 prerd co ia ‘eget si a 
_ gt wae ma i gh pas me % te390.13 | TLE oo Spee |e 
aes " « +8 a ao toy, ee 


‘ brepatW sbi 4 i 
Aha a see ct 


. . 


a 


Hing ng a 2) 
7’ re 


aa 


, ; o 
vv 


we ‘s. ih 
: 


y hs *., 
7 if 7 
mo iy 
SD TH. cee gg of > 
wre ina 


yl VE je, Pe ee 
7 , ' 7 


i: ‘ties tts oe | 


e om Sait Taw a La a os 


2.5 W.L. Nighan, W.J. Wiegand, R.A. Haas, “Ionization Instability in 


200) 1M 


CHAPTER 


or ici. 


SO ee 


Ss | Re 


cia Sm oe 


oro | fF 


S20) “Gr. 


CHAPTER 


Gl FBS: 
4.2 G. 


fee ey 


4.4 N.A. 


Ho. (Ge 


4.6 S. 


CO>eLaser Discharges’ Appl .. Physsy Lett: ,7Vol.. 22,°No. TI, 
Dabo) Ful os oF 


.L. Nighan and W.J. Wiegand, "Causes of Arcing in CW C09 


Convection Laser Discharges", Appl. Phys. Lett., 
Vol-<S255-No;: 145: pia6s3 .fd974" 


Holland, W. Steckelmacher, and J. Yarwood, “Vacuum Manual", 
Spon itd, 1974. 


Nasser, “Fundamentals of Gaseous Ionization and Plasma 
Electronics", Wiley, New York, 1971. 


- Nam, Ph.D. Thesis, University of Alberta, Dept. of Electrical 


Engineering, 1978. 


. Thomas, "Software for Numerical Mathematics", Loughborough 
Leic. Eng., Academic Press, p. 315, 1974. 


-A. Frungel, "High Speed Pulse Technology", Vol. 1, Academic 
Press, New York, 1965. 


. Glasoe and J.V. Lebacqz, "Pulse Generators", McGraw-Hill, 
New York, 1948. 


Tanenbaum, “Plasma Physics", McGraw-Hill, New York, 1967. 


. Sutton and A. Sherman, "Engineering Magnetohydrodynamics", 
McGraw-Hill, New York, 1965. 


. Landau and E.M. Lifshitz, “Fluid Mechanics", Pergamon 
Press, London, 1959. 


Krall and A.W. Trivelpiece, "Principles of Plasma Physics", 
McGraw-Hill, New York, 1973. 


. Capjack, D.M. Antoniuk, and H.J.J. Seguin, "Dynamics of 
Magnetically Stabilized Laser Discharge", J. Appl. Phys., 
Vode oe ONO wer De 40175) 1OGk. 


Chapman and T.G. Cowling, "Mathematical Theory of Non-Uni form 


Gases", Cambridge University Press, London, 1970, 3rd edition. 


231 


wt spitidaaiel rohaveaaat saat 
“ghT OR ae POE yt eed 7 


“ 


OOS WI wh entontt’ on sntanth ubeuge tw “ | rn st ! 
2 TRS RN a. ‘aponedioe rt. 79864, MOTT: 
PART bar * a _ | 


"Tagagh woesel” , Deore AY bee per an emen: ww “stant. a 
| BOT - bid oR? 

3 ei | pa pesy Le Ay a 

gia Ft ‘ie HOF tus biol auosene rv alstneaabaud" a 

» TYet Tahal wont et 0 eo 


| ae P 
*o  dqut ,atiedtA te ‘poesia ty 0.08 ih 


fears oai4 
treont a 
Tek ee oe 
yuo Warinuad 5, "aptt sonar seart inal a ar6us9 02", oes a 
7 RT DHE ge dant, 2 awe, ent. | 
TRAE 4A HON “op toa nate’ et gti: aoe Ms 
| | / i aa Bn at "gh 3 ; 
ET a “eer erning® sate sso Wb bys so2et. 
a on ia . 


habe bal 
cf vied 


Stas te wis +. sata 


Any ow 


*o' eo tiage .” beng ‘Laie 
Basil a ea 


ri rte eo i” 


A : 
4a ot wf 
ae we a fe 6! 4 CET - 


en) 


7 ee St Be 


Meda, Pe 


CHAPTER 


Se OR: 


Be2 30): 


aor “DE 


2 a a 


5.09. &. 


5.07 WwW, 


SEW Rm Lc 


ee eles 


CHAPTER 


Gl ash. 


Gua cit. 


.O. Hirschfelder, C.F. Curtis, and R.B. Bird, "Molecular 


Theory of Gases and Liquids", Wiley, New York, 1954. 


J. Elliot, 0.P. Judd, A.M. Lockett, and S.D. Rockwood, 
"Electron Transport Coefficients and Vibrational Excitation 
Rates for Electrically Excited C05 Gas Lasers", Report 
#LA-5562-MS, Los Alamos Scientific Laboratory, 1974. 


M. Morse and H. Feshback, "Methods of Theoretical Physics", 
McGraw-Hill, New York, 1953. 


H. Pennington, "Introductory Computer Methods and Numerical 
Analysis", Macmillan, New. York, 1965. 


W. Peaceman and H.H. Rachford Jr., "The Numerical Solution 
of Parabolic and Elliptic Differential Equations", J. Soc. 
IndusteeApp |]? Mathes Wo lef3s }pst2854 19552 


Potter, "Computational Physics", Wiley, New York, 1973. 


F. Chen, "Introduction to Plasma Physics", Plenum Press, 
New York, 1974. 


H. Finan, Ph.D. Thesis, University of California, Davis, 
Department of Applied Science, 1980. 


R. Briley and H. McDonald, "Solution of the Multidimensional 
Compressible Navier-Stokes Equations by a Generalized 
implicit Method” mJ. So§ Comps Phys. Vol: 24, p. 372,, 1977. 


D. Richtmyer and K.W. Morton, "Difference Methods for Initial 
Value Problems", Interscience, New York, 1967, 2nd edition. 


R. Lindemuth, Ph.D.Thesis, University of California, Davis, 
Department of Applied Science, 1971. 


S. Craxton and R.L. McCrory, "Two-Dimensional Calculations of 
Non-Spherical Laser Fusion Implosion", Report #99, 
Laboratory for Laser Energetics, University of Rochester, 
1980. 


J.d. Seguin, C.E. Capjack, D.M. Antoniuk, V.A. Seguin, "A 
Magnetically Stabilized Radial Discharge for a High-Powered 
Waser }ersppl ganvs i eeett.., Vol...39, No. 35. p. 203, 1981. 


coe 


seetuslaleiade aroha kha 
Jha ret work. Poe phage tanta, 
boowtoa? 0.2 bs ‘\ttesend Ma Com 2.0 

tient toe? Fares dared W Sine. zanetottsageay 


prdes ."ovenel ap) 609 bartoxd yf fisokres 
NCR woreda. | at tt gas toe aome tA 201 ‘i 


tee heal sineadT Yo sbovtem ,tosddesd Wer 
Caer oY, wail tinea 2 


— ar wars 
aS al Fe soaghew oa 
) mt *% : ala) 


ll bins power ies eitioationsal® tort 199 7 
CRRT ea welt “a FO Si lena 


fie ni Sule? | ga) yoru eet” . : ois bot 
t aint SUeS to itnavaetrd ots 
eet tg | at ‘te 


oki i 


atone bietatitiet: to ton 8 
ae esas sail gak 


ie rhe -g co "a VF Re: 
COAT SKE. RS co. Ne 


fatstal yo? ehodte gan 
AED this has ORT tit 


eh atest we | 


i Nh , { F be, if : ee i walt ae P "Se 
‘ ‘ T ‘ i h ey + fer , 4 a em a ’ 
we ai 7 a U eur td i r A rn ve ea le , 


tw not sehupisd Taco] aan 0-o4 


gO 
contents to \ {gral ‘ 


233 


6.3 C.E. Capjack, H.J.J. Seguin, D.M. Antoniuk, V.A. Seguin, 
"A Magnetically Stabilized Coaxial Laser Discharge", 
Apple Pnys= bs.2V0l. 26, p. 161, 1981. 


6.4 J.-L. Crandall and A.W. Cooper, "Acoustic Generation of Moving 
Striations in Neon", Proceedings of the 8th International 
Conference of the Phenomena of Ionized Gases, Springer, 
Vienna, p. 133, 1967. 


6.5 L. Pekarek and V. Krejci, "The Physical Nature of the Production 
of Moving Striations in a DC Discharge Plasma", Czech. J. 
PhySi.VB.. VOlwoN In. De ico, 196). 


6.6 A.K. Brewer and J.W. Westhaver, "The Cathode Region in the Glow 
DVSChardesws 70 sae CNYS., VOlMeG. ID. . 779.4 1937. 


6.7. M. Cloutier, "Computed Physical and Thermodynamic Properties 
of Various He/CO2/No Mixtures", DREV Technical Notes, 
Quebec, Canada, 1970. 


6.8 M.L.E. Oliphant, "Liberation of Electrons from Metal Surfaces 
by-Positive tons ; Proc. Royal’ soc. , A=127, pe 3735. 1930. 


6.9 0. Almen and G. Bruce, "Collection and Sputtering Experiments 
with Noble Gas Ions", Nucl. Instr. and Methods, Vol. 11, 
Diy 2ovig) bo ks 


6.10 V.E. Golant, A.P. Zhilinsky, I.E. Sakharov, "Fundamentals of 
Plasma Physics", Wiley, New York, 1977. 


6.11 S.N. Sen and R.N. Gupta, "Variations of Discharge Current in 
a Transverse Magnetic Field in a Glow Discharge", J. Appl. 
PHYSS Dia VOdeeo Deol coy k, 


CHAPTER 7 


7.1 T.N. An, E. Marode, P.C. Johnson, "Monte Carlo Simulation of 
Electrons within the Cathode Fall of a Glow Discharge in 
He lium’; dus PNYSeLU,- VON lO. porcoldseloy ee 


7.2 M. Wehrli, "“Funkenpotentiale im Transversalen Magnetfelde", 
Ann. der Physik. Vol.: 69, Sers 4, pb; 285, 1922. 


7.3. H.A. Blevin and S.C. Haydon, "The Townsend Ionization Coef- 
ficients in Crossed Electric and Magnetic Fields", Aust. 
UeePAVSSS VOleesiae NOs TO, Dd. Joo. 


7.4 M.J. Bernstein, "Townsend Ionization Coefficient for Hydrogen 
in a Transverse Strong Magnetic Field", Phys. Rev., 
VOIR ler Demodes 902. 


een AM ‘estima net acca os 
omretag?G wasd Melead® bast hidekey ih 
BRT fot ct 88 OV 6 


“GA EVOPE ho ntsauanle ‘oKiauoon" eso) ih tne ioe ney. ine 7 
fendtsenvotnt nab edt to ‘apntbeanor’ ‘nash at; ane an ba. 
vapid? .asees bak twol to enemonor® ant bit: sipgerdose Hi : 
; ay a : 
notpowioyd ont TH: ouial taokeyd ont” “Patent oe V bag ones ‘ 
ob .. does . ane? agpedoeld 20.6-af ‘enor dstys2: pnivom, t 
. - faRk , 28h ae i: FOV 2.8 Lay | 
i S ~ . a RD Oe m : _ 
HOfg af? at norpen abotte) oft” ‘ae. We Lb bas ‘qewod aA 
NERF RAED «8 ov, RYRT se Fagan hy “wren? |e 
a Sak A y 
aT syaaats jaws naboto et “Ons Teste bosudind" -wohu01) a 
soto IndhnitgaT ¥IA0 ea Beat cabo auo te¥ to. 
Arca pie | We 
gant facem Wer it anys teeta: to fof revedt" < ghadgt tar’ a. J ) , 8s 
pee re 2 Vee. one toys 904 4 “enol. mete ea ' 


{ a. 


zsnemiveqad parasatuge brs eottaet too” soured . sta rent, 2. ; 
ET ote sebodism bee ean, hi enol 2a 61606 { oe < 


io awaits” -vOnBRR: Be: ro ¢ es 4k 4 
tet aot wan . taxi) 


gy” wi Wai ba on 
Adal pics oo 


nt peeve? ooxeiaetl te ‘aa £ 


a 


Firat “Ly “ep reroate wore b “ tenga sevav a | 

a ) ‘aera = ee es a a a if 
Wh rh : a ? 
| Viens 


i 
iy 1% 


gu9 ‘oor, 3B Met ‘a ¢ ae 


to eo ofGo aantok? tt | RE oR 
vr baie a & to: 7 hel, Short, 7 oad ntdtiw enongoal3! bal, va 
Ver WRB nt if FOK: 4a sett sia "mute Ha i 


7 ke gras — ts gremarmiet <Ftot aM an 


| ss Patt cihewg wl 


234 


7.5 A.R. Strilchuk and A.A. Offenberger, "High Temperature Absorption 
ih COp-at NUSGe une nPop. Optics, Vols 13.0No. oll; 
p. 2643, 1974. 


7.6 A. Yariv, “Quantum Electronics", Wiley, New York, 1977, 2nd edition. 
7.7 M.d.W. Boness and G.J. Schulz, "Vibrational Excitation of C05 by 


Electron Impact", Phys. Rev. Letters, Vol. 21, p. 1031, 
1968. 


nob droeda sear i Hott 
u “* nie A, ; 


nol tbe ons. ree Avy ‘a At saaecoats 


d 609 0 ooraegtaia fanobrevare® extuded. 
OGr.g SF. enatte]. VR aii 


Re 
ffs 
’ 4 
“ i . 
See orl io © 
~~ 
- / : 
i) 
os 
“a 
. 
t 
oa 
” 
1 
1 
a ty 
cc * pee Sted 
me ar) i a: , - 
a ; Oo —— ae mn : ve “ 
’ 5 : A j 
ay Dail @. Nee OP an, CACO onl ee ye 
- i 0 ean ’ “YT 


- 7 ; % 4) ' ae A i ve a F. 5 ' + he f 
i , ‘ i fj : 6b “i : : 
_ a ¢ 1 An ee ay a K, x aah r) 


APPENDIX A 


In order to utilize the conductivity tensor co in a cylindrical 
coordinate system, it is necessary to perform a similarity transforma- 
tion. This new conductivity tensor would then be valid for any 
orientation of the nonlinear magnetic field depicted in Fig. 3.7. 

The first step in this coordinate transformation requires a change in 
basis vector from local magnetic field (Cartesian) coordinates to 


local cylindrical coordinates. This can be represented by 
x = Rg! (A .1) 


ee a ° ° . . 
where x = } “Ze. is a vector representing cylindrical coordinates and 
L 


z' = we, is a vector representing the local Cartesian coordinates 
- 

along the magnetic field. The matrix R is an orthogonal trans forma- 

tion matrix formed by the change of basis vectors from Cartesian to 

cylindrical coordinates. From Fig. A.1 the new cylindrical basis 


vectors are written as 
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By using Eqns. A.]l and A.2, the transformation matrix becomes 
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Fig. A.1 Transformation of Coordinate System Basis Vectors. 
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where & is an orthogonal matrix such that R =R =R. 


In local Cartesian coordinates, Ohm's law may be written as 


> 


J's gte«EF , (A.4) 


Now, using the transformation Eqn. A.1 for both the current density 


and electric field vectors yields 


or 
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where 5 = Rott} is the similarity transformation. 
The Cartesian conductivity tensor o' is normally aligned with 


the 2 component of the magnetic field. Written in matrix form, o! 


becomes 
FF, ee 
oh hee ce | (A.6) 
0 0 oa 


Performing the necessary matrix multiplications as indicated by 
Eqn. A.5, the similarity transformation produces the new cylindircal 


conductivity tensor 


nae 2 = i -o )si co 
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